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0.1 Lecture 1 


0.1. 


1 Frenet Triad 


Parametric Curves are functions which are defined in terms of useful 


and related quantities, often in téme or arc length. Before going on 


much further, it is necessary to define a few expressions: 


The velocity (or tangent vector) can be thought of as the change 


in position with respect to arc length, or T' = Ae = fi(s). W 


also want to define the unit tangent vector such that T = va 
where its magnitude is always equal to 1. 

yon the tangent vector we can define the normal vector as 
iF a or in other terms, T’ = kN. k here is called the curvature; 
to be specific it is the extrinsic curvature. 

The Binormal vector is defined as B = T x N. The derivative 
of the this vector can be related to the normal vector N by 
B' =—7TN using a quantity called torsion, represented by T. 


Since N = B x T, 
Neo B x PLB xP =—1(N xP) 4i(Bx NM) =1B a 


Together, these quantities form the Frenet Formulas: 


T! 0 k O T 
eS | Se Oe N 
B' 0 —7r 0 B 


Example: Helix of Unit Velocity 


Let us look at the function 
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where a,b and ¢ are arbitrary constants and s is the length of the 
curve itself. It could equivalently be written by its components 


x(s) = acos (=) , Yy(s) = asin (=) and 2z(s)= a 


Cc 


Plotted, this curve gives us a circular spiral that rises or falls de- 
pending on the ratio of b to c. Note that if b = 0 then we have 
a two-dimensional curve (f’(s)) that traces out a circle on the xy- 
plane. The tangent vector can be found to be 


= A a b 7 
s\= (-=) sin (=) be (=) cos (=) 3 + (*) k. 
C G C C C 
By analyzing the tangent vector, we can guess the following orthog- 
onal vector in the (x, y) plane: 


N(s) = (<) COS (=) i+ (<) sin (=) j+0k 
a € a c 
(observe that we chose |N| = =, From this, we can find the deriva- 
a 


tive to be 


N"(s) = (--) sin (=) i+ (=) cos (=) 5+ 0k . 


With the tangent and normal vectors, we can also find the binormal 
and its derivative: 


B(s)=Tx N= (-2) sin (=)i+ (*) cos (=)i- 1k , and 
B(s) = (-=) cos (=) i= (<) sin (=)i+ Ok . 
ca C ca C 


Comparing with the equation for N : 


C2 C2 az+bh?- 


If b = 0, then there is no torsion (7 = 0). 
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Observe that if T is a unit vector, then a? + b? = ¢. 


From the equations for T” and N, 


The parametric function of s can be Taylor expanded with a 
differential As. This expansion immediately reveals the quantities 
we defined above: 


As)? As)3 

F(s) = F(0) + £70) s+ "0 A + pry) CP 4. 
—— es 
To ko. No No’ko 


and, if we only look at a small portion of the curve, we can assume 
that k = constant. Then, 


(As) | ko(™Bo — koT)) 9) ae 


f(s) = f(0) + ToAs + koNo 5 6 


In the last term we neglect a second order one (k§7j), and the 
final form of f(s) for small increments is: 


(As)? 


f(s) = f(0) + ToAs + koNo + koTBo oa (0.1.1) 


Oscular Plane and Geodesic 


The helix is an interesting case that can be used to introduce the 
concepts of the oscular plane and the geodesic. 

In the case discussed above, the oscular plane is the one formed 
by the tangent T and the normal N. The vector B is perpendicular 
to this plane. In this plane, we may draw a circle that is at the 
point in consideration to the tangent T, and that also contains the 
differential segment of the helix that encloses the point. Thus, the 
circle has the same curvature R as the helix (R = k). 
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Helix 


Next, we will show that if the helix is the curve that joins two 
points A and B on the surface of a cylinder, then it is the shortest 
distance between the two points, just like a straight line is the short- 
est distance between two points in a plane. The curve that fulfills 
these properties is called the geodesic. 

Assume a point moving on a cylinder from A to B. We will show 
that if the distance between A and B as measured by the length of 
the curve is the shortest possible, then it is a helix. If we take s as 
the parameter on the curve, then the coordinates of a point on the 


curve are 
x = acos(s) 
y = asin(s) 
£=% 
or 


dz = (—asin s)ds 
dy = (acoss)ds 
dz=dz. 


The function z(s) is to be determined. The differential length dl is 
such that 


(dl)? = (dx)* + (dy)? + (dz)? = (ads)? + (dz)? . 
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Thus, we want to minimize the following integral: 
B B B 
f= | dl = | a2(ds)? + (dz)? = 1, ds a2 + 2” 
A A 


d 
where 2’ = zee The problem of minimizing an integral belongs 


to the gibieet of Variational Calculus. This branch of calculus is 
discussed in Lecture 7: Review of Analytic Dynamics. It may be 
presented in the following way: if the function f(s, y(s), y’(s)) is to 
be integrated between the point A and B while s is a parameter (it 
could be, for instance, length or time), then the integral 


r= [° f(s.u(s).u(s)) ds 


will have a minimum or maximum when y(s) fulfills the following 
differential equation (named after Euler, the mathematician who 
found the solution in 1700): 


of d_af 


a) 
Oy as oy ) 
In our case, 
a? + 2” 
d 
and z is the function y and thus 2’ = y/ = 7 
S 
The Euler equation says 
O f(s) _ 
Os Oz! 
or 
d € 27) ) 7 
ds\2 Ja? + 27) — 
Thus, 
2! 
ee 
where c is any constant. And so, 
se ure _ p2 


Sipe 
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or 
z=Ds. 


where D is a constant; so our curve with the minimum length be- 
tween two points is a helix. Since D could be any constant, there 
are an infinite number of helices. 

The fact that the helix is a geodesic is easily seen if the cylinder 
is unfolded and converted into a plane with width 27r. The helical 
coils become straight lines, and obviously the shortest path between 
two points A and B is a segment of a coil. We note that on the 
cylinder, as mentioned before, we may create an infinite number of 
geodesics, each characterized by a slope a. The slope of the helix 
passing through two points A and B is predetermined by those two 
points. 

Another interesting example that helps to visualize the meanings 
of the oscular plane and the geodesic is the case of a thread on a 
sphere; the only force acting on the thread is the normal reaction 
from the sphere. The normal reaction passes through the center of 
the sphere at any point on the thread. So, the plane containing 
the normal is the oscular plane, which cuts the sphere into a grand 
circle. The grand circle through the equator and prime meridian is 
the geodesic of the sphere. If the arc of that circle is greater than 
180 degrees, then the distance between two point A and B on this 
arc is not the shortest distance, because the shortest distance is the 
supplement arc. So, in this case we have a geodesic which is not 
necessarily the shortest distance between two points. 
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0.1.2 Connecting Functions 


As in the case of Frenet frames, when the frame moves, the deriva- 
tives of T’ : N, and B give the intrinsic properties of the curve — its 
curvature (/) and its torsion (7). These derivatives with respect to 
the curve parameter are called “covariant derivatives.” 

It’s interesting to note that those derivatives (T", N’, and BY) can 
be expressed as functions of the frame itself (I, N, and B). We can 
think of these spatially as the variation from some point P to some 
point P + oP. 

We can define this variation a bit more formally if we start in 
some general frame: €}, €,...,€, such that €;-€; = g,. In this 
frame, we can define the variation of the vector dé; as 


dé; = Wi; ‘ej 4 


Here, W;,; is a function that depends on the point P in our gener- 
alized coordinate space (qi, q2,---;Qn) Where the frame €, @),..., En 
is located. It is closely related to the matrix form for Frenet’s for- 
mulas previously shown: in the simple case of the Frenet frame, the 
derivatives of T' ; N , and B with respect to the parameter s which 
varies along the tangent T are the covariant derivatives. As such, 
a=7, aN, Sas, and 


dé; = Wi5e; . 


The matrix W,; is 


T’ 0 k 0 
Wi, « | N’ —k 0 
B' 0 —r O 


Considering the special case of an orthogonal frame (€; - €} = 6;;, 
where 6;; is the Kronecker-Delta Function), we can use the product 
rule to expand that relation to (dé;)é; + €;(dé;) = 0. Using this, we 
can determine 


Wi; = = Wy 5 and Wir =0. 
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W,; is a connection function which depends on point P and it 
expresses the differential variation of e€; in the direction of e;. If 
we act under the assumption that W;,; is a linear combination of 
coordinate differentials dq,,dq,...,dq;,...,dq; then it can be re- 
written as 


Wi = dvi) = i -dq* , 


where the coefficient re was first introduced by (and is named after) 


E.B. Christoffel in 1869. 
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0.2 Lecture 2 


0.2.1 Metric Tensor 


A point M is determined in the Euclidian space by its coordinates 
with respect to a reference point 0. The vector OM in rectilinear 
coordinates is fixed by the n coordinates M(x, £2,...,2n)- 

The idea is to find a set of curvilinear coordinates with center 
in M that allows us to describe the geometry of the space around 
M. Let us call y', y?,...,y” these coordinates. The n y* coordinates 
are called “normal” or “natural” coordinates at M. At M + 5M this 


system may be different. 


As M moves by 6M it generates a set of n vectors @1, €2,..., €n, 
such that 
_, OM 
C4 at 
Oy* 


When introducing a metric g;; in the space, the “surface” where 
M exists becomes a “geometric surface”. The differential length 
ds? = (6M) is given by 

ds? = gigdy' dy’ , 


where g;; now may depend on the coordinates y's. 
If the curve described by M is parametrized with parameter tf, 
then the length between two points a and 6 is 


_ b | dy’ dy! 
iM = f Fi at at dt , 


and the volume dv formed by @jdy' with origin at M is 


dV Saf lg) dyjaus,dy™ 
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Example 1 


The curvilinear coordinates at M are the spherical one with the 
vectors @, along the OM direction, @, on the parallel around the 
x® axis and @9 in the meridian passing by M. Thus, if 7, and @ 
are the unit vectors in these three orthogonal directions, then 
ee=7 C= rsin(O) Cg=rb. 
The matrix element is then associated with the differential length 


ds? =dM x dM by its Pythagoric form 
ds? = (dr)? + (rsin0)?(dw)? + r*(d6)? . 


So 
p w 6 
r fil 0 0 
M= 7] 0 (rsin@)? 0 
A \0 0 r? 


The change in variables from the rectilinear coordinates to the 


curvilinear ones and viceversa are: 


(21, 2X2, £3) —_ (y', ie y°) 


(a? + 23 
r= yajta34+ 3 w= tan! (=) 6é=tan! oe 
x 3 


1 
(yu, 9°) —> (@1, #2, £3) 
x =r sin(0) cos(w) y = rsin(@) sin(w) a= 7 cos(G) 


(Where the y’ are understood to be, in this context, r,0,~). 


0.2.2 Euclidian Space 


Assume a Euclidian Space with coordinates p', p?,...,p". A point 
M(p',...,p”) in the space moves in the direction of p' generating a 
vector 

+, OM 
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Thus the differential of IZ equal dM is a vector 
dM = @idp' . 


(The summation on the repeated indices is understood.) 
If the motion of M is on a curve, the element ds of the curve is 
given by 


@s=dM -dM = @;- @;dp'dp’ . 


The scalar product is called g;;, the metric tensor of the space. 
That is, 
ei: € 


Tip = Ca gs 


Simple Examples: Vector as Tensor 


Let us assume a plane in the Euclidian space. The axes are given 
by the unit vectors € , and @». A vector V has components 2! and 
2 
z*. So 


ad > > 
V Soe) +2 €x.. 


The components x! and 2? are the contravariant components 
of V. The projection of V = OM on @; is V- €; = x; and the 
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projection of @» is ae €2 =X. These are called the covariant 


components of V. 


1 2 1 2 
1=0guit2' ge =x +cos(a)z* , 


where a is the angle between the @; and @2. This equality may 
be verified geometrically by observing that the two shaded areas are 
equal. 

Let us show in a simple form that the scalar product of two 
vectors V and W may be written as 


V-W=ViW, or V-W=ViWigi; . 


In fact, 


W=V-@ , W=V-e , VV et+V ee 


and similarly for W. Thus, 


V1W, + V2We = VW - 21) + V2(W - So) 
= viw'e, + We.) = ey + cae 
=VIW'gn t+ VW gis + V°V" gar + V?W G00 - 


In a matrix form, the transformation from contravariant to co- 
variant of the vector OM of the figure could be written as : 


ti\ (gu ge\ (z'\  f 1. cosa\ [a' 
to}  \gor goo) \e?} \cosa 1 i 


The inverse of the matrix (g;;) is (g‘7) such that 
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aaa ame a (a 
g aos 1 


Then 
i_ ij, 
C= G°xj, 


which is a simple case of raising indices by means of the fundamen- 

tal tensor. In the Euclidean space where €;- @; = €j- @i, gij is 

symmetric. Observe that the determinants of g’? and g;; are equal. 
Let us see the value of g;; in other curvilinear coordinates. 


1. Polar coordinates 


The point M is given by 


(p,0) : ds? = dp” + p*dé? 
@y=pd and €3;=0 
where 6 and p are the unit vectors. So 


G00 = pO - pO = p Joop = P*P=1 
Jop = Goo = 9 


0.2. LECTURE 2 21 


2. Spherical Coordinates. 
Then @9 = 6, €, =p, @y = psindy 


where 0, p and w are the orthogonal unit vectors. 


Joo = C9: C9 =P Iop = 1 Jus = p” sin? 0 


1 O 0 
g=|0 p? 0 
0 0 p?sin?@ 


0.2.3. Tensors 


Let us change the coordinates such that the space will be described 
by a new set of n coordinates q!,...q”. The coordinates p’ of M are 
given by the following functions: 


where the index 7 runs from 1 to n, and the transformation is 
biunique and regular. That is, the Jacobian of the transformation 
is different from zero (see note at the end of this section). Let 
us call Q a system of several functions that we call “elements” of 
Q: Ey, Ey,...,E,. Assume that EF is just a function F(p',...p”) 
that changes to Q(p'(q',...,q”),-..p"(q',...q")) in the change of 
variables. 

Obviously 


So a simple substitution keeps the “element” E invariant. 
Now assume that the element is the derivative of a function. So 


_ OF (p', p*,-.- 4B") 
Opt 


E,(p', p’, eee Pp") 
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By changing variables E;(p',...,p") —> Q;(q',...q”) such that 
j 


OFDM G hd eh O lO ne xO") (2") 


5 = Op" Og 


we find that a simple substitution of the variables will not preserve 
the element. It will give the first factor only. The second factor that 
is the ratio of the old variables p',...,p” to the new ones, q!,...,q", 
is necessary to preserve the character of the element, the geometric 
meaning of the derivative. 

Elements like Q; are called covariant elements. We identify it by 
a subscript and they transform by having a factor which is the ratio 
of the old coordinates to the new ones, symbolically written as Be 

The metric g;; is also a covariant element. Since there are two 
indexes 2 and 7 in the transformation there will be two ratio-factors. 


In fact, 


ae On' Op! 
ds” = g,dp'dp! = a 7 dq’ = p,sdq"dq? , 


Prs = Gij qr ag F 


So gj; is a covariant with two indices (or of order 2). Observe that 


where 


despite the fact that g;; may be a function of the coordinates, our 
derivation is correct since we just expressed dp’ and dp’ as a function 
of the q'’s. Let us assume now that the element is a factor A” 
which forms a linear combination of the derivatives of a function of 
(p',...,p"). That is, the linear combination: 


OF Diep") 


A’ 
Op” 


which, after the change of coordinates f(p!,...,p”) > h(q',...,q”) 


Ar OF (pty P") = ORG, 9") _ psOF(t,---P") Op" 
Op” Og? Op” Og? 
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So the new factor B® is related to A” by 


Oq* 


Be = A’ : 
Op” 


The element A” is said to be a contravariant tensor. It changes 
to B* by the ratio of the new coordinates to the old ones. In gen- 
eral, if P is the element with n number of covariant indices, the 
corresponding Q after the transformation p — q is written as 


Op*! Op*” 
Sn Oqr aie aqrn : 


Pry soot = Vey 


And for the contravariant system 


wg” Oq’” 
po Bom 


PTeotn = Qe 
Every index r; or s; takes the n values of the coordinates: 


ere) Or - +p -s.%00 


in the second side of the equation. Repeated indexes imply ad- 
dition. 
0.2.4 Operations with Elements 
Addition 


Two covariant or contravariant elements of the same order may be 
summed to result in another element of the same type, with the 
same order, 


Sr tn = Aby,.itn a Bry tn : 


This means that (it can be easily proved) as the coordinates 
change p —> q, 
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Op*! Ope 


Sra (P) ae Oe agri wees Og’ , 


These elements may be also multiplied with a result of another 
element whose order is the sum of the order of the factors 


7, x Bayes = te ee eee : 


For contravariant elements we have similar relationships. 

It also can be shown that when we multiply elements with co- 
variant and contravariant indices the result is an element with less 
indices, as in the next example: 


S1,...8 a 
A ; 9 5 aera pe Tip) sok $ 


Changing from Covariant to Contravariant and Vice- Versa 


From the previous equation that expresses the contraction or satura- 
tion of indices we may see that the following equation is an identity: 


go = IG" gi ; 
We observe that by definition 
gv = a 
l9| 


where C“” is the cofactor of g’”. Thus 


l9| TE jue = Gw9"” |g| 


and the following must therefore be true: 


Iwo tll Ue 
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Another way of writing this is 


where 67; = 1 if a = yp and is zero otherwise, and OS 41 it COSY 
and is zero otherwise. 

In general, we can raise or lower indices, i.e. go from contravari- 
ant to covariant indices and viceversa, by using the fundamental 
tensor g;;. For instance, 


TLyeeln — AT181 Tn8n 
R™ = wl Oe Ties ova 


For a contravariant element, 


= S1,...8 
Fae = Grist ++: Qrgiht ars 


This and previous sections have been based on the works of M. 
Christoffel (1859) and that of M.M.G. Ricci and T. Levi-Civita 
(1898, Methodes de Calcul Differential Absolu et leurs Application). 
What we have called an “element” here is in today’s language a 
“tensor.” 


0.2.5 Tensor Density 


In a change of reference from p — q, é; > €7,, and giyj — fiy; then 
(ds)? = gidp'dp’ = frsdq"dq? , 


where 


feces _ Op! Op! 
rs — Yij dq qs : 


We put this into a matrix form by introducing 
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thus, 


Taking the determinant of the above matrix equation, 


feel = lar lle2l gig] = F719 
where J is the Jacobian of the transformation. Thus, 
lfrs] = Jy/loigl - 


In the four dimensional space of Minkowsky, 


The element under coordinate transformation is related to the 


light path: 


(ds)* = e*(dt)* — (dx)* — (dy)” — (dz)’ . 


The fundamental matrix in this space is then 


10 0 0 
_fo -1 0 0 
aaa ea: es 
00 0 -1 


The element differential “volume” in 4-D space with respect to the 
old variables is 


dv = dp' dp? dp dp* . 
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In the new variables, the differential “volume” (which is different) is 
du’ = dq' dq? dq? dq* 
and the equation that relates these two is the Jacobian, 

dv = Jdv'. 


Thus, if we are working with the integral of a tensor T' in the 
space of p’’s, {Tdv may be directly related to the space of g’s if 
instead of T we use the tensor t, such that ,/|g|t = 7. In fact, 


[tev =f Viorel ta 
= | it 
= f Vital Goa! 
= f /\fsl tae’. 


Obviously the tensor TJ’ has changed to T”’, and so t to ?’ after the 
change of variables. 


0.2.6 Notes on the Jacobian 


In changing n of the coordinates from p — q, the variation of the 
point M changes from dp',...,dp” to dq',...,dq”. Let us call F 
the transformation which relates both coordinate systems, 


Pag tay) 
In order to ensure the inverse relationship and its uniqueness, the 
Jacobian must be different from zero. 


1 1 1 
ee as 
Po, Ea tag: 
J= q qd q 
iT n nr 
Pies oe 
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The partial derivative of F’ with respect to q/ is indicated by F';. 
The differential change of the coordinates will be related by 


dp! | ie nay ce OO we ele 
dp” Ee, at Ley NAG 
Let us see this in two dimensional spaces 
p=F\(qjq) and pe=F*(q',¢’). (0.2.1) 


The functions F! and F? create an affine transformation or a 
mapping which assigns every point (q',q’) in the region R% to a 
point (p', p*) in the region R?. 


RP Ra 


In matrix form this is 


dp" — {411 G12 dq’ 
bea 7 2 4) 
where the matrix (a,;) expresses the linear transformation shown in 
(0.2.1) . Thus, R? is the image of R? through the matrix (a;;). In 
order to assume that the inverse transformation exists and is unique, 
the determinant of the matrix (a;;) should be non-zero. 
Thus, in the affine transformation a line is transformed into a 


line and a surface into another surface. 
In a differential form, 
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where 


Thus, the inverse transformation exists and is unique if 
(7) 40. 


A classical application of the Jacobian is done in the calculation 
of multiple integrals when we change the variables in the integrand. 
Specifically, the integral defined in the space R? is to be transformed 
to the region R%. For simplicity, let us assume two dimensional 
regions, 


=i F(p',p”)dp' dp? . 
Rp 


Next we show that in the R? region, 


= Ve F(p'(q',¢@),p(¢',¢)) J dq’ dq? 


where J is the Jacobian of the coordinate transformation. 
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Rd 


The element of area of R? is (hk). The corresponding element of 
area in R?% is twice the shaded triangle. We recall that the area of 
the triangle of coordinates (0,0), (a,b), (c,d) is $(ad — bc). In the 
above figure the area of the shaded triangle is 


1 (Op! Op? 7 Op! Op? bh 
2\0q! 0q? dq? Oq! 


The area of the parallelogram is twice this, so the ratio between 


both differential areas is just the quantity inside the parentheses, 
and this is the Jacobian of the transformation, 


dp' dp? a J dq'dq? . 


Differential Volume 


1 2 
(pl + kp, p? + kos) 


RP 


As asimple generalization of the previous equations, the four-dimensional] amy 


volume is 
dv = dp'dp*dp*dp’ . 
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In the new variables, it is 
du’ = dq'dq’dq?dq* 
and the equation that relates the two is the Jacobian, 
dv = Jdv’ . 


We have also seen that the metric tensor in both spaces is related 


by 
V | tral =J | 93 : 


Then, if we are working with the integral of a tensor T’ in the space 
of p’s, [ T'dvu may be directly related to the space of q’s by using the 
tensor t instead of T, such that ,/|g|t = 7. In fact, f Tdv is 


/ | gig lt dv =f lfsl5 av =f VlfalS dv’ = [ lifesleedue 


So, the tensor T has changed to 7’, and t to ¢’ after the change of 
coordinates. d 

If the p coordinates are chosen to be orthogonal, then their metric 
tensor g;; has a determinant equal to 1. so, the metric tensor f/f; of 
the curvilinear coordinate q’s has a determinant such that 


Vfrsl = J 


where J is the Jacobian that expresses the transformation from p’s 
to q’s. Then, the element of volumen built by the hyper-parallepiped 
€; dq‘, é5 dq?,... is given by 


dv = 4/|frs| dg’ dq?... . 
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0.3 Lecture 3 


0.3.1 Absolute Tensor Derivative 


Most differential equations in physics contain first or second deriva- 
tives of some tensors. We must be sure that the derivative does not 
change the nature of the tensor. In other words, the derivative must 
also be a tensor. For instance, a vector, when moved differentially 
to a close point, must still be a vector. Thus if a vector V(q!,...) 
is written as 
V=evi 

where é; is the unit vector, and v' its contravariant component, its 
total derivative will depend on the point where it exists. When it 
moves to the next point, both v' and &; may vary. Using the letter 
D for the ”total” or ”absolute differential”, 


DV = déV' + dVié; 
where V is a function of g‘,q?,... etc and also &;. The differential 
variation of €;,dé;, may be written in a linear differential form 


de; = Tyg’ é; 

dé; _ 

dgk — 

The connecting factor I’, is called the Christoffel symbol. 
Thus, 


jo 


DV = dVie, + VTE dg*e; . 
Changing 7 for 7 in the first term, 
DV = (AVI + ViTY dq" )e; . 
In general, a vector W is expressed as W = W'eé;;, where W* are 


its contravariant components. Thus, the contravariant component 
of DV is the expression within the parentheses. We indicate it as 


DV? =dVi + VTi dg? . 
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Since 
- OV 
dV) = ——dq* 
Ogk ; 
then . 
ov! : 
DV! = (Fy + VT aa" 


So the absolute derivative of V expressed in its contravariant 
components is given by: 


As we will see next, the absolute derivative of V expressed by its 
contravariant components V; along the coordinate q” is: 


To prove the above, we assume that 


OV; 
DV; = age + BEV; 


where Bi. is to be determined. 

Let us take a scalar made up from the contravariant component 
of W, W*, and the covariant component of V,V;. That is, the scalar 
S is defined as 


SaW Ve. 


Since S is a scalar, the absolute derivative is equal to its func- 
tional derivative. Thus, 


V,+W' (0.3.1) 


On the other hand, taking the absolute derivative of each com- 
ponent, 
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D,(W'‘V,) = (D,W*)V, + W?(D;V,) . 


Introducing [ and B in the last equation, 


D,(W'V;) = Vi0,.W? + ViIT;W? + W'O,Vi + W' BEV; . (0.3.2) 


Comparing (0.3.1) and (0.3.2), we see that the following sum 
must be 0: 
VTi W? + W' BL; =O: 


Changing the repeated indices, one finds that 
W'V; (Ti, + Bi) =0. 


That is, Bj, = =P 
Thus, the absolute derivative of the vector V in its covariant form 
is: 


DV; = OnV; —THVi - (0.3.3) 
Absolute Derivative of the Metric Tensor, g;; 


Gig = Gi Cj 


The variation of the unit vector é€; is given by the Christoffel symbol 
as: 
dé; =Tme7 og" . 


Thus, the variation of gj; is: 
dgis = (Tyida* ei) - & + Tida*em) - & , 


and since g;; is a scalar, 


Agi; 
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From these last two equations we find: 


ts = Tigi + Pig Gim - 
We could rewrite this equation by using the other form of Christoffel 
symbols: 
Veit =Tyj and [yi gim = Vrij - 
Thus, 
we = Deje + Vrag - (0.3.4) 


Christoffel Symbols as Functions of the Derivative of the 
Metric Tensor 


Written ‘i 

dgk as GJij,k 
we rewrite three equations like (0.3.4) with different orders of the 
indices. 2, k, 7: 

Gijk = Veig + Dige 


Gni,j = Ujre + Pray 


Gina = Tage + pe 
Then, 
Gaet ous — bie = Tepe = 20 ae 
Thus, 


1 
Die = Peay = 5 (Yiask OGG = GiRa) 


We recall the symmetric properties of g and I’, as well as their 


definitions: 
Gig = Dj 
dé; = T,e:dq” 
Dijx = Pegi 


TF 91 =i; 
Vigng’ = Vig 
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Absolute Derivative of a Tensor of the Second Order (Two 
Indices), Tensor ¢;;: 
_ Oi; 
ijk = Bqk 
This equation may be derived by using the expression of the 
absolute derivative of a vector W; such that W; = tigV? , where V4 
is an arbitrary vector in its contravariant form and uniform in the 


space. By uniform, we mean 


and so its absolute derivative is given by 
Dey al 
(The proof is shown in the next note.) 


The equation above could be applied to the second order metric 
tensor g;;. Thus, 


O9ij an O49; 
Digij = Ogk —VRGng — D591 = Ogk =D lees 
According to the equation derived earlier, 
Ogi; 
Age = Digi + ray 
and 
_ 99% = = 
Digiz = aq* — Vegi — Prag = Pages + Prag — Tags — Thay = 0 - 


This conclusion is known as the Ricci Theorem or Ricci Identity. 


Note on the Absolute Derivative of a Second Order Tensor 
bij: 
Let us multiply it by the uniform vector V’. Since uniform means 


that the vector has the same component in every point in space, i.e., 
ov) 
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DVS Sl ev 


Now the vector tizV? = W; has an absolute derivative 


On the other hand, the absolute derivative of a product follows 
the same rules as the contravariant one 


Dz (tiiV”) = V! Dati = tijD.V? = Vi Dati + ihe ‘ 


Equating these two previous expressions for D;(t;;V7) we get 


OW, ; 
age Dig Wi = Vi Dytig + tigVimV™ - 
Since . 
OW; _ OV? = Ot; vi 
Oq® Og* Og" 
and 


s 
Wi =TieV* , 
we rewrite the previous equation as 


Ot; 
Ogk 


Vi -Tt.V? = VDaty + tT V™ - 


Now changing the index s for 7 in the left hand side of the equa- 
tion, and also 7 for m in the last term of the right hand side, we 
obtain 


Obs 


Dg (tig) = Oak 


— (UfRtmag —Tfitim) - 
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Summary of the Absolute Tensor Derivative of the Previous 
Examples 


1. The scalar tensor is a simple function F (q',...,q”"). Its ten- 
sor derivative is just the common derivative that we learn in 
Calculus. Calling Q the tensor and D;@ the tensor derivative 
with respect to the variable q*: 


OF tay, .090") 
Og* ; 


Q=F and D,F= 


2. The next order of tensor is a vector V. We say that the order is 
equal to the number of subindices. If we called this number m, 
then the scalar corresponds to m = 0. For this case, m = 1 and 
Q will have only one index. When it is up (down) it corresponds 
to the covariant (contravariant) components of the vector. In 
fact, we have seen that 

OV; ov? 


pie T3,V; and D,V’=—+T% 


DiVi = ag* j 


The element or tensor in question is Q = V, Q; = V; and 
Qi=V?. 


3. Finally, a second order element @ was discussed, the tensor ¢;; 
of two covariant indices. We found that 


Ole ie 
Og* 


In general if Q is a covariant tensor of order m, that is, if Q is 
expressed as Q,,....r,,, then its absolute tensor derivative is 


aQ rst 
_ Tighe 


PO in: Sars gk - ¥ ee oe Ome ene ee ee . 


N= 
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For the contravariant form of Q, 


T1=Tm 
= OO ies ao os Oe ae 
Og* ‘ 


MSS 


DeQire™ 


The summation in both equations with respect to s from 1 ton 
is understood and there are m terms with Christoffel symbols. 


These equations for the absolute derivative were introduced by 
M. Christoffel in 1869. 


4. Relation between TI and the second derivative of the variation 
of the point M in the space. 


dM = @,dq 


So dq’ are the contravariant components of dM. Similarly 


dé; = (dT) @ 


so (d¢'Ty,;) are the contravariant components of d@;. From 
those equations 


aM PM 
dqgi  Og’dq* a” 


and so the Christoffel symbol is directly related to the second 
derivative of the point M as it moves in the space. 


5. The absolute derivative implies the variation of the unit vector 
é; and it is indicated by D. The small d in dV’ is the partial 
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derivative of the component V* of V multiplied by the variation 
of the coordinate 


NOTE: Advanced Conclusions on Parallel Transport as it 
will be Discussed in the Next Section 


We may see that the Ricci relationship, D,g;; = 0, implies that the 
length of a vector that moves parallel to itself does not change. In 
fact, the length / of a vector V is given by 


— VV" 955 
sO 
dl = (DV*) Vigi; + Vi (DV!) giz + V'V3 gi 


We indicate the tensor or absolute derivative by D. In the case of 
1, which is a scalar, D is the same as d, which indicates the simple 
derivative of a scalar function. 

Now according to what will be discussed in the next section, the 
absolute derivative of V as it is displaced parallel to itself is null. 
Thus, the two first terms are zero and 


Thus, if Dg;; = 0 then di = 0 as well. 

The Ricci relationship applies to Riemann space and of course 
to Euclidean space, but for the Weyl space there is a more general 
relationship. In fact, Weyl space allows for a variation of length as V 
moves parallel to itself. The “connecting function”, or equivalently, 
the new Christoffel symbols, are different from the one discussed 
previously. 
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In fact, we will see that in any displacement 


Ogi; 
Oq™ 


Dn Giz = + DP nig +P img = GigPm - 


The last equation express the linearization of dl, which according to 
Weyl, is assumed when the vector moves parallel to itself 


. = ,,dq” (0.3.5) 


The term I is equal to the Christofel symbol (I’) plus a term that 
includes ®,,. The relationship, as we will see, is 


= 1 
Th =T§ — 5 (Gig + Dia - Dea) - 


We will prove this in the next section. 

We see that D,,g;; = 0 means that ®,, = 0 and so dl = 0. The 
vector moves without changing its length. This conclusion assumes 
the definition of ®,,,, see (0.3.5). Thus, in the Weyl space, the vari- 
ation of length of a vector which moves parallel to itself may be 
non-zero, and the following equation is valid: 


This may be seen as a generalization of the Ricci Theorem. 


Summary: Relationship Between the Christoffel Symbols 
and the Metric Tensor 
1 
Daji = 5 (Ges + Gik — Gui.g) 
Dagis = Dijn 


ee = oT ni 
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cofactor(gjn) 
i) 


D = Determinant of g;n 


jh 


ies (ge ae 


eae Gea 2Gtm 
gh hj ’ Im,n Oqn 


Examples of Christoffel Values for Curvilinear Coordinates 


1. Polar Coordinates (p, 0) 


Pe 20) 

9po = 0 0 

1 O 
(0 3) 

Py 


so the square of the path of an element in the space is 


As it was before. 


and 


ds? = dp’ + p*de? 
We will calculate [;;;, (¢,7,k being p or @). That is, we will find 


the Christoffel values for D po, P60, V0, le, P09. We know 
that I is symmetric in the extreme indices: [p49 = T'99,. Thus, 


iI 1 
Dopp = 5 (Yoo.r + Gop,e — Gpp.p) = 9 Jeep =0 


1 
Ppp = 5 (90.0 + Gop. — Ippo) = 9 


1 1 0p? 
I .06 = 5 (90.0 T 966,p — 9 00,0) = 2 Op =P?p 
1 1 0p” 
V66p = 5 ($00.0 T 9p6,0 — Jop,0) = 2 ap =p 
1 Lop. 


Taoe = 5 (ona + ) = 5.900 = 
600 = 5) 9Jo0,0 + Jo0,e — 9Joo,0) = 9 966, RO 
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OR 
20p 


1 
D609 = 5 (Yoo. + 90,0 — 960,p) = 


1 
ers = 5 (on.0 te Ipp,0 — I00,p) = 0 


2. Spherical Coordinates 


1 0 0 
Jos =| 0 p? 0 where i,7 = p,6,W 
0 0 p?sin?@ 


ds* = dp” + p*dé? 4+ sin? Ady? 


1 1 
Dopp = 5 (Yoo.r + Gop,e — Gpp.p) = 9 Jeep =0 


1 1 0p? 
Do6p = 5 ($00.0 + 96,6 — 90,0) = 2 Op =p 
DP ywe = po sin 7 cos 6 
Pupp = 0 
ear =/p sin? 6 

p00 = p 
Christoffel Symbols in Polar Coordinates and Their Rela- 
tionship to the Centripetal Force 


From the values of T';;,, where 7,7,k are p or 6, we determine Thy 
such that 


Pea =Piyga or Pay = oP ag - 


1 - 
We recall that g°? = 1, g°? = -, and that otherwise g’’ = 0. Let us 
p 


determine the following: T'%,, '%,, Ve 
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Mo = FT ppp + oT 6p = 0 


ey = 9° bp a7 IT 66 = 0 


To = 9 Tee + 9°’ Tape = —p 


Calling y the acceleration of a particle around a curve whose 
parameter is t, its contravariant components are given by: 


where ¢' = “ and 7 is 0 or p. So ¢? = 6 and ¢? = p. Assuming 


i = p, then 
VP = PHT ECE = P+ T596%0 
Where I, = 0 if and only if [4 =0. Thus, 


2 

cs : . U 

A pee 
p 


Where v = p$. The last term corresponds to the centripetal 
acceleration. 
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0.4 Lecture 4 


0.4.1 Riemann Space 


Let us look into the absolute derivative in relation to parallel dis- 
placement in the Riemann Space. 

Let us take one point M in the space that will be a function of 
q',q’,q°,... For simplicity, we use only q to represent q1,qQ2,... A 
vector applied at that point may keep the same component as it 
moves parallel to itself. That is, 


V'(q+6q) = V'(q) . 


Although this may be true for a given system of coordinates, it is 
no necessarily so for all coordinate systems. Let us assume that it is 
correct for the system g. This system is called a system of geodesic 
coordinates. Thus, 

V'(G + 6g) = V"(q) . 


Now, in a coordinate transformation from the system g to q, 


V'(q + 5q) = Ay(G + GV" + 59) 
where Aj, is the matrix transformation equal to 


Oq' 
Ogk © 


Neglecting second terms we obtain: 


acho OA os is 
Vi(q+q) = (A + ah") v"(@) - 
In the last term we used the previous equality V*(¢+6q) = V*(q). 
Now, since 
Vi(q) = V"(MA®) 
then 


aA\ 


Vi(q + 6g) — V'(q) = ag V*(q)oq" . 
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Thus, the vector V that moved parallel to itself has equal compo- 
nents in the geodesic, but different ones in another curve given by 
the previous equation. So the vector V moves parallel to itself as 
long as the difference of its components is equal to (replacing A‘, by 
Og’ 

agk) 


A2qi 
= V*(9)6¢" . 
Dakagh (q)dq 
Now let us express the difference in the g system of coordinates 
using the inverse transformation 


il og 
A'n(d) = dq 
So, 
VAG) = Al(QV"(4) 
and 
bq" = A’ (q)da 
Then, 
Vi(qg+6q) —Vi(@) = OF ark ih) yma) 
q+ oq q) Dqkaget mt q . 
By definition, aq! 
i tk ath 
P= —aqkage es 


I’, is symmetric with respect to m and Il. if the space admits 
curves with torsion, this not true anymore and we are in the space 
of Cartan. 

In summary, a parallel displacement of the contravariant compo- 
nents of the vector are not necessarily preserved, but they change 
as 

Vig t+ 6q) =V'(q) —TinV™5¢' « 
That is, the variation of the components between two parallel vectors 
is not zero anymore, but rather is equal to 
ov’ 
Od! 


bg = —TyV Og « 
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In other words, the total differential DV’ is 


: Vi . 
DV'(g) = or OF +TuV Og = 0. 


This result agrees with the previous section, where the Christoffel 
symbols were introduced to account for the variation of the unit 
vector é;. It tells us that the absolute derivative DV of the vector 
V is zero. 

Let us introduce unit vectors again: 


> 


V=6eV' 
The total or absolute derivative, DV (its component, ak? is called 
qd 
Dg), is: 
D,V = age + V D;,.€é; 


Since, in parallel transport, DV = 0, and we define I? as equal 
to the variation of the unit vector between two near points, D,é; = 
Prem, we again find that 

OV? 2 


Now, changing 7 for m in the first term, we get, 


ov™ 
— +VTi=0 
Og" tk 
which is the form for the derivative of the contravariant compo- 
nents of the vector V equated to zero, that we previously found by 
transporting V parallel to itself. 
We summarize out previous pages by saying: 


1. Parallel displacement or transport of a vector means that its 
absolute derivative (tensor derivative) is equal to zero. 


2. Parallel transport does not mean that the component of the 
vector is invariant except when it moves on a geodesic. 
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3. The variation of €; from a point p (of coordinates q', q’,...,q") 
to a close point of coordinates q' + 6q',...,q" + 6q”, is given 
by the expression dé; = 'e7, as q* varies. We found that as 
V =V'é& moves parallel to itself, it obeys an equation which 
is consistent with this expression for de;. Thus, the Christoffel 
symbols are related to the second derivative of the variable qs 
with respect to the geodesic coordinates. 


Example of a Geodesic and Parallel Transport 


Assume a sphere and let us move a vector V on the surface if such 
sphere: 
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Let us start at 1 with V tangent to the equator, which is a 
geodesic. Let us move V to the North Pole perpendicularly to the 
meridian a. Thus, V is always perpendicular to a and will now form 
an angle with the tangent to any parallel on the b meridian. There- 
fore, when we reach the equation at 6, V is no longer tangent to 
the equator, but rather forms an angle 6 with the initial position 1. 
So, in parallel transport from 1 to 6, V at the 6 position does not 
coincide with itself at the the 1 position. 

In a plane or flat space, things are different. If we transport 
V ona circle with a constant angle 6 with respect to the tangent, 
we always come back to the same initial vector. The difference is 
based on the fact that this space is flat while in the precious case 
it was curved. In a flat space, we may use curvilinear coordinates 
that we may transform into rectilinear coordinates where V may 
run parallel to itself with changing its components. They are the 
geodesic coordinates mentioned above. 


Curvature of the Space and Riemann Space 


As we discussed previously in the first section of 0.3.1, the variation 
of the unit vector dé; may be written as: 


dé; = (Ti,dq*)é; (0.4.1) 
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where ra are the Christoffel symbols. This was discussed in the 
framework of a Euclidean space. The Riemann Space is more gen- 
eral. In this space, there is also a metric tensor g;; and so a differ- 
ential quadratic ds? which is just a Pythagorean form between two 
close points: 

ds” = gi,dq'dq 
where (q!,q?,...,q”) are the contravariant coordinates. 

The coefficients of g;; are now any arbitrary functions of g', con- 
strained to the fact that they are continuous differentiable functions 
up to a high order. There may not exist, however, any curvilinear 
coordinates that will generate the functions g;;. 

We may still work as we have done so far, assuming a Euclidean 
space at least locally. That is, we assume that it is always possible 
to have a Euclidean plane at the point M(q',...,q") of the Riemann 
Space, which is "tangent” to it at IZ and where the Euclidean metric 
gi; at that point has the same value as that of g;;(J/) of the Riemann 
Space. 

Every tensor at M has a corresponding one in the Euclidean 
*tangent” Space. The tensor algebra, which implies derivatives, is 
developed in this Euclidean space which is taken as representative 
of a local surroundings of the Riemann Space. At the point WM, the 
algebra which characterizes the local geometry is the same in the 
Riemann Space as well as in the Euclidean ”tnagent” plane. We say 
that the metrics are osculatory. The osculatrix plane will contain 
the normal vector to the tangent plane. A curve C’ passing through 
the point M in the Riemann Space will have been reproduced locally 
in the osculatrix plane with the same curvature that C’ has at M in 
the Riemann Space. 

When moving to another point M’, another Euclidean ”tangent” 
space is considered with a different gj;, such that 


9ij3(M) = Gi; - 


The algebraic connection between a point M and another M’ = 
M-+dM, which carried out the absolute derivative, is done as before 
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by the Christoffel symbols of the Euclidean Space tangent at M. 

The curvature of the Riemann Space at each point may be calcu- 
lated by considering an elemental parallelogram based on M,M + 
dM, and M+6M: 


M+o6M <d(M+6M) x 


6(M +dM) 


M dM M+dM 


We displace the unit vector at M, &;, from M to M’ via two paths. 
Path 1 goes from M to M+6M to M' and Path 2 goes from M to 
M+dM and then to M’. If the space is curved, the variations of 
e; will be different on each path. If the difference is zero, then the 
space is flat. 

The differences are taken in the Euclidean ”tangent plane.” Call- 
ing é;' the "unit vector” at M’ and &; at M, then: 


(é' — &)o = dé + 6€} + ddé; 


(Gi — G&)1 = 66; + d& + ddg; 
where the subscripts 1 and 2 stand for path 1 and path 2, respec- 
tively. Thus, 

(Jo — (G%)1 = dde; — ddé; 
When these two elementary vector variations d and 6 commute, we 


have a perfect “plane”; otherwise these differences give rise to a 
tensor as follows. Since 


dé =T,dq'e and 66 =T4,6q'6 


52 CONTENTS 
then 
6dé; = Ti, 6d(q*)e& + 1,dq*5e + (0,T%,)dq" dq" . 
Similarly for déé; (we change r for s) 
déé, = T,d6(q")e; + 1i,5q*de; + (0,1, )6g*dq?é 
In the difference of 
ddé; — dd(é;) = (dé — dd)é; 
the first terms cancel out and we have 
(dd—dd)éj = (14, 6q*dée;—Ty,dq*de) + (0,(P4,)dq*dq?—0,. (TY, oq’ dq") ¢j 


In the last term of the last parentheses, we change k for s and r for 
k, so it becomes: 
(O03, — OT},)dq* dg? é; , 


while the first parentheses becomes, after introducing 
dé; = Tdq' em and de; = Tidq'em ; 
Cry a 40% )dq'dq*ein 


In the second term inside the parentheses we changed k for / and vice 
versa. If in the previous expression, which contains the derivative of 
the Christoffel symbols, we change 7 for m and s for /, we find that 


(45 — dd)é = [ATE — ALT) + CEG — TaD) lemdg'da® = ier 
where 

OP = [CAPE — wT) + TEL — PAT Ge) Ag’ dq" 
The curvature tensor is the expression in brackets: 


Riva -_ (Orin -_ On) + (MP ~~ rr) 
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For R7,, = 0, the operators d and 6 commute, such that (dd—dd)e; = 
0. 

By making equal the indices m and k we reduce the four -index 
tensor to a two-index one. The result is the tensor Rj, called the 
Ricci Tensor. 

It is not difficult to show directly, using the coordinate trans- 
formation, that 92 is a mixed tensor and that Ri is as well. By 
simply recalling the definition of a tensor with regards to the coor- 
dinate transformation, you will find that both elements fulfill those 
conditions. 

The tensor R77, is called the Riemann-Christoffel Tensor, and 
may be written with four covariant indexes as Rpj;4, by multiplying 
Rei, BY Gmn- Let us rewrite the Riemann-Christoffel Tensor in the 
following way: 


Rin = (OVE + TEP G) — (Ue k) 


where the last parentheses signify the first expression in parentheses 
with the indexes / and k interchanged. Thus, 


Rnitk = Gmn((Ol'R +140) — (Le ky] 


We recall that 
Ogenki al tah 


Therefore, 
OL Imn) = (Pink) = ALE) Imn + VOImn 
So, 
Raiae = [O(Vine) — 0% (Vinm + Damn) +00 jn] — (Lo ¥] 
which simplifies to 


Rrite = [OV ing — TET mma] — [1 & &] 
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where 


Dink = UG¢Ghm, ete. 


A T with three sub-indexes (TI...) is called a Christoffel symbol of 
the first kind, while a [ with one super-index and two sub-indexes 
(I™.) is called a Christoffel symbol of the second kind. 


Space With Normal Coordinates 


As was mentioned when the concept of the geodesic was discussed, it 
is always possible to generate in the Euclidean Space a local rectilin- 
ear system of coordinates tangent at M called normal coordinates. 
In that system all metric tensor components g;; are constants, and 
therefore the Christoffel symbols at that point are zero. In this 
system the curvature tensor Rp; 1, is reduced to 


1 
Raith = 5 (AGignk + AWOngia — NOnGik — OO:Gni) - (0.4.2) 


In this equation, whose origin is the variation from point to point 
of the Christoffel symbols, we have used the expression for those 
Christoffel symbols as a function of the metric tensor. 

The symmetric relationships can be easily seen: 


Rritk = Riz,ni 


Rhitk = —Rrie = —Rinar 


and by circular permutation of the indexes 2,1,k, it results that 


Raith + Rniki + Raza = 0 - (0.4.3) 


Ricci Tensor 


At this point, it will be illuminating to review tensor contraction 
that reduces a tensor to that of two indexes. For any mixed tensor, 
equating a contravariant index to a covariant one means that we are 
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contracting those indexes or reducing the tensor to two orders. For 
example, 


Tt =O) He (0.4.4) 
k=1 


The parentheses emphasize that repeated indexes imply summation 
from 1 to n, where n is the number of dimensions in the Riemann 
Space. In relativity, n = 4. 


J 


The Ricci Tensor R,:, whose source is the curvature tensor R;,,,, 


is defined as 
Ry = 0 Re ay (0.4.5) 
j 


R,. is equal to Rj, because 


Ry = GOR = Ore” Riva = GP Ria = Go Ring = Rr 
(0.4.6) 
In the third step, we exchanged ml — hk, ml > Im, and hr > rh, 
and then used the symmetric properties of Rag ys. 
The Ricci Tensor may be reduced further to a simple scalar R by 
multiplying it by g”: 
g" Rn = R 


R is called the scalar Riemann curvature. 


Bianchi Identity 


The equation with which we previously defined the curvature tensor 
Ri, could also be written as 


where D, indicates the absolute tensor or covariant derivative with 
respect to g!. Now introducing the operator D,. in (0.4.7), we have 


If t;; = eje;, the external product of the vectors é; and €;, then 


DDg(tiz) = (DiDr ei) ej + ei (DiDpej) + (Drei) (Diej) + (Drei) (Drej) 
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Changing | + k, and the last two terms cancel out by subtraction. 
Thus, 


(DD; = D,D))ti; = €;(D,D;, — D,D))é; + €;(D, Dz = D,,D))é; 
Now, using (0.4.7), 


Operating on the left-hand side with the differential operator D,, 
we get 


In (0.4.9), we assume t;; to be a tensor equal to the covariant deriva- 
tive of e;. Thus, t;; = Djé;, and (0.4.9) may be written as 


We may eliminate the second term by cyclic permutation of 7,1, k, 
by summing all three equations, and then by using the property that 


Riin + Re + Pikg = 9 (0.4.11) 


The left-hand side of (0.4.11) is identical to (0.4.8) if we change r 
by j in (0.4.8). So, the left-hand side of (0.4.11) becomes 


Dj (Rep )em + Riin(Djem) + cye. perm. of 7,1, k (0.4.12) 


= Ri,(Djem + eye. perm. of j, 1, k 


The second term of (0.4.12) plus its corresponding cyclic permuta- 
tion cancel with the term on the right side of (0.4.12). What remains 
is 

Dy (Rein) em oF eyes perm. of 7,16 


Since e;, appears in all the terms, then 


D,.R? 


alm 


+ DR mp + Dini; = 0 (0.4.13) 
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This is called the Bianchi Identity. Thus, the curvature tensor R?;; 
fulfills the identity of (0.4.11) as well as that of (0.4.13). This is a 
third order derivative of the vector €;, whose variation gives rise to 
the curvature tensor. 
By substituting m = j and contracting the index i with /, we 
obtain 
2D, R! —D,R=0, 


which becomes, after using the Ricci Theorem, 


1 
D(Ri.— 5 IR) =0. 


By introducing the symmetric tensor S;,,, 
1 
Srl = Ry = Grif ’ 
2 
we obtain a conservative equation for this tensor: 
Disa) je. 


This equation plays an important role in the Einstein interpretation 
of the gravitational fields of General Relativity. 

Einstein looked into the metric tensor in the 4-dimensional Rie- 
mann Space, gi;, which is a function of q°,q',q?,q¢°’, as the source 
of the gravitational potential. He proposed an equation similar to 
the Laplace and Poisson Equations which expresses the conserva- 
tion of gravitational energy. Thus, the tensor of the second order, 
say S),,, that Einstein uses as representing the gravitational energy 
must obey the Laplace type of equation 


Say =0 


if no other actions exists on the particle but the gravitational one. 
This equation implies that the second order derivatives of the metric 
tensor are equal to zero (see (0.4.2)). Thus, this is a type of Laplace 
equation for the g,,, viewed as potential functions. 
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With other energies besides the gravitational, like electromag- 
netic or microstresses of the given fluid in question, the equation 
must be a Poisson type, or 


Diy = Kies ) 


where the tensor ()),, represents all other actions and y is a universal 
constant. Einstein uses the following expression for 5},,: 
1 
ine = fe Ink 
which we show holds the Laplace Equation. For Q),,, Einstein visu- 
alized the medium as a continuous one with charges particles. Thus 
the energy is represented by the sum of the ”mechanical” impulse - 
energy tensor P™ and the electromagnetic impulse - energy tensor 


M yy 


1 
P™ = purul + a (0.4.14) 
and 
1 sl ap . 
My = 5 (Jou pgs FraFe ; (0.4.15) 
Thus 
Qau _ Mu = Pha : (0.4.16) 


The tensor T*” represents the mechanical flow of energy transmitted 
to the surface of an element of volumen of the continuous medium 
which surrounds the point P that moves with velocity u* in the 
4-dimensional space. In the expression of M),,, F “8 ig the electro- 
magnetic field. 

Thus, the field equation is written as 


1 
Ryu _ 5 Imukt = XQDaru . 


Since DS), = 0, then D,,Qa, = 0, expressing the conservation 
law for the other fields. 
Particular values of the tensor Q),,: 


0.4. 


LECTURE 4 59 


. No electromagnetic field and no stress in the medium: the only 


energy in consideration is that of the mass motion 
Q™ = pu*utt (0.4.17) 
d dq” . . i 
where u* = —— and p is the density of the matter with respect 


Ss 
to the frame at rest. 


Case where the medium is a perfect one; that is, there is no 
shear, only pressure P: 


if Xr ie Xr E Xr 
a Pee od + Qt ur 
or P P 
OS (0 - 3) — a (0.4.18) 


. Only electromagnetic field, no matter: 


Qu. = My, (0.4.19) 
In the first case, D,,Q*’ = 0 means 


D,(pu")u* + puD,u* = 0 


d r 
Since ud = it is (uAu“ = 1) and tangent to the curse s, 


S 
whose differential is ds. This curve is by definition the geodesic 
of the point whose mass is p. By being u* unitary, 


u"D,u* = 0 
and therefore, by calling J = pu the "mass current,” 
Dy Ie=0 
Multiplying by wu, > 
i Di hy= Dis =O (0.4.20) 


which tells us that the divergency of the current is zero. 
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0.5 Lecture 5 


0.5.1 Weyl Space 


Variation of the Length of a Vector as it Moves Parallel to 
Itself: Weyl Geometry 


In the previous sections we found that a parallel transport of a vector 
in Riemann Space was associated with the nullity of the absolute or 
tensor derivative. For the contravariant component this means that 


Our dees 
Ogk +U Vk = 0 
The Christoffel symbol T°, may be associated with the variation of 
the unit vector €; as it moves from P to P+ 6P, two close points in 
the space. Thus, 


dé, = Ti ,é;dq* 


As we remarked before, I, is not a tensor but rather a function that 
links two close points, like a bridge in the space. As expressed in the 
previous section, [ represents the linearization of de; in the space. 

Weyl introduced another function that represents the lineariza- 
tion of the variation of | (the length of a vector); explicitly: 


dl = ©,,dq" 


The existence of the vector function ®,, (a.k.a., the differential po- 
tential) implies the possibility that 1 = vv, varies from point to 
point. The Christoffel symbols will be different from those in the 
Riemann Space. They have to include the differential potential ©,,. 


Transport of a scalar / around a closed curve 


Let us assume that the scalar ! moves from P to Q via the paths a 
and b. We will look for the variation of | as the paths vary. 
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The differential change of | may be written as the scalar product 
of two vectors: 


Mato) KD) _ t= ode 
The circulation around the curve will induce a work of ®,. We 
ssume that this vector ©, moves parallel to itself. 


The circulation from P to Q via a and back to P via b is given 


by 
b,dq" +f &, dq! 
foo Bulg en 


By the Stokes Theorem, the circulation may be expressed as the flux 
of the rotational of ©, on the surface inclosed by the curve 


f Pag" = [fo x ®,) ds 


If the rotational is zero, the function /(q) is simple and there is 
integrability. 


If the rotational of ®,, is different from zero, then 


(la(Q) — la(P)) + (lo(P) — b(Q)) # 0 


In this case, / at Q is different from / at P, and thus the difference 
of the circulation of ®, (the differential potential) will depend on 
the path. We say that ©, is not integrable. That is, 

0%,  O®, 

One Ox” 
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The scalar | may be written as the scalar product of two vectors, 
u and v, in the space whose metric is gj: 


= gla 
Lew ga 


and so 
[+ dl = (u* + du*)(v” + dv”) (Gur + Agu) 


As we have shown previously, the Christoffel symbols are related to 
Juv and its derivative, so ®, also has to be related to them. Let us 
see how: neglecting second-order terms and indicating with ”d” the 
derivative of a scalar (not a tensor derivative), 


a (dar gig re” gag OO dG (0.5.1) 


dl 
Now, using races ®,dq’, 


dl = 1@,dq? = uv" g,Pp,dq? (0:52) 


Before equating these equations for dl, let us introduce the general- 
ized Christoffel symbols I’. by the following: 


du, = dq? uy, 


and 
dy= Ty dq?e" 


These equations imply that the absolute tensor derivative D,u = 
u,v" = 0; in other words, the parallel transport of wu and v. 


Since dl is a scalar, dl = it and 


di = org dur ur Gade? or ud gay 
From the previous definition of I, 


du" = Ty ,v°dg’ and d,ul = —Tygu" dq? 
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sO 


ct =U" Gl guedg? — Gwt"T vd’ awards (Gita) 
On the other hand, by linearizing dl as “ = ®(q)dq, 
dl = Idg*®, = vu" Gu Phage 
Comparing these two equalities for dl, 
—v"u? g.T gq? — Guw'v°T, dg? ur ord a Gi) = v’uH Gur Phage 


In the first two terms, we use the symbols I... as 


V — 


Gil a= we and cro Dee ere 


Substituting into the previous equation, changing k by p in the 
second member, and then dividing by dq?: 


ha ae O 
<0 wT og — uv"T ap + eer = VU Gur ®p 


Now by changing 8 by yw and a by v in the first two terms (they are 
repeated indices) and dropping v’u", we get: 


= a Og . 
Vv — Pupp + ae = Iwp (0.5.3) 


If ®, = 0, this equation is the same as that which relates the orig- 
inal Christoffel symbols to the variation of the metric tensor. In the 
case dl! = 0, that is, the scalar vector of u and v as they move parallel 
to themselves does not change. If &, 4 0, the generalized Christoffel 
symbols I will depend on the variation of the metric tensor and the 
differential potential @,. Soon, we will see that the metric tensor 
may change with the gauge, but I’ will remain constant. 
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Generalization of the Ricci Theorem 


Let us differentiate the length | of a vector €, such that | = €'€/g;;: 
dint = (Din€ VE? + €9:;(Dm€é") a5 C&D (Giz) 


By definition of I, and assuming that the vector & moves parallel to 
itself, 


and similarly for Dé. So, dnl = €'€/Dm(gij) is zero if Din(gi;) is 
zero. 
Now, taking g;; as a tensor of the second order, 


O Gig =l =k 


Din( Gj) = ag 1 mi Gil — Ding Gik 


and according to (0.5.3), 
Dr (Gig) = Gig Pm 


We may say that in the Weyl picture, this is a generalization of the 
Ricci Theorem (for parallel transport only, but Ricci is general in 
the Riemann Space). 


Invariance of the Angle Between Two Vectors as They Move 
Parallel to Themselves 


From the previous equation, we have found a useful assertion: the 
absolute or total derivative D of the covariant component of € is not 
zero as Dé’ = 0, but rather 


Dingi = Din (E" 9:3) = (Din€) (9:3) ate @ (Dini; 


The first term is zero, and for the second we use the generalized 
Ricci Equation: 


Dae: =© (Digs) =O Ons 
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and thus the absolute derivative of the covariant component is 
Dag On Gij =0 


Let us apply this theorem to find the angle between two vectors 
in parallel displacement: 


where J; and /, are the lengths of € and 7. The total derivative of 
cos @, which is a scalar, is 


Vi ly Viyls lily 


The second terms is zero. The first is 


(Dn&i)t — 'E Omgij 


Vile Vlils 


D,, cos @ = 


The third term is 


i Lie (1) dh) + (dla) 
gin Pa =1 &( 5) (ita)s 


® dl 
Since d,,pJy = poe = 1,9,, and dmlg = lo®, where dl = —., 
etc., the third term becomes 


Ph oa 


(habVhb vVhb  vVJhh 


which cancels with the first term. Thus, D,, cos? = 0, and @ is 
invariant. 


Relationship Between Weyl I and the Christoffel Symbols 


0 ss = 
Oger i = Gig Pe + Vagi + Vix 
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and similarly, 
= x O = = 
agi ok = gi2j5 + ine +T ja, and agi sk = gies + DU jus + Page 


Summing these three equations and dividing by 2, 


(ngs a O59 = O:gni) = : (gi; Ph, + Ii j = Gjk®i) + Vyin + Uji a P jit 
2 2 2 
The first term is, in the Riemann Space, the value of Ij, of the old 
Christoffel symbols. Thus, the new Christoffel symbols and the old 
ones are related by 


1 
Dik = Djie — 5 (ii Pe + griPj — 9jn®i) 


Curvature and Generalized Christoffel Symbols 


Let us define a vector v that moves parallel to itself in a closed curve, 
that may lie on a curved surface, and such that it is related to the 
connecting function by 

du' = ®idq" 
We add 7 to the set of ©, to identify in it the given contravariant 
component of v. Since it moves parallel to itself, 


dy’ = Twdg* 


The circulation of v on the closed curve will, in general, be different 
from zero. So the final variation of v is 


: : 1 . ; 
dv! = f Ojdq' = 5 |] (i. — Pi,)ds" 


This variation will be, as done before in the Riemann Space, pro- 
portional to the curvature of the surface: 


: 1 ; : 
dv' = 5 if Fi gwids 


where 5 5 
i Aue pata 


m 4 m 
pie al — Vril gh 
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0.6 Lecture 6 


0.6.1 Gauge Theory 


Introducing Different Gauges, or ”Etalon,” in the Same 
Space 


As we will see in this section Weyl assumed the possibility of chang- 
ing the length of the vector by changing the unit of measure, or 
”etalon” in French. This will correspond to two different Riemann 
spaces superimposed in one. Thus, when a vector moves from one 
point P to a closed one Q while remaining parallel to itself, the 
metric tensor g;;(q) changes because the coordinates g change. This 
change corresponds, as we have seen in the previous section, to the 
curvature of the space. But if we change the ”etalon” this brings a 
new reason for the variation of g;;(q). If g;; changes by a factor A as 
the etalon” changes, then gi; > gj; = Agi. If I(q) is the length of a 
vector € as it moves parallel to itself I(q) changes to I/(q) = Al(q). So 
the connecting function that we have called 6 = — in the previous 


section and which is directly related to the generalized Christoffel 
symbols I also changes because 


di! = \dl + 1dd 
os ait. ar ax dd 
ge aaa eas 


Thus the question regarding this general parallel displacement dis- 
cussed by Weyl is how the whole approach changes as g,,, is multi- 
plied by a positive factor A, which may depend on the coordinates. 

The element of space at the point ds? changes to ds? = \g,,dq’ dq" 
The factor A is an indication of a change in the length of the unit of 
measurement, “etalon”. It is as if the ruler was changing its “length” 
as we go from one point to the next, like it was made of metal and 
at each point of space-time the temperature was different, causing 
dilation and contraction to occur. (The fact that the physics does 
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not change as we change g — Ag is obvious for the light transport 
from a point P to P + dP where ds? is always equal to zero.) 

Let us discuss gauge change from another point of view. Let us 
call EF the “etalon” at P on the geodesic path: with its help we 
measure the length of €;. Now call E’ the etalon” at P on the 
nongeodesic path: with its help we measure the length of &;’. E is 


called the ”geodesic etalon” and is different for different points in 
/ 


E | ae 
the space. Thus, |é’] = leila and gi, = & - 6} = (G- a(z) 
/ 
If we say that \’ = (Z). then g;; — gij’ and the length of & 
changes from 1 = €'€g;; > NE"E%g;; = Xl’. When we move from 
P to Q, ra will conserve its components, and of course its length, as 
long as the “etalon” is invariant. That is, 


(Q)=UP) and €(Q) =€(P) 


When A is assumed to be a function of the point in the space-time, 
from P tp Q, then A will vary by d\. On the other hand, by changing 
the "etalon” from FE to E’ at P itself, the length will vary at P by 


7 


E'\2 
a factor of (=) = 2 that is 


i= DS SP 


Then, by moving to Q (a close point to P), l’ will vary by 


dl’ = Id + Ad] 
SO 
a tas 
and a = “ -- “ is the variation of /’ when P > Q per unit of I’. 


dl 
As in the previous section, let us call — = d®, where ® is as- 
sumed to be a continuous function. As P moves to Q, this function 
varies as 


do! = db + “ = d® + d(log X) 
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The vector ®,, may be seen as a field that is generated as the “etalon” 
varies in the space-time worlds. Thus, in a parallel displacement, 
if gj; changes to Agij, the length changes and at the same time 
the differential of the connecting function d® is changed by the 
differential d(log A). The differential function d® may be expressed 
in the generalized coordinates by its components ®,,, such that 


dé = @,,dq" 


If the circulation of ®, around a closed path is not null, then the 
rotational of ®,, generates a new field which is independent of its co- 
ordinate system, that is, invariant under coordinate transformation. 
®,, depends on the change of the “etalon” as well as the transport of 
ra from P to Q. Thus, the change of ”etalon” changes the quadratic 
form ds?, but at the same time, in a parallel transport the differen- 


dA 
tial connection function d®’ changes by the term —. 


d 
The equation d®’ = d® + “ = d@ + d(log A) may be rewritten 
as 
@ dg" = &,dq" + d(log A) 
or 
/ d(log A 
@, = OT (0.6.1) 


If instead of \ we use e%, then this equation becomes 


On 


fe a 
Out Ga 


be 


Introducing the variation of g, we see that 


= + 
Gij A 


do’ 


and 


O(log X) 


o! - 10g}; 
Oq" 


Mg Og" 
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Next we will show that ine are invariant under changes of “etalon” 
Before, we found that g,, T° 


wp, and the potential ®, are linked dur- 
ing parallel transport of a vector by the expression 


OGiiw = = 
Bg? = Ip + Vor + Pupp (0.6.2) 
Now, if the ”etalon” changes, then 
Duy Duis = Nw 
and a 
®, => D, = ®, + dg? 
Then (0.6.2) changes to 
OG wy dX 1d LT LT 
A Age + Iu op = Agu(®, 4 x ? T* pup ' ~ vup 
or 3 
Dw =/ =/ 
A age AIuv Pp + Dou + Dyp 
Dividing by 4, 
Ow Rivas by 
oe = Gyn By + Ee (0.6.3) 


= =k =k = =i 
Substituting Le, > bey = Bs Gey, and since Pu; = [;,;9i, the 
last term of (0.6.3) becomes 


1 xk —rk 
xr Mme’ wade) 


Similarly, the last term of (0.6.2) becomes 
= = 
(Poi Skv se DP Gku) 


Comparing (0.6.2) and (0.6.3), they are identical if ee 


=T, px BbG. 
So, I. is invariant under gauge change. 
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Invariance of the Absolute Derivative of £' 


We have seen that g;; gj; leaves ne invariant. Thus, if the coor- 
dinates do not change (only the gauge), then D,,€' remains constant 
and in parallel transport is equal to zero. 
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Invariance of Curvature as Gauge Changes 


Let us show that the curvature tensor defined from the geometric 
idea of parallel transport is invariant under gauge change. 


P, Pre 


dq? 


In the above infinitesimal parallelogram, we may start at point 
P with our transport of a vector v and reach P,, via the path a 
or b. At Piz, we compare the two vectors: in general, they will 
not coincide except if the parallelogram is on a plane. We take the 
difference between the two vectors as the element proportional to 
the curvature: 
Va(Pi2) — Ue( Pig) = Av 


Then, if Av may be written as the product of the curvature tensor 
per unit of surface AR’, such that 


Av" = ARfv" 


In general, AR# is proportional to the surface of the parallelogram, 
and so 


AR’ = RY dq’dq? 


vpo 


Since the contravariant components of the vector v, v“ and Av", 


do not change under a gauge change, then Rf and Ri, are also 


invariant. 


The tensor R!,, is antisymmetric in the interchange of the indices 


pando. This is so because the element (Aq)?? = dq?dq?’ refers to an 
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oriented surface, and so changes sign when it goes counterclockwise 
rather than clockwise. Thus, 


Aq’? a —Aq’” 
and 
iui = =i 


Thus, if we use the subscripts p and o only once and ignore the 
order, then 


The curvature tensor R? po 18, in this approach, generalized from 
the gravitational field Pf,, as well as from the electromagnetic field 
vector F’,,. It may be seen that the additive contribution of both 
fields may be expressed as 

a, 1 
ew ae 5 ov Fp 


vpo vpo 


where 6! is the Kronecker Function: 1 or 0. The tensor Pf, is 
antisymmetric in both pairs: pu,yv and p,o. 

Obviously, when the electromagnetic field is not present, Riis 
corresponds to the classical Ricci Tensor and parallel displacement 
preserves the length of the vector. However, when F’,, is introduced, 
we do not have a unique geodesic and the integral will depend on 
the path taken. In any case, the change in gauge or etalon” does 


not change the curvature tensor. 
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Appendix 


1. Some Remarks as to the Differences Between Euclidean, ag 
Riemann, and Weyl Geometries 


A space is an object that becomes expressed when a set of 
numbers is associated to it, and consequently becomes a nest 
of functions. Thus, if the set R of real numbers is associated 
with a 3D space, then we have the usual FE? space, and a point 
in the space is a set of real number (x,y,z). A function of 
f(x,y, z) will be defined as a mapping of a region A of the 
space, which corresponds to a subset of R that goes over to 
another region B of R, which is occupised by the numbers r; = 
f(r1, y1, 21), T2 = f (2, ya, 22), etc. This mapping is expressed 
as 
FLA SE” 5. 2, 9%5eA 


A very close definition to this is that of a patch, a concept which 
is important in topology. In this case, a 2D subspace, described 
by the two real numbers (2, y), is mapped over a space E?. This 
final mapping is called a surface in EL. It is expressed similarly 
as 


f: A> EB 


These definitions are independent of the space itself, and may 
be used in any other space. 


A space with geometry in it is a space where distances and 
angles may be measured. Because of this, it is traditional to use 
the dot product as the fundamental tool. Euclidean Geometry 
is built with the dot product between two vectors: 


> > 


ae 
Ul + V2 = Gijt eo 


The Euclidean Geometry preserves the elements of the space 
as curves, triangles, surfaces, etc. Those that are congruent to 
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themselves are isometries, or mappings that preserve distances 
between points. That is, if F’ is the mapping and d indicates 
distance, 

fi. SE 


d(F(p), F(q)) = d(p,@) 


for all points p and q of E%. Translation, rotation, etc. or 
a rigid body are examples of these isometries. The figures 
are called congruents because one can be superimposed on the 
other, and they are identical although located at different parts 
of the space. Among the properties preserved are the geodesic 
and the curvature; both concepts were used by Einstein in the 
Minkowsky Space of four dimensions in order to express the 
gravitational laws in General Relativity. 


2. Riemann Space and Riemann Geometry 


As in the Euclidean Geometry, the concept of elements is pre- 
served by isometries. From a point of view of differentiability, 
the Riemann Space can be visualized as an infinite number of 
Euclidean Spaces; that is, a space where each point point has 
a neighborhood that may be taken as a Euclidean Space. In 
addition, the passage from on region to the next is smooth, so 
differentiability exists in the local region as well as in the over- 
lapping of these regions. Another way to express this fact is by 
referring to Riemann space as a Compact C'® Manifold with 
a metric g;; that is greater than zero, differentiable, symmet- 
ric, and bilinear. ”C'°” expresses the fact that it is infinitely 
differentiable. 


In the four-dimensional Riemann Space, the geodesics are also 
world curves which represent the history of the particle in con- 
sideration. Those curves are tangent to the velocity vector in 
4D, and they are also the stationary path between two points. 


The algebra of the Riemann Space V,, is carries out in the 
tangent or osculatrix plane. At each point of V,,, we create this 
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plane where the tensor algebra brings the same results as if we 
had used the same narrow interval around the point of contact 
of the space. If we use a set of coordinates y', y?,y?,...,y” at 
the point M of V,,, then there is a point m (the contact point) 
in the tangent plane that follows a Euclidean algebra such that 


(a) The fundamental tensor g,,, at M in V,, is the same g;; at 
the point m in the Euclidean plane €,. 


(b) The first order variations of these metric tensors, when we 
move from M to M+AM and therefore from m to m+Am, 
are the same up to the second order in the variation of y’. 
That is, we neglect (y’(m) — y’(m + Am)). 


(c) The coordinates y“ in V, at M become y' at m in the 
Euclidean tangent space. (We use the Greek scripts for 
the Riemann Space and the Latin scripts for the Euclidean 
Space.) The y’ in €, become curvilinear coordinates. 


Thus, according to (a) and (b), we assume 


(9i3)m = (Qu) 


(m= CGE) 
dy* J m dy* / 


This assumes that the Christoffel symbols are the same at M 


and 


and m, and also that the parallel transport of a vector M to 
AM agrees with the parallel transport of a vector from m to 
Am in €,. Thus, 
. Ov? : 
(Dyv") us = (Dyv") mn = et ry 
OYk 
If the point M describes a curve C’ in V,,, there is an image 
curve y in €,, where at every point M g,,(C) is equal to gi;(7) 
at m over y. This implies a multiple of Euclidean planes at the 
different point M of C. 
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The fundamental portion of arc at C is the same at M and m. 
That is, 
ds” = gy,(M)dy"dy’ = gi;(m)dy'dy? 


The curvature of the Riemann Space at M in C is found from 
the Christoffel symbols at m going from m to Am in the Eu- 
clidean Space ¢€,. This assumes that the second derivative of 
the fundamental tensor is the same at M in V,, and m in én. 


3. Summary of Weyl Geometry Approach to Gauge Change 


We start by scaling the metric tensor g;; by a factor of A, such 
that gj; = Agij. So, I’ = Al, where / is the square of a vector 
v — dl = work of field = ©,,dq". In order to preserve, during a 
parallel displacement which yields a change in the magnitude 
of the vector, the same field values and the same geodesic and 
curvature, the variation of | is written as the work of a poten- 
tial field ®, that does not affect the original field if it always 


d 
expresses the change per unit length of TT: 
Weyl Space is a generalization of Riemann Space. Now, be- 
tween two points in the space there is not a unique geodesic. 


Taking a vector v at P a moving it parallel to itself, we reach 
the point Q via path a. 


This path is a geodesic where, if 7 is a parameter describing the 


vector v4 = i: tangent to it, the equation for the geodesic is 
3 
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such that the absolute derivative D,v" = 0, or 


= d? qt dq’ dq? 
Pare = 

Beeler Cdr 
But if we now allow the length of the vector V" to vary, we 
may choose another path b that goes from P to Q and still 


0 


have a geodesic which holds similar equations. In this case, the 
integration of v’-V,, around the loop a+ 6 starting at point P 
and returning to P is different from zero and equal to the work 
of a function ® which in principal is different for each loop and 
is part of the new Christoffel symbols. Next, as we change gy, 
to Agu,, some tensors may vary, while others may not. We have 
seen that the curvature tensors as well as the electromagnetic 
tensor do not change under this gauge change. 


In the development of Weyl, the parallel displacement of a vec- 
tor as that given by the absolute derivative of the contravariant 
components is again equal to zero, such that Du“ = 0. But now, 
the I’ symbols are not equal to the Christoffel symbols I, but 
instead a similar equation holds, namely, 


OUp a — 

Bge 1 mee" 
for the covariant, and 

du’ 

eae Ppt 

dg?" se 


for the contravariant components of vector v. This is a simple 
extension of Riemann Geometry, but I, are not only functions 
of the metric tensor g;;, as before, but also of the potential 
function ®. 


In order to preserve the ”connecting” function I’), in the above 
formula as we change the gauge or etalon, ®,, must vary with 
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the change of etalon as follows: if g., — Agu, then 


O log » 


sa any a 


This will preserve the objects in the 4D geometry, like the cur- 
vature and the circulation of a vector parallel to itself, which 
eventually will be used to express the physical laws. 


According to Mie, all the laws of nature correspond to an in- 
tegral of the ”action density,” which must be specified for each 
phenomenon. This integral must be stationary as the poten- 
tials (or state functions) of the system vary. It is assumed that 
the 4D domain where the law applied is surrounded by a sur- 
face where the potentials vanish. The action is a simple scalar 
expressed as 


Action Density = S = / W dQ 


where 2) is the element of volumen in the space and W is the 
action density. Since dQ x ds? & gisgyda'da*da! and gu, > 
Agu, then dQ must change proportionally to A? and Q must 
change proportionally to A~? in order for S to be invariant. 


We are going to pursue this line of work in future chapters. 
But before we leave this discussion, we would like to outline a 
parallel between the Weyl approach (or the gauge variation es- 
tablished earlier) and the Theory of Electromagnetism. In clas- 
sical Maxwell Theory, the relevant objects are the field strength 
E, B, and the charge e. However the four potentials vectors 
(po, &;, By, 3) are introduced to make the theory mathemat- 
ically sound. 
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0.7 Lecture 7 


0.7.1 Analytic Dynamics 
Variational Calculus and the Euler Differential Equation 


The extreme of a function belongs to Ordinary Calculus. The ex- 
treme of an integral belongs to Variational Calculus. 

Virtual changes of a variable of a function, which represents a 
real situation, are different from the actual physical variation. 

For example, a ball at the bottom of a well with no external 
forces but gravity acting on it does not move, but it may suffer a 
virtual displacement du from B to A and we see it as a mathematical 
variation which implies a variation on the function F'(w) such that 


ou 


B 


The following is an example of a stationary problem involving 
a finite integral: the J. Bernoulli (1667-1748) curve of minimum 
descent time. 

We look for the curve y = f(x) such that the time t of descent is 
minimum. From simple mechanics, 


i= foe _ 


The heights a = f(a) and 6 = f(b) are given. This is an example 
of a variational problem. We look for y(x) that makes t a minimum. 
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Euler (1707-1783) found a differential equation to solve the ex- 
treme of a integral between two point. Let’s look for y(x) such that 


P. 
[, fe-ula).y'(@)) ae 
presents a minimum (or in general, an extreme) between two fixed 
point P, and P». 

Solution: Assume yo(x) is the solution. Take any close-by curve 
y(x) between P; and P. It is assumed that y(x) = yo(x) + € d(x), 
where 6(x) is an arbitrary function that is zero at P, and P», and € 
is a small parameter. The functions 


He) = f° Fle, wole) +6 6(e), vite) +e 8x) de 


should represent a minimum (or in general, an extreme). Let us 
make J’(e) = 0. Since P, and P are fixed, then 


ae | 7 [fvd(x) + fy5'(x)| da =0. 


P; 


Integrating by part the 2" term and knowing that d(x) = 0 at 
P, and Py, we obtain: 


d 
[fu — (fu) ]6(@) = 0. 
Since dz is any function, then 
d 
ty = qt) =0, 


but 


d 
aa (@ Y; y’) = Fut a figy a feat : 


Thus, the Euler equation may be written as: 


fy — Fue + fy’ + fyyy") = 0. (0.7.1) 
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Let us generalize this to the case where the function may depend 


on several functions y(x), say yi, yo, ..., etc, and the corresponding 
first derivatives y}, yb, ..., ete. 
The function J(e) will appear: y; > y?+e d(x), yi > y,ot+e 6'(2), 
etc, with i =1,2,... and y? being the assumed solution. 
P2 / / 
J(e) = of f(z, yite d(x), ygte d(x), ..., yP+e d(x), yo te O(a), ... 
Now, 
aI(e)_ 
( de = ~ 0 Gy 
Pot of Of Of Of 
d(x) 4 a(x) + ( d(x) 4 s(2)), ...] dx=0, 
Jos ((Syp0)* 5yg® OD) + (Spel) + Seg 


Each parentheses is like the previous case, so this will mean that 
the condition implies that 


Ps (fyj..0 oe ta Be fut V1) =0, 


ia (fos 0 fytaoVo ao fut, Yo) =0, etc., 


or in general: 


Fu = (Sut oT fut aydi a fyi) =0. (0.7.2) 


This formula may be applied to the motion of a particle in an 
n-dimensional space. In particular, it will be applied to a four- 
dimensional space to deal with the general relativity of Einstein. 
We will obtain the geodesic and also bring in an interesting rela- 
tionship of the Christoffel symbol versus the derivatives of the fun- 
damental metric tensor g;;(a) as a function of the four-dimensional 
coordinates. 
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Application of the Euler Differential Equation: Geodesic 
from the Variational Point of View 


A geodesic line is by definition a line whose integral between two 


point P, and Py is stationary: 


P2 
) ds =0 


Pi 


Assume a particle moving along the geodesic and consider the 
coordinates of the particle x’ as a function of one parameter 4, 
x =x'(\). And instead of applying variational calculus to 


d*s = gi;(x*)da'de! , 
we will be doing it to 


Bec si tay a) (a2) | 


In order to apply the Euler equation, we look at d?w as a function 
of \ where the functions y;(\) become x‘(A) and its derivative y/(A) 


is a etc. The parameter \ plays the role of x in the previous 
derivation. The function f(A) is identified with d’w. Thus, 

dax* (A 
=F OC; a )) From this, 


eS) 
dak FEV) Van]? 


a (oF) = “(2 5.) 
Aarts ak ea 


The factor of 2 is due to the fact that i and 7 are exchangeable, 


and 


since we assume a symmetric tensor g;;. 
Now, 


d dx’ dx? dz? dri 
5 (29 rh) = 2955.0 Gy Gy ost t G4 yas 
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The Euler equation reads: 


jk + Gij a2 


ode (4, da? det , toy) 
Faia, ah. NO a jk 


By changing o for 7, and exchanging j for k in the first term of 
the right-hand side, we get: 


da" = datdz’ 1 - da’ dx =f 
Ht go Day. aye Gh ay gx. 


or 
da! 1 da’ dx 
Gu que + (Gis ~ g946) ay ay = 


Since j and 7 are interchangeable, gj.,; = gjx,;, and 
1 
Gik,j = 5 (Yik.s TeGghays 


Thus, the Euler equation becomes: 


d’zi 1 dx? dxi _ 


Gij ae 5 (Yik.s b Ojki — 93k) ay ay = 0 hOst3) 


The term in parentheses is the Chirstoffel symbol Tix; .. 


1 
Ding = 5 loins Ong Oia 
Multiplying by g'*, we have a new form of Christoffel symbol, 
Vins = ie . 


By direct analysis and since gj; = 9ji,x, it is easy to see that Ti,; 
is symmetric in 2 and 7. That is, 


Dikg = Djni - 


We must observe that Christoffel symbols are not tensors, be- 
cause under a change of variables they do not change, as a tensor 
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does (called in our previous lecture an element). They are not at- 
tached at a point in the space, but rather they indicate a bridge 
between two close points. However, we may lower or raise their in- 
dices, as we do with tensors, by using the metric tensor g,;. Also, if 
a tensor varies in a coordinate system it does in any other; this is 
not so for Christoffel symbols. 

Thus, the equation of motion of a point on a geodesic may be 
written as: 

dx | dada’ _ 
ogee PY aa gg. 2 

where we use ds instead of dA. 


Changing 7 by / and 7 by k, in the first term it may be rewritten 
as: 


da! LT da’ da? _ 
Fk gg2 ds ds 
Multiplying by g”, 
dz! 4 gD, dat da? 
de 9 Fas ds — 


By using the new form of the Christoffel symbol, the equation of 
motion of a particle on a geodesic is: 


don =p dat de 
ds? - 


a =O. 0.7.4 
ode as ( ) 


L _ )kl 
where Di; = 9 Tiny: 
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Geodesic and absolute derivatives 


Having explained the equation of a geodesic using the Lagrange 
equation from variational calculus, next we will derive it from the 
concept of the absolute derivative of a vector developed previously. 

Let us use for s, the length of the curve, time as the parameter. 
The velocity in contravariant form is then 


i dq' 
dt 
Let us call ¢’ the contravariant component of the acceleration 
. dv? 
— ae The absolute derivative of this contravariant vector is 


Dé! = de’ + C/T, ;dq" 


Dividing the differential by dt, we get 


Dé dc? jess dg* 
= Tt 
ae ae 1S as 
, = ee ee OS 
Now returning the original variables qg’ and calling y’ = de 
,_ aq, dg’, dak 
Tae "at at 


If this acceleration is equal to zero, then we have the equation 
for the geodesic; namely, 


agi 
dt? 


5 gt 
dq’ dg" _ 


L~ 


~~ *I dt dt 


The parameter t may be changed to s, and we have the same 
equation that was obtained before using variational calculus. 
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We may say that the velocity v has the same componenet at 
the point q and at its neighbor q+ dq. In other words, v moves 
parallel to itself, preserving the same contravariant components v’ 
and obviously tangent to the geodesic. 


Also, 
Du _ Dv dq* 
ds dg* ds 
where the parameter s represents the length on the geodesic. Since 
Du _ 
ds 
it means that 
Dv dq" - = 


— — =D,v-r* =0 

dq* ds 

Sd ; Du 

where r* is the unit vector along the tangent. In other words, gk 

qd 

is a vector perpendicular to the curve. As we will see later, Weyl 

generalizes the study of parallel displacement to the Riemann space 

and other paths which are not geodesic. And again, this approach 

is equivalent to the absolute derivative of the statement that the 
vector is null. 


Geodesic: Summary 


Let us summarize some of the main features of a geodesic line. 


1. From the variational calculus point of view, the geodesic is 
the curve that represents the minimum distance between two 
points in the space. In a flat space it is just a straight line. 
Einstein associates gravitation with local or proper acceleration 
in empty space, but only if the space is curved may the geodesic 
be generated by a gravitational field. 


In a spherical surface the equator is a geodesic and the mobile 
will move with a velocity tangent to it. 
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2. A parallel transport of a vector is equivalent to having its ab- 


solute derivative equal zero. In this displacement, the vector 
may or may not preserve the values of its contravariant compo- 
nents. If it does, then the transport is done using coordinates 
that we will call geodesic coordinates. If it does not preserve 
the contravariant component, then the variation of the compo- 
nent from one point to another next to it is associated with the 
Christoffel symbol or connecting function. 


. The geodesic equation establishes that a particle in a Riemann 


Space will have absolute acceleration equal to zero. That is, the 
velocity on the path moves parallel to itself and its components 
are kept constant. Thus, if the velocity of the mobile is zero at 
the beginning of the path, it will continue to be zero in its ”par- 
allel” transport. We say that the geodesic is a null one. This is 
a new way to express Newton’s First Law in Riemann Space: 
a particle with a constant velocity or at rest will continue to 
move at the same velocity or remain at rest. The particle is as- 
sumed to be in the gravitational field, but not under any other 
type of field. 


. The geodesic thus obtained has been derived from a variational 


principle, or from the tensor derivative introduced: one by Eu- 
ler and the other by Christoffel. Neither approach makes use of 
any system of coordinates. This independence of the differen- 
tial equation from the system of coordinates is much in agree- 
ment with the Principle of Relativity. Equations in physics 
derivable from a variational principal are called ”self-adjoint.” 
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D’Alembert Principle 


D’Alembert (1771-1783) proposed to look at the dynamic of a point 
as a static case where the forces on the particle are of three kinds: 


1. External forces F’ 


2. Reaction force of the constraint R 


and 


dx 


3. The inertial forces (“may =a): 

Thus, if [ is the vector representing the inertial forces, and F 
represents the external forces F + R, the reaction from the con- 
straint, then m is in perfect static equilibrium. When moving m by 
some arbitrary displacement (which we call virtual displacement), 
the work of [+ F is 0. This virtual displacement is done at a given 
time. They may or may not coincide with the real displacement. For 
instance, if the point moves on a curve c, the virtual displacement 
is A—B. But if the constraint depends on time and c moves to c’ in 
At, the real displacement is A — C’, which is different from A — B. 

As we refer to virtual displacement we are considering small dif- 
ferential displacement 6x, dy, dz, or in generalized coordinates, dq. 

If the displacement is compatible with the constraint, and this is 
frictionless, the virtual work due to the external forces F’ and the 
inertial forces is null. 

For instance, in the x direction 

d?x 
(Fe - maz) be =—0. 

Between the external force and the inertial force there are, from 
the mathematical point of view, great differences: F’ may be derived 
from a potential V(x,y,z), so Fda = —AV; however, the inertial 
forces are not derivable from a function and we may have to treat 
each point of the system individually. 
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In order to make the inertial force derivable from a function, Eu- 
ler (1707-1783), Lagrange (1736-1813), and finally Hamilton (1805- 
1865), transformed the approach of D’Alembert into a variational 
problem. We have seen that Jakob Bernoulli (1657-1705) and his 
brother (1667-1748) tried the problem of finding the path that will 
take the least time for a particle which travels from A to B under 
the action of gravity. 

This problem implies the variation of an integral, and so falls into 
what we call variational calculus. Lagrange changed this problem of 
minimizing an integral into a differential equation. As we saw, the 
differential equation that solves the function f(x, y(x),y/(x)) such 
that [7 f(a, y(), y'(x))dx is a minimum (or maximum) is 


Of dof 
dy dx dy” 
where 
d Of —(Of\Oy' | (Of\ oy | Of 
dx Oy! (5) Ox | (Si) ae Oy!’ 
since 


5a = alesule) ve) 
In order to bring in the variational aspect of the dynamic equa- 
tion, we integrate between two times t,; and fz, and the virtual dif- 
ferential works as expressed by D’Alembert: 


If F' is derivable from a potential V 6V = —F'dz, then the first 
term of the integral is: 


te to 
i dio (vas 


ty ty 
The second term may be integrated by part: 
i) to 


Bid po ; od 
a gy tox di: = [so(ma)] = mé_ on des 
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Taking 6x2 = 0 at t; and to, the first term is zero, and the second, 
noting that the derivative of the variation dx is equal to the variation 


of the derivative, is: < (60) = 52) = 04. 


te ta 
7 ma(ea) at =— ["5(™=) ar = sf T(t) dt, 
ti ty 2 


where T is the kinetic energy. So: 
r= sf" -~V +T) 


We indicate the derivative of any variable with respect to t by 
an upper dot (-). By calling £ = (1 — V) (the excess of kinetic 
energy over potential energy) the Lagrange function, the problem is 
reduced to finding the £ that make he £ dt a minimum or stationary 
between two instants t; and ty. T is called the action integral. Using 
the Euler differential equations, we find the function x(t) from 


ac dae _ 
Ox dt0t 

If we are working with generalized coordinates q’ , q?...q", then 
ddL OL 
dt og — ag 0, (0.7.5) 


and these are n equations (i = 1...n). 
The equation is called the Euler-Lagrange equation. 
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Problems Using the D’Alembert Principle and the Euler- 
Lagrange Equation 


Four masses slide down on the surface of the triangle without 
friction. They are connected by a massless thread. Assume 
they slide positively from B to A through O. The coordinates 
will be called g, which is the distance from O to the point 
masses. The virtual work will be due to a variation dq which 
in this case is similar to a real one (a slide down the surface). 


The work of the inertial force is: 


2 2 2q 2 
dg — M3—504¢ — Ma, 09 


while the work of the external force (gravity) is 
myg sin adg + mag sin adg — m3 sin Bdg — m4 sin Bdg 


We assume that the motion is direction downwards for m,; and 
mg, as are the forces acting on these objects. For m3 and 
m4, the force is opposite in direction to the displacement, so 
according to the D’Alembert Principal, the sum of these two 
works is zero. That is, 

Oq 
Ot? 
Only if the left side of the equation is null will the motion be 


gsina(m,+ mg) —gsin B(m3+ma4) = (My +m2+m3 +m) 


constant with uniform acceleration. 
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Now let us solve the same problem using the Euler-Lagrange 
Equation, £ = T — V, where T' = i> mq?. Assuming zero 
potential at the base of the triangle ABO and taking h as the 
height, 


Vi=(h-qsina)mg , Vo=(h-@gsina)mag , V3 = (h+q3sinB)msg 


VY, = (h+ qsina)mag | 


Now 
Mies —(sina)m deme —(sina)m 
Oa 19; go 29 
Fe = (sin B)anag St = (sin B)mag 


Since dq; = dqz = —dq3 = —dqa, 
= = -=- =g{(m + m2) sina — (mg + m4) sin §] 


and 


OL 
where M = ¥>m,. Thus, the Euler-Lagrange Equation es 
q 


M@g = gl(m, + ma) sin a — (m3 + m4) sin J] 
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The point MM moves on a horizontal circumference, varying the 
angle @. The circumference rotates at the same time around 
a point O fixed to it with angular velocity w. There are two 
variables: g' = 6 and q? = wt. We assume that w is constant. 
So, 


g =6 and @=w 


Since M moves on a horizontal plane, there is no potential 


variation. Now, T = g(t + y*) where m is the mass of M. 
Also, 


z£=OM-6&= (OC + CM) -€; = Rcos(wt) + Rcos(@ + wt) 
where R is the radius of the circle, and 


y ZOMe = (OC +OM)-& = Rsin(wt) + Rsin(wt + 8) 
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Thus, 


é = —Rwsin(wt) — Rsin(wt + 0)(w + 6) 
y = Rw cos(wt) + Rcos(wt + 0)(w +8) 


f= 5m( Feu? + R?(w + 0)? + 2R?w(w + 6)[sin(wt) sin(wt + 0) + cos(wt) cos(wt + 4)]) 
_ 1 ne 2 42 5; 
T= 5mk [w~ + (w+ 0)° + 2w(w + A) cos 6] ea el) 


Now, since v is a constant, the Euler-Lagrange Equation is 


OE: Oly SOs OL =. 
dq dtdq Oq dtdg — 


0 


There are two equations, one for g = @ and the other for g = wt: 


1 ; ; 
a = sR (—2u(w + 0) sin 0) = —mR?w(w + @) sind 
d OT = 2 d * = 277 ns 
aon lw + 8) + wcos 6] = mR (0 — w6 sin 0) 
Ol. aor a & - ae 
ke intial ceca _ _6 = 
a0 dt Ob w* sin 8 — wé sin 0 + wé sind = 0 


or 
6 = —w’ sind 


This is the equation of a pendulum as seen from the center of 
the circle. 
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EKuler—Lagrange—Poisson—Jacobi-Hamilton— 
Canonical Equations 


The Euler-Lagrange equation is: 


dor Or 
dtOq, Oq 


=Q, , 


dq 
where Q, is the force and q, = aes 


CONTENTS 


Assume that q1, G2, dv, g: are the degrees of freedom of the system 


of n points where there are h equations of constraint such that: 


k=3n—-h. 


In the following we simplify our equations by eliminating the 


subscript v and using such g as a general variable. With the same 


ead : a ; sath hee aie OE. 
criteria we will use a comma to indicate a partial derivative like i 
qd 


as Tq. 


The kinetic energy T is given by T = a’”"q,q,, where the ma- 


trix a”” is a function of the coordinates q, q2, .... Assuming JT’ = 


T(q,q), its variation is given by: 
of = 1.9 69+ T,¢-00 | 
By defining p as: 
p= 0G ’ 
oT =p oqg+T,¢ oq . 
By summing and subtracting q dp, then: 
6(qp —T) = 5K =q4 dp—T,q 4g, 
where we have defined K as: 


K=qp-T 
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If K is considered a function of g and p only: 


OK _ OT 

Og (OG 
and 

OR «4 

one 


From the Euler-Lagrange equation written above and the defini- 
tion of K, 


dp OK 
at aq 
or aK 
Dae Og =Q 
For @ = 0, 
2, BOIS 
If Q is derivable from a potential V such that 
iN. 
=a 
where V = V(q,t), then by defining H = (kK + V), 
oH 
Op i 
and 


OH OK OV 0 
dq oq | aq ee et 


By using the Euler-Lagrange equation and the definition of p, 
oT 


P= 8q’ 


oH _ 
oq 


These two equations, where p and qg play symmetric roles, are called 


the Canonical Equations, or Hamilton Equations: 
OH. OH : 
De =q and on Spin 

The variables are qi, g2,.--, Ge and pi, p2,..., Dk- 
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Lagrange—Poisson—Jacobi—Maupertuis 


After the Lagrangian £, other interesting functions were created 
whose first derivatives will give rise to the equation of motion. 
They assume different sets of variables. The generalized variables 
1, 42,;---,t that we represent by (q,t), were used by Lagrange in his 
function £(q,t) and Jacobi in his function S(q, t). 

As we have seen, £ is (T’— V), the excess of kinetic energy over 
the potential, and it holds the equation 


dol OL 
dt0g Oq — 


The function S' of Jacobi is given by 


56H - [cat de 


If £ = (T —V) is written as £ = —(T + V — pq), where pq 
indicates p,q; = 2T (repeated indices implies summation), then: 


stat) = [UP +V) - pal at 


By saying that H = T’+V, we then rewrite the previous equation 
as 


t q 
s(a.t)= [ Hat— [pag 
0 qo 
which generates two first order fundamental derivatives: 


HOD goat VOR et. 
oe Se ge 


H 
These derivatives resemble the Schrodinger recipe for quantization: 


0 0 
aan and any 


The Jacobi Equation that will give q = q(t) is 


2 = H(ap.1) where pees 
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Thus, the Jacobi Equation reads: 


as. AS | 


HG — ao }= Ot a 


0 (0.7.6) 


where S(q,t) will be equivalent to solve the particle motion. 

To show this, let us see that S(q,t) generates the Canonical 
Equations. In fact, a solution of the Jacobi equation will give 
S = S(q,t,a) where a is a constant of integration (one for each 
coordinate of q). By deriving S with respect to a we generate a new 
constant b such that 


a8 _ 
da 
In S(q,t,a) we have q = q(t) for each a. Then from 


b (0.7.7) 


Os 
~—=- 0.7.8 
pg ee (0.7.8) 
we get p(q,t). 

Let us differentiate (0.7.7) with respect to t and a. From (0.7.7), 
since b is a constant, we have 


Os... O78... 


plea abo. 


By differentiating (0.7.6) with respect to a, which is implicit in 


ee and a 
ot Og ’ 
PS OH PS _ 
dadt ~=,OS, Aadq 
oe 
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But we saw that 


sO 


which is one of the Canonical Equations. 
To obtain the second Canonical Equation, we differentiate (0.7.8) 
with respect to t and (0.7.6) with respect to g and compare: 


OS O'S 
Os 
as an oq) an _, 
otog OS. &q | dq 


Equating, after a change of sign in the last equation, the second 
term gives rise to the previous Canonical Equation and the third 
gives 


OH 
oq 
which is the second Canonical Equation. 
In this form, we have shown that the Jacobi Equation creates 
S(q,t) that solves the dynamical problem in question. 
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Maxwell Equations and Their Invariance Under Gauge Change gy 


The Maxwell Equations are based on the following four fundamental 
equations: 


1. Faraday: 
OB 
Vx #4+—=0 
+ Bt 
2. Gauss: 
pa 
€0 
3. Ampere: 
VxBe=J 
4. No free magnetic pole: 
V:-B=0 
where é i 
H=— , D=eok , /poto=— 
Ho c 


Since V- Vx = 0, Ampere would say that V- J = 0. But the 
OE 
continuity of charge, V- J+ Bem 0, must be fulfilled. So Maxwell 


changed Ampere’s equation to the Ampere-Maxwell Equation: 


vx Bp 2E 
C 


= pod 
Ot Lo 


Now, if we apply V- to both sides we generate the continuity equa- 


tion. 
So the Maxwell Equations are then: 


Vx B=-— (0.7.10) 


Va (0.7.11) 
€0 
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E 
V x B= pod + g (0.7.12) 


c20t 


V-B=0 (0.7.13) 


These four equations can be simplified by introducing the EM 
Potentials: A and ®: 


B=VxA (0.7.14) 
OA 
&=-VO—-— — ata 
V y (0.7.15) 
From (0.7.13) and (0.9.51), 
V-VxA=0 
From (0.7.15) and (0.7.11), 
A 
Vib= yew A jesyes a=? 
ot € 
or 
2 0 p 
Vo4+2(V-A)=—-? (0.7.16) 
Ot € 
Substituting (0.9.51) and (0.7.15) in (0.7.10), 
A 
V x (-V®—- = Sv x A — all cancel since V x V® = 0 


Substituting (0.9.51) and (0.7.15) in (0.7.12), 


A 
VxVx A= pot + 2 (v6 _ ae (0.7.17) 
ct 
Observe that 
VxVxA=V(V-A)-V°A 


or 


a Vort a) = Sig. (0.7.18) 
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(0.7.16) and (0.7.18) ar coupled equations for A and ®. 
A and ® may be changed without varying Maxwell’s Equations 
(0.7.16) and (0.7.18): 


A+A'=A+VA 


and ay 
&® > 6’ = 6 — — 
st Ot 


where 4 is an arbitrary differential function. These gauge transfor- 
mations leave (0.9.51) and (0.7.15) unchanged, and thus B and E 
are invariant. In four-dimensional tensor notation, these transfor- 
mations may written as A“ > A” — OFX. 

Lorentz fixed the gauge by requiring that 


or in four-dimensional tensor notation, 0,A” = 0. This choice de- 
couples (0.7.16) and (0.7.18) to the following two Poisson Equations 
in four-dimensional space: 


®=-—p and A=-J 


or 
OPA 
=— d WA- = =-poJ . 
p an V Ct Lo 
This is the so-called Coulomb Gauge, given by the condition that 
V-A=0. 


Tensor Form of the Maxwell Equations 


From the potential equations ((0.9.51) and (0.7.15)), an antisym- 


metric tensor F,,, is generated: assuming the coordinates x, : («° = 


te Sa = ae = 2 and 001,01 = 0, Or = Ops =O), 


Fac Oy Ay = OpAm 
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where y and y run over the indices 0, 1, 2, 3 that represent the 
four dimensional space with the tensor potential A,: Ag = ® and 
A, = A;, Ag = Ay, and A3 = A,. Thus, F,,, is given according 
to (0.9.51) and (0.7.15) by the covariant tensor F,,,, which may be 
expressed in its contravariant components by means of the metric 


tensor 
1 0 O 0 
en OCeE: e <0 
cad) Div ey es 
0 0 O -1 
As was shown in Lecture 2, 
Gg 95=1 


and thus 


It results as 


0 —-E' —-E? —F° 
E' 0 Be —B? 
E? —-Be 0 B!} 
Ee B? —B 0 


Fu = 


Weyl deduced this tensor from the connecting function ®,,, such that 
dl = ®,,dq", which has a non-zero rotational. So the circulation of 
the scalar 1 around a closed curve is non-zero. Thus, 


Fw ~~ (Puy = ®,,.) 
has the following property: 
Oy Opt ay Onl ya Ol 


From the two-index tensor F“”, we may generate a one-index tensor 
(or a vector) as 


jE = 4a” 
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where J“ is the four-dimensional current J“ = (p,j). Thus, by 
identifying ®, with A, we regain the Maxwell Equations. 

As we will see, Weyl makes a connection between matter (not 
gravitation as said in an earlier work) and the EM fields which move 
the particle out of the gravitational geodesic. 

Let us assume a potential A, and shift it via the gradient of a 


function WU: 
Ay(z) > A(x) = A, (x) + 0,V(z) . 


Since the relevant property of A, that generates the neutral fields E 
and B, is its circulation, the gradient does not make any contribution 
and leaves the fields EF and B invariant. 
This transformation of the vector potential A,, which is a non- 
observable quantity, may be seen as a gauge transformation. (Lorentz, 
in his famous work on EM, assumes that 0, U(x) = 0). That is, a 
change in the connecting function ®,, or a gauge transformation, 
does not change the Maxwell Equations. In other words, the tensor 
F has the same form as functions of the field E and B. 
The tensor Ff” gives rise to some invariants: 


Invariant (1): Fy, 


Invariant (2) : Pata Oe 


where e*’™ is the four-rank antisymmetric unit tensor. Invariant 


(1) gives rise to: 2(B? — E”) and Invariant (2) to: (E- B). 
Note: Calling 2 the four-potential vector (®o, A;, Ay, Az), the 
Maxwell Equations may be reduced to two equations: 


Q = —pol 


and 


V0 S06 


where 


PS Cp; dedi ds 
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Review of Relativity 


1. Concept of Simultaneity: Transmission of Information 
Between Clocks 


Let us carry out the following experiment: A and B are two 
walls and a ball goes from A to B with velocity v. 


ball 


At time to, the ball left A, and at time t,, according to B, 
reached B. As the ball reaches B there is a sound that reflect 
back to A. The sound travels with velocity c and reaches A at 
time ty (according to A). 


Time tz must depend on the velocity as 


ad ad Ss Sk 
tg —to) = ae 
(t2 = to) c ov ‘ 2 
If c >> v, then 
tg —to & - 
v 


But if v is comparable with c, then the problem of time will 
v 
depend on c in the ration —. In other words, the sound is the 


Cc 
element to determine the simultaneity of an event observed by 
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A and also by B. If B moves away from A at the same time 
that the experiment is carried out, the sound speed for A is 
different from the one measured by an observer at B. 


Einstein took light as the element for the correlation of time, 
and postulated that this velocity is the same for A as it is for 
B: a radical new concept. 


2. Simple Derivation of the Lorentz-Einstein Transforma- 


tion 
A A 
S(a, Y, zy t) S(E, 1, C; 7) 
O > oO! ° > 


The system S’ moves with respect to S with a constant veloc- 
ity v. Let us take a point P in the space of S’ with coordi- 
nates (€,0,0,7) and assume that S finds P to have coordinates 
(x,0,0,t). Our next assumptions are 


(a) The transformation from the S coordinates to the S’ coor- 
dinates is a linear one. 


(b) At the limit of low velocity v, the transformation is equal 
to the Galilean one. 
(c) The space is isotropic. 


(d) The velocity of light is equal for both observers, even when 
the light was emitted when O and O’ were together and 
the light was moving away from O’. 
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Thus, from (a), assuming a, b, a, and £ are constants, 
€=ar+bt and r=ar+$t 
From (b), b = —aw, or 
€ = a(z — vt) (0.7.19) 
From (c), we may exchange z for €, t for 7, and v for —v: 
x =a(€+uvT) (0.7.20) 
From (d), we may write 
geek and! 6 er (0.7.21) 


From (0.7.19), (0.7.20), and (0.7.21), we may determine the 
value of a: 


cT = €=a(ct—vt) =at(c—v) and ct=2 =a(cr+v7T) = ar(H) 
By multiplying these two equations, we obtain 
(tr) = a*(rt)(c— v)(c+ 0) 


or 
3 eo 1 if 

a = = 5} — 2 

c — Vv ieee 1-8 


2 


( 
where 6 = as (In the usual notation, a is called 7). 
c 


By using (0.7.19) and (0.7.20) again, we obtain 


1 
&= 1p — vt) (0.7.22) 


1 
i Vio-e + vr) (0.7.23) 


and 
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From these equations, we solve for 7 by substituting € into the 
second equation, 


So 
= | = 90 + ot] = (yaa 
or 
UL 
os ft - -) (0.7.24) 
By changing t for 7, x for €, and v for —v, 
= a(r + =) (0.7.25) 
C 


Equations (0.7.22), (0.7.23), (0.7.24) and (0.7.25) express the 
Lorentz-Einstein Transformation. In matrix form, we write 


or 


The product of the two matrices, is, of course, equal to the unit 
matrix. 


3. Variation of Length and Time from the 5S’ Frame to the 
S Frame 


From small increments of the coordinates, we may write the 
following equation: 


Aé = y(Az — vAt) 


When the variation of the length is taken at a given time, At = 


0 and 
_ AS 


Y 


Ar 
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mince 7 S.No Ae, 


We say that, as A€ moves away, its length looks contracted 
from the point of view of S. Similarly, as 


MS (Ar + “ AC) 
C 


at one point in the space, the clock at S’ moves slowly. That 
15,68 Ae =), 
At =7yAr 


and so At > Ar. Thus, there is a time dilation from Ar = At 
when the two clocks were at rest to a At larger than Ar as O' 
moves away. In other words, the time interval measured by the 
moving clock will be small compared with that of the standard 
clock at rest, in spite of the fact that they were synchronized 
when O and O’ were at rest. The two phenomena observed 
at S, time dilation and length contraction, are just velocity 
effects which become more significant when v approaches c (the 
velocity of light) in magnitude. 


. The Differential Length of the Four-Dimensional Space 


as an Invariant 


Let us write ds = drc, where 7 is the ”proper time” and c the 
speed of light. As we have shown, 


v 
T=7 aaa 


We assume that the events that we measure are located in the 
moving frame where 7 and € are coordinates. We see that 


an(i- 28) o(t-8)=3-! 
a ed) Ce) 7 4 


dr v? 5 3. de 
a = 1-— 2 = dr = dt = D 
cdr? = (edt)? — dz? = gijdxr'da? 


— 
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where 7,7 = 0,1, dx°® = cdt, and dz! = dz. The 2 x 2 matrix 


in Gij is 
1 0 
ames es 


As we see, ds is not the module of the vector (cdt, dz) under 
the normal scalar product, which will correspond to a g;; equal 


(on) 


We may see it, however, as the module of the vector (cdt, idx) 
with an imaginary x component. 


to the unit matrix 
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Relativity According to Minkowsky 
The principals of Minkowsky Space are: 


1. Any object or ”substance” in space and time is a point in the 
4D world. The evolution of that point describes a world-line 
in the space. Objects that are stationary or at rest in the 
(x,y, z,t) coordinate system have a world line parallel to Ot. 
Objects with constant velocity v have a world-line forming a 


dx 
constant angle with Ot, such that — = v indicates the slope 
of a straight line. A world-line parallel to Ox or perpendicular 
to Ot indicate an instant view of the world. 


2. A very important scalar is the magnitude of a vector, which 
represents the distance between events. In orthogonal coordi- 
nates, this is expressed as 


et = to)? a Ca Dy) aes Yo)” ae ee zo) = he (0.7.26) 


The value of k? is independent of the coordinate system that 
may be changed by any linear transformation, or equivalently 
by any rotation of the 4D system. With k? = 0, the scalar 
represents the transmission of light with velocity c from Pp to 
P or vice versa. In general, the scalar product of a vector in 
the four dimensions is 


=>. => 


WW =2'y' Qi; (0.7.27) 


Comparing (0.7.26) and (0.7.27), the matrix g;; for orthogonal 
coordinates is 


i Oe. fe 
Got 20 -20 
22" eps Shy 36 
0. O° 0 1 


In other systems of reference (for instance, in spherical coordi- 
nates), gi; depends on the coordinates. 
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A point in the space-time P(x, y,z,t) may describe a line that 
we call a ”world-line.” As it moves from P to P+ 6P (a point near 
P) with velocity v < c, the value of the differential distance in the 
space-time is (assuming two-dimensions only) 


ds? = c?dt? — dx? 


If ds = 0, the substance described by P(x, t) isa photon moving with 


dy 
velocity c; otherwise, the substance moves with velocity v = — <c. 
The increment ds is greater than zero and invariant with respect to 
a change in the coordinates, say t > t’, or x > a’. 


”Proper time” 7 is given by the following equation: 
cdr? = ds? = c’dt® — dx” 
For stationary points, dr = 0 and 
ec Re cea 


That is, the proper time is the same as t. If the substance is at rest 
in the system (t’, x’) that moves with velocity v, the proper time is 
i 

The invariant ds may be seen as the distance in space-time (clock 
and ruler) of two close events. and if the two events are simultaneous 
in the reference system (x,t), ds? = 0 and it will be so in any other 
system related to (x,t) by linear transformation. 


When the system of coordinates is not orthogonal, then 
ds”? = qijdx'dx! 


where gj; is a function of the coordinates, 7,7 = 1, 4. 
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Plges) 


The point P is a point in the space-time coordinates. If a photon 
leaves O and reaches P in the direction OP, it will be at P at the 
(x,t) space-time point. The coordinates of P are such that 


(c¢)?-— 2° =0 


A mobile with velocity v < c follows the same world trajectory OP. 
The coordinates at P will fulfill an equation such that OP-OP =i 
or 


(f= =a 


where a is a constant greater than zero. Taking a = 1 and allowing 
P to move around O, OP will describe all the objects that leave 


d 
O and reach P with velocity 7 = v. The point P describes a 
hyperbola c?t? = 27+ 1. If we choose OP aa the divection of the 


new time coordinates, then those objects will be at rest in the new 
system with time coordinate 1’. 
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ct 


Once we choose the axis Ot’ we choose an orthogonal axis, O2’, 
such that the inner product OP -OP’ = 0. The dot product in this 
case will have the following g matrix: 


(0 4) 


tt’ — ar’ =0 


==> 
Since OP = (ct, x), then OP’ = (x, V14+ 22). 

It is simple to see that the tangent to the hyperbola at P is also 
perpendicular to OP. In fact, the tangent at P has a slope a such 
that 


sO 


t 
peepee 2 


dx JV1+ 22 


because ct = /1+ x2. Then 
1 J1+2? 


cos a = = 
V1l+tan?a V1+ 22? 


and 
tana x 


V1 + tan? a 7 V1 + 22? 


sina = 


. . = =? 
may be taken as the components of the direction OP’ along Oz’. 


Aaa visa x 
OP’ = (cosa,sina) = ( : ) 
V1 + 2x?" V1 + 2x? 
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5 
And thus, OP’, which is parallel to the tangent, is perpendicular to 
OP, which was taken as the time ¢’ axis. Explicitly, OP = (ct,z) 


and 
xcosa = V1+ 2x?%sina 


or 


—aVl+n? alte? | 


OP'.OP _ = 
V1l+2272 V14+2r? 
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Lorentz Length Contraction From the Geometry of World-Lines 


ct 


Let us return to the hyperbola in orthogonal coordinates lying in 
the (t,x) plane and given by the equation c?t? — x? = 1 or 


1 
#=+-V1+4+ 2? 
c 
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The curve consists of two different hyperbolas, one for t > 0 and the 
other for t < 0. They are both symmetric with respect to the origin 
O. The asymptote passes through B(1,1) and OD is taken as 1. In 
fact, the asymptotes have an equation given by 
dt 
ree =+1 


Cc 
dx V14 7 4 400 


Thus, the equation of the asymptote is 


1 
t= +-2 
é 
and the so the point B is on the asymptote. 
The tangents of the angles a, 3, and ¥y are such that 
tana = co HAL ly tany =< 
c v 
tany = ” is the slope of the straight line Ot’. By choosing the new 
@ 
time coordinate of Ot’ passing through A’, the point A’ will be at 


rest in this system of (2’,t’). Next, we choose Oz’ parallel to the 
tangent at A’ that has the slope 


(=).- Gee) 
dx A’ 7 1 + v2 Al 
dt 
As A’ goes to A on the world curve, (=) is a horizontal line 
A! 


t 
parallel to Ox. Also, as c > ov, (=) is zero, so in both curves 
x) a 
dt ) 1 dt 


—} =-,orc—= 
dx c dx 
1, which is the slope of the asymptote (the angle £). 
The rhomboid A’B’D'O becomes OABD as A’ becomes A, B’ 
becomes B, and D’ becomes D. OD" has the same length as OC’ as 
/ 


Ox’ becomes identical to Ox. As x4 > 00, 


seen by the observer at the (x,t) system. The ratio indicates 
the difference of length of a ruler OD as (2’, t’) moves with velocity 


v. So the same ruler is seen by the observer at (x,t) once at rest 
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(OD) and then moving with velocity v as OD’. If the ruler is at 


1 
rest on Oz’, its world-line is parallel to Ot’ and its slope is —. So, 
v 


1 d 
tana = —. Since = = v, v is the tangent of the angle between Ot 
v 
and Ot’. 
To find the ratio between OD’ and OD, we use the following 


analytic procedure, our first goal being to find the coordinates of 
/ 


tana = (c<*) = (=) 
dx / a V1+2? Al 


tar 1 
Since A’ is on OF, nee —, and since it is also on the hyperbola 
x Al U 
Ct, = sea + 1. 
1 
La = ——— and ty = 
2 2 
care 1 
ve Ve @ 


Substituting x4 into the equation for tana we find that 


Vv 
tana = — 
C 


1 
The point B’ is the intersection of the asymptote tg, = —xg and the 
c 


eee 
line D’B’, whose slope is tany = © So the value of the coordinates 
v 


at B’ are 1 
tp = —(xpr = OD") 
U 


and so we get 


or 


— 

We see that as v > c, OD! > 0. 
Also, since B’ is on the tangent at A’, whose slope is tana = —: 

C 


c(t — tar) = u(xpr = LA’) 
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Substituting tp, ta, and x4 into the equation with the values pre- 
viously found, 


Uv 
LB = c 
st 
cC 
Thus, 
2 
Uv 
/ 
OD = /1—- 2 
and since OD was taken as 1, 
OD | vu? 
OD — Ce 


This means that as the ruler moves with velocity v on the Ox axis, 
2 
v 
it is seen as shorter by a factor of 4/1 — —. This was first discovered 
c 
by Lorentz while he was working on Electron Theory, although he 
assumed a unique and absolute time for every event in the world 
and for every observer. Einstein was the first to take time as a new 
variable: a different time-clock for every observer. 
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Proper Time as the Parameter to Create Velocity and Accelera- 
tion in the Space-Time System of Coordinates 


to =cta, 


Since x” — a = 1, 
dx® da! 
27° — 27'— = 0 
“ds ds 
dx® dz} 
Calling 7 = and Te ae ra 
oe? — ae! . 


Thus, (x°, x1) is orthogonal to («°,z1) according to the new in- 
ner product discussed above. Taking the derivative of the previous 
equation, 


~02 a .412 ss 
e° + ¢°r° — ah — Ze’ =0 


or 
ce ee - 02 .72 
#o° —#'a' +27 —¢' =0 
oe arg Se =e: 
~—Q2 -12 * : < 
Let us call «2° — # =v? the magnitude of the time-space velocity 
squared. In this simple case, it is (c? — v). The vector acceleration 


122 CONTENTS 


in this space is (#°, 71) — A such that 


Ce ears 


Taking the magnitude of these vectors, the acceleration is 


v? vu? 


Ae eos = Se, 
eee jn°,e"| pp 


where p is the magnitude of the vector OP or the radius of curvature. 
For any world-line passing through P, this is a hyperbola with 
the same tangent at P (or with 3 points in common). The curvature 
of that hyperbola is then equal to that of the world-line at P. 
If the acceleration (#°, #1) is zero, then the radius of curvature is 
oo and the world-curvature is a straight line. 
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Relativistic Dynamics: Application of Tensor Theory to 
Special Relativity 


Any meaningful quantity in the theory of relativity must vary covari- 
antly under coordinate transformation in a geometry generated by 
the inner product consistent with the Lorentz transformation. The 
fundamental tensor of such a geometry, g,,,, in the four-dimensional 


space x” (where wy = 0,1,2,3 and 2° =cl.a! = 2,27 S42" =z in 
rectilinear coordinates) is such that 


ds? = Oye da 


is Invariant. 
Again for rectilinear coordinates, 


ds SCR ade"? dy? = d2" 


or 
1 O 0 0 
0 -l1 O 0 
Iuv = 
0 0 -1 O 
0 O 0 -l 


To make our discussion simple, let us work in rectilinear coordinates 
and two dimensions only, so z° = ct and x! = x. Thus, if the velocity 


dx 
of the particle aoe is constant, then 
ds* = (c? — v”)dt? 


or 


ds = Vc? — v2dt 


After integrating and squaring again with s(t = 0) = 0, we have 
sis, geese 
Differentiating with respect to t, we find that 


O-e 0 
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where z is the time-space position vector of the point M in study, 
and v is the velocity in the O frame where the observer measured 
t. Thus, in general we may say that the four vector position and 
its four vector velocity in the S frame are orthogonal, by the inner 
product determined by gy. 

By taking the second derivative with respect to t of s? = go — gl? 


we find that 


C2 


where @ and v are the acceleration and velocity of the point M in 
the O frame, where the observer measures the time t. Thus, 


@=(#,#') and @=(2°,2') 


where the dot means the derivative with respect to t. 

We observe that the vectors @ and v depend on the frame of refer- 
ence. However, the following definitions of velocity and acceleration 
are tensors that depend on the particle itself only. 


3: ¢" as a function 


Given the four-vector displacement x°, x, x”, x 
of the proper time 7 (which is measured by the clock at the particle 
in question), the following derivatives form a vector that we call the 
velocity of the particle: 

_ daxt 
~ dr 


In the space-time, x“(7) describes a curve with a differential length 


U" 


given by ds = cdr where c is the velocity of light. If we wish to relate 
U¥ to the three component velocity w determined by the variation 
of x', x7, and x° with respect to t, we must remember that there is 
a time dilation such that 


dt = y(u)dr 
Thus, 


dat _ dxdt i bf 
de Sipe 
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where i = 1,2,3 and y(u) = 


We have shown that - = (G,.14),80 
U" = y(u)(e, u) 
and the covariant component of U is 
Uy = y(u)(c, -%) 
The magnitude of U? is 
U? =U4U, =P (u(e —w) = 


Thus, U? is invariant. This was not the case of v? previously dis- 
cussed. 
The corresponding acceleration is defined as 


dU" 
Bia coe 
dr 
Since u also depends on 7, 
d AE is ty. 0 24 
— 7)| = sae hs 
AY = Sh(u\(@] = we, Dat ra 
where : ; 
oT eG 8g d 
d= Fat 2,2") qe 6 92) 


which is the three-vector acceleration measured by the observer 
whose clock reads time t. 
The covariant A,, is equal to 


_ y ch Geer cc eae 
Ay = SwA = (u) [eo (pita) 


The expression of A“ may be simplified in the following cases: 


126 


1, 


CONTENTS 


If u = 0, ie., we are measuring A in the instantaneous rest 
frame, then 


=] oo 
(0) and a 0 


so 
A* =(0,a) and A, = A”’g,, = (0,—a) 


The proper acceleration vector a = a, and the scalar A? is 


A? = A,A* = —a? 


. If the magnitude of w does not vary with time, such as in uni- 


2 
(ale then oa = 0 and 
r dt 


AY = 4(u)(0, Ya) 


form circular motion where a = 


and the scalar A? = A“A,, = y?(u)(0, —7a7), so 


2ha4|2 
Perea | 
r 
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Relativistic Forces 


The four-vector force is defined as 


d 


_— 
ae 


MoU) 


where mo is the mass of the particle in its own frame. Let us relate 
F to the classic Newtonian force f: 


P= (u) = 


on q 7 d 2 
yer movi) = yu) ¢ (me, mu) = y(u)( cdt ’ dt! ) 


= 109 SG at) = 10 Sat) 


Using tensor notation, we may find the invariant J = F*U"g,,, 


where 
P= (PPP, PS) = oy) ( FPP) 
dx 
and 
UY = (U°,U',U*, 0%) = y(u)(c, @) . 
Then, 


1dE 


[=P U = FU" Gy = 9(a) (2, F) WeDo 


= w)(S-Fa) 


This invariant J, at the rest frame where m = mo, y(u) = 1, and 
(60), is 


fig d(mo) @ (=) 
dr 


since dt = dr for the rest frame. Equating both Is, 
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A pure force means that there was no change in the rest mass mo. 


Thus, for a pure force (cme = 0), then J = 0 and 


and so [ = F'.-U = 0. The force F’, when it is a pure force, can be 


Faw (£27) , 


Cc 


written as 
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The following examples may prove helpful in understanding the 
difference between the three-vector and the four-vector velocity and 
acceleration: 

Let us look at a high-velocity particle in an accelerator where the 
path is a circular one with radius of curvature R at the moment when 
we study the state of the particle. If the instantaneous three-vector 
velocity is along the path, then there will be a radial acceleration 
equal to “ that is normal to the velocity. All of this is from the 
reference frame of the observer at the accelerator. We then ask 
what the acceleration is for the observer at the particle, that is, the 
proper acceleration A” in the four time-space coordinates, and a, 
the three-vector proper acceleration. 


A? = AMA, = 0? . 


w of) aea()). 


d 
If we assume that |u| is constant in time, then 1 (which implies 


du 


—) will be zero, and 


dt 
At = ¥(u)|[0, ya] 
and 
A? = AMA, = —y*(u)a? . 
So 


a=7(uja. 


For example, if wu = 0.888c, then 


a 5 x (0.888)? =, = 4.2, 


~ 1 — (0.888) 
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It is common to express the Lorentz transformation using the hy- 
perbolic functions. For instance, if we define the angle ® (sometimes 
called the rapidity) as 


tanh 6 = is 
é 
where u is the magnitude of w, then, 
d ldu a 
—(tanh 6) = —— = — 
a ae 


But the first member is 


1 \d® d® u\d® 1d® 
a pee ho) = (1- ) 7 
(—) dt testes ) dt c2/ dt +? dt 


If the frame (7,&,7,¢) sees a particle in rectilinear motion, then its 
/ 


U u 
acceleration a as seen by S,, (—-), is related to a = —— by a= a, 


z 
where the primes refer to the coordinates in the moving frame. Thus, 


1d® al 
ydt ci 
and 
ae 
vat 
1 
but ae so 
ae 
dr cc’ 


If a is constant with adequate boundary conditions, then 


TO 
CS 4 
c 


In this particular case of a being constant, we find interesting rela- 
tionships: 


1. Since tanh ® = oe then 
C 


a tanh (=) ’ 
Cc Cc 
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a il TO 
y(u) = ——$———— = cosh ® = cosh (=) 


Vv1l— tanh? ® 


3. In a previous discussion, we found that 


di’ d - 
er qy (14) 


which is also true for the magnitude of @ and wt’, so we may 


a= 


write q 
a= S(r(u)u) 


and since a is a constant, 
ta = y(u)u (0.7.28) 


with an adequate boundary condition of t = 0. Thus, using the 
previous derivation for y(u), 


a 7(u)— = cosh(®) tanh ® = sinh ® = sinh (=) ; 
c C 


This gives an equation that relates the times t and 7 and the 


value of a: be ee 
— =sinh (=) (0.7.29) 
Cc Cc 
4. Since u = then 
dt 
ch =u=ctanh®. 
dt 


Integrating and assuming some 29 when t = 0 (also assuming 
that =O torty=0), 


t T 
tanh(®) dt = c| ytanh ® dr . 
0 


t=0 


| de =(¢=75) = c 

a0) 

From the expression of y previously derived, 
@ p72 


z—ay=c [ sinh® dr = | © sinh ® d®& 
0 


Bo 
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or 
2 ®& 


L—-—Xp= < | cosh 
a ®o 


Since 7 = 0 for t = 0, ®) = 0. Then, 


2 
X=2X+—cosh®. 
ay 


Let us define another property of ®, the concept of ”rapidity” of 
a particle moving with respect to O; with velocity 


where wu’ is the relative velocity of Oj that moves with respect to 
O; with velocity v. wphe “rapidity” of ® is such that its hyperbolic 
tangent is equal to —: 

c 


tanh 6 = _ 
C 
If the particle is at the observer Og, u’ = 0 and u = v, thus 


tanh ® = e , 
Cc 


Assume now that there are three observers: O;,O2,O3, and three 
particles located at the observers’ positions. If v, is the relative 
velocity of O2 with respect to O,, then the ”rapidity” of the particle 
at Og with respect to O, is ® such that 


C+Uy, 
ae 
C— VU, 


This equation is easily derivable using the expression for sinh ® and 
cosh ®, and by recalling that 


sinh ® + cosh ® =e” . 
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133 
If W is the rapidity of the particle at O3 with respect to the observer 
at Oo, then 


C+ Ug 
ev = ,/—_= 
C— V2 


where v2 is the velocity of the particle at O3 with respect to that at 
Op. 


The velocity of the particle at O3 with respect to that of the 
particle at Oj is 


where wu’ is the relative velocity of the particle at O3 with respect to 
Oz, and is equal to vo. Thus, 


9 Ca 


c—U 


i 
z Sx yt |. te V1 + Ve = c2 + c(vy + U1) + UV 2 
i (oe 
c 


c? — c(vy + v2) 


U1{V2 


As we can easily see, this is equal to the product 


C+Uzy [C+ V2 
ee eV = 
SUT 


C— va 


So the rapidity @ of the particle at O3 with respect to the observer 


at O, is the sum of the rapidity of the particle at O3 with respect 
to Og and the particle at Op with respect to O;. Thus, 


d6=O+UV. 
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If instead of three particles we have n particles with rapidities with 
respect to each other, ®,, &2, 3,...,®,, the rapidity of the nth 
particle with respect to O, is @ such that 


Calling v the velocity of the nth particle with respect to O,, we may 


write 
c+v 
ef = ; 
C—Uv 
Assuming ®; = ®2 =---=%®, = 9, 
e = er? 
6) 
29 ane _ © + U 
Cc—v 
Calling2" =e", 
Vv zy —)l 
@ gre 


If we call vo the velocity corresponding to the rapidity ®, 


C+ Vo 
peP re 


C — V9 


n 
wa (=~) 
C— Uo 


1 
Assuming that vg = 5° then 


and 


OS ead. Ay al 
Cc gm4+] 3h) 


For n = 5, the particle in the proper system will acquire a speed 


243 — 1 
: 31 _ p99 
Cc 


~ 94341 
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Relativistic Lagrangian and the Lorentz Force 


The momentum of a free particle is (we work in Cartesian coordi- 
nates): 
m da 
p= met = —"—i where ¢= — (0.7.30) 


J1- dt 


The Lagrangian is related to p by 


p=-— =~ where £, is the free particle Lagrangian 


By equating both expressions, we derive the following for the free 
particle Lagrangian: 


5 a 
p= [mort = 5 dé = —moe 1 — &? (0.7.31) 


If we now assume an electromagnetic field interacting on the par- 
ticle, then (the charge q is assumed > 0, the electron charge e = —q): 


= —mc*\/1 — 62 — [o-2.a] 7.32 
£ Moc B?-—@q ; (0.7.32) 


where ®, the scalar potential, and A, the vector potential, are both 
functions of x, y, z, t. 
The Lagrange Equation is: 


dol OL 
dt0z Or 
Let us separate the Lagrangian in the free particle and the interact- 


ing one: 


Note: In Classical Mechanics, the Lagrange is £ = JT — V, where 
T is the kinetic energy and V is the potential energy. For the free 
particle Lagrangian, £, will be equal to T which is > 0. In relativis- 
tic cases L, = —moc?/1 — 8? is negative. To see that these facts 
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are consistent, we expand L for low values of v < c. Let us assume 
6 to be the smallest increment in the expansion: 


ai 
Loy = —mgc? — moc?( 2p BF A nae utr sand 


Ignoring the rest energy, which plays no role in Classical Mechanics, 
ai at aa We 
L, = T, the classical kinetic energy equal to a Mov 


End of Note 
As is well known the Coulomb force is 


F.=qE 


and the Lorentz force due to the magnetic field B is 


where U is the velocity of the particle and qg the charge particle. 
Next, we will derive this last force from the previous relativistic 
Lagrangian. 

The interacting Lagrangian, when an EM field is present, is given 
by 


> U 
£,=-9(@-A.*) 
Cc 
Next, we will see that this expression for the interacting Lagrangian 
gives rise to the Lorentz Force. 
If Q; is the external force, and since 


OL, 
is 
then 
dC. OL. «dl, 


dt Ot Ox dt Oz = 9 
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Now 
Ly = —moc?\/1 — B? 
and 
poe 
Oa er sees” 
since ; : ; 
Vv we tye + 
a —— 
Cc 
and 
a 
Ot 8B? 
Thus 
OL, Mok 
Ox V1 — Bp? 
and 
OL OL; = O® 0/2 JY Zz 
Ox Ox “es meee re =A.) | 
Also, 
ddl; ddO/fq, \ ada qd 
dt 0% ware )= cdi Og | Aer yey? ce ee ae 


So, the Lagrange Equation becomes 


d Mo& 4 OA, OA,. OA,. OA,;. 
ered | Ge "On Oy? a 

O® q/OAz.  OAr. | OAx. 
~ ( Ox ” Ox : Ox :) 


lon" ¢ 


or 


(i) =the Pt af the (Oe) 
dt\./1—B?/  ~—c ‘Ot Ton cl7\ Ox dy/ “\ dr Oz /\" 
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q OA, O® gq 
cot Oa | pA Ale 


From the Maxwell Equations 


Ob OA, 
BS els oy, and B=VxA 
so finally 
d Mox qd, = 
& 2 Ses vA 
eid) Ter ee) eS) 


The last term represents the x-component of the Lorentz Force and 
the first component of the Coulomb Force. Thus, 


Q2 = qEy a “(8 B). 


and d er 
(7p) = 
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Hamiltonian for a Charged Particle 


The Hamiltonian is 


where (assuming Cartesian coordinates and dealing with x only) 


OL Mox qd : qd q 
Ot V1-£B? ¢ faa Per 
px is the free particle (/1 — 62) momentum. Thus, the first term of 
H is 


Since = 
L= —me*y/1 — B2- g(@ — 2 - A). 
Moc B qg( és ) 


The reason for this is that A is the generator of B (the magnetic 
field), and it does not induce any work. Thus, H is the total energy, 
the sum of the kinetic mc? and the potential ©. 
The first term of H is just 
2 


Mov 
FoF + moc?1/1 — 8? 


or 
> mec? — mace +t meu? — mee? 
i, = z— 2 
Vv v 
C2 c2 
moc*(e =) mgcru* 2.4 22 
| — 
= ON OE = (rnb! +?) 
ee = 
Cc Cc 
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MgU D 
where p = mv = ——— . So H may be expressed as a function of 
P /1 — B2 y P 


the kinetic momentum p = mv such that 


H = E = \/(pe)? + mgct + q® 


H = cy/mbc? + p2 + p2 + p2 + g® 


where moc? is the rest-mass energy, the addition to the classic Hamil- 


or 


tonian. 
Now, using the definition of impulse of the coordinate x, 


OL 
Px Sata 
Ox 
we find that it is the sum of the kinetic momentum and potential, 
OL o 
Ox c 
From this we find that 
Cc 


Substituting this into the expression for H, 


H = eye? + (P— LA)? + gb = cynic? + p+ a8 


Cc 


2 
Pp 
D2 
Mpc 


For small velocities v << c, or << 1, the Hamiltonian is 


2 
H— moe? = =— + qo 
mo 


The classical Hamiltonian is H — moc? (the constant moc? does not 
show in its definition). Changing p for P as above, the classical 
Hamiltonian becomes 

1 
2m 


q 2 
t= (p - £4) + q@ 
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Note on Maupertuis and Fermat 


Fermat (1601-1665) stated that an optical path is a minimum. If 
n is the index of refraction and ds the differential light path, then 
We n ds is a minimum between the point Mp and M,. 

Maupertuis (1698-1759) stated a similar principle in mechanics. 
It established that the following integral is a minimum: 

M, 

pigi At 

Mo 


0 


or in other words, the variation of that integral is zero, so 
M, 
é a pigs dt = 0 
Mo 


Again, repeated indexes implies summation. 
For a free particle, this means that 
M, 


mv-ds 
Mo 


m Cages, 
is minimized. (In relativity, m > FP v-ds = aw = 
1, 2,3) in Cartesian coordinates. 
As we have seen, the generalized momentum is defined from the 


Lagrangian as 


~ OL 
(OG; 
The total Lagrange function for a charged particle in an electromag- 
netic field is given by the free-particle Lagrangian and the interactive 


one: 


P, 


LS Lh PLS 
where 
Ly = —moc\/ 1- Be 
and Ls 
La = —q® + gA- . 
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(q is the particle charge). 
If we work with Cartesian coordinates, then the kinetic plus the 
momentum due to the potential A is 


Mot 


Py =... 


The second term of P may be written as a function of the Coulomb 
potential U if we assume that the particle’s energy is conserved. So, 
the Hamiltonian may be expressed by 


H = me? 4+ qu = mc? + 


Therefore, 
(m — mo)c" = q(Up — U) 
but 
™o 
ae 2 
v 
aes 
Thus, 
1 
MoC 5 -1| = q(Up) — U) 
v 
a 


From this we obtain ape, which gives rise to the isotropic part 
Vv 
e 
of the path, while the unisotropy is due to the term ey After 
C 


simple algebra, we obtain from the previous equation, 


(fe +2) 


mo 
Fisekee ic = V)amo( 
3 ee MoC2 
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Canonical Equations of Hamilton—Poisson Brackets and Ja- 
cobi Algebra 


From the Lagrangian we have defined the Hamiltonian H and the 
generalized momentum P: 


P= ae and H = Pig; —L 
Odi 
Assuming P; and q; are variables (called canonical variables by Ja- 
cobi), 


OH | ; and OH. _ oO __d (<) _ _ OP; 
op, % °° oq dq  dt\dq/ ot 
These are the Canonical Equations of Hamilton; that is, 
OH. . OH. 
Gi= OP, and P; = da, 


Given a function of q, P;, and t, F'(q,Pi,t), its variation with 
respect to time is, using the Canonical Equations, 


dF OF OF, OF, OF ,dFOH OF OH 
dé Ot | dq" OP, ' Ot | dqG%OP, OP. dq 
OF 


PF 
ee +{F,H}= a + Poisson Bracket of F’ and H 


Jacobi defined the function S(q, t) as 


Stat) =— f £(a,1) at = 


From this we found that 


Since H. = (Pig; — £), then substituting the two previous equations 
into this expression we obtain the Jacobi Equation: 


OS OS , 


a “(ay % +2) 


144 CONTENTS 


Observe that if we abstract the function S from these expressions, 
assuming H and P to be operators, we observe that H is associated 


O 
with c and P with Bq" These correspondences are the basis of 


qd 
the Schrodinger Equations in Quantum Mechanics. 


Poisson Brackets 
Given a function F'(t,q,Pi) in the phase space, i.e., the space of 


Canonical variables q,P;, its total derivative with respect to t is 


dF OF | OF dq , OF dP; 
dt Ot | Oqdt° OP; dt © 


If the point (q,P;) describes a dynamic imposed by a Hamiltonian 
H, then as we have seen in lecture on Analytic Dynamics, the Hamil- 
tonian Equations are: 


dP, OH |, dq_ OH 
di dg - dt OP,’ 


In this case, : 

dF PF 

—_ = —_ PH 

Go oe 
where are called Poisson Brackets; that is, the Poisson Brackets 
are: 


OF OH OF OH 
The previous equation may be seen as an operator O acting on F, 
that is 


{FH} = 


where d 5 


The Poisson Brackets have the following properties: 


1. Anticommutative: 


{F,G} = —{G, F} 
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2. Distributive: 


(PEC =the FH 


3. Product: 
{FG, J} = F{G,J}4+ {F, J}G 


4. Jacobi Identity: 


ean 1G JT} } a { J, ee Gh} 5 5 {G{J, F}} = 0 
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Mathematicians That Contributed to Analytic Mechanics 


Pierre de Fermat (1601?-1665, French) 


e Pierre-Louis Moreau de Maupertuis (1698-1759, French) 


Leonhard Euler (1707-1783, Swiss) 


Jean le Rond d’Alembert (1717-1783, French) 


Giuseppe Lagrangia ((Joseph Lagrange), 1736-1813, Italian) 


Simeon Denis Poisson (1781-1840, French) 


Carl Gustav Jacob Jacobi (1804-1851, German) 


e William Rowan Hamilton (1805-1865, English) 
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0.8 Lecture 8 


0.8.1 Matrices, Operators, and Elementary Group 
Theory 


Any operator acting on a vector % of the space L will create a map- 
ping in the same space between the vector # and the vector y such 
that 


> 


Of =a. 


If the mapping is one to one, then the operator O may have an 
inverse such that 


or 
00° =0-70= 


where J is the unit operator. 
The operator O is said to be linear if 


O(#+ 7) = O#+07 
and 
O(az) = aOz 


where a is any number. 

If the vectors # and y are given by their components in a frame 
that expands the space L, we may represent the operator O by a 
matrix such that the columns representing # and y are related by 


p\ (Ob OF --) fa 
y? — Q3 eee eee ia 
or 
y' = Ota! : 


The necessary and sufficient condition for the matrix O% to have an 
inverse is that the Jacobian of the previous equation be different 
from zero. 
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~f 3 


If we change the basis from (€7, €3,---) to (€1’, €’,---) by a linear 
transformation A, then the vectors and y change as 


g = Ax 


and 
yf = Ay 


where A is the matrix for the linear transformation. Now, 
OF =y 


SO 
AO# = AG=7 . 


Introducing between O and # the unit matrix A~‘A, 


AOA 'A# = y! 
or 
O'?7' = y’ 
where 
O' = AOA“! 


which is the equivalent operator of O in the new space described by 
the new set of coordinates. 

It is important to see that (a) operating on a vector # of the 
space (which is written as a single column matrix) is equivalent to 
& operating as a single row matrix (which is the transpose of the 
column) on the transposed matrix (@) where (@;; = a;;). 

Example in two-dimensions: 


Q11 42 UT \ ff G11%1 + 1202 
@21 422 v2 A211 + A222 
is equal to 


( a a ) G11 G4q1 \  f ©1411 + L242 
1 2 = 
a12 22 ©1421 1 L2Q99 
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and thus 


Cage) 


In general, (ab) = (b)(@). That is, the transpose of the product of 
two matrices is equal to the product of the transpose of each matrix 
in reverse order. 


The Adjoint Matrix (a*) 


By definition, a? is the transpose matrix with complex conjugate 
elements. That is, 
+ * 


ee 


where * indicates the complex conjugate. 


The following are properties of the adjoint operation: 
Let us (a) and (3) be two matrices or operators; then 


(aB)" = Brat. 
Let us a see this with a simple example, a and ( both being 2 x 2 
matrices: 
“(3 2) ats 2) 
a= ’ B a, 
a2, 492 bai boo 
Then 


a* as ok * 
ees u 21 and Bt = a 21 
Aig G99 12 2 


af = ( 41041 + Gy2b21 4112 + Ay2b22 


@21b11 + A22b21  daibi2 + aeabo2 


So 


and 
* o* ok * ok * * * 
Brat ~ ( 11071 + by1@5 5} 145, + 05,495 
* * * * 2 * * * 
BigQ}1 + DjgQtg 051.051 + d59050 


which is equal to (a3). 
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The Hermitian matrix is one that is equal to its adjoint. Thus, a 
Hermitian matrix (a or 3) or operator is such that at = a and Bt = 
B, so that (aB)t = Ba. The anticommutator and the commutator 
are: 

(a6 + Ba) = (a8) + (aB)* = real matrix operator 


and 
(a8 — Ba) = a6 — (a8)* = pure imaginary matrix operator 


As an example, let us take two Pauli matrices: 


a=(4 A and £6 


Then the commutator is 


w-80)= (10) (5 


while the anticommutator is 


wora0d= (To) (Go 


“(oa 


A unit matrix or operator leaves the scalar product invariant. By 


definition, the matrix (a) is unitary if (a)(at) = J. 

When we deal with complex spaces where vectors themselves may 
be complex functions, like in Quantum Mechanics, the scalar prod- 
uct of two vectors is defined as follows: 


* 


where * indicates the complex conjugate. The product itself is a 
complex number. Like in any vector space, the following properties 


must hold true: 


> 


(7, cy) = c(#,y) and (c#,¥) = e(#,y) 
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and 
(Z1 + Zo, ¥) = (£1, Y) + (Zo, 9) - 
The number c is, in general, a complex one. 
Next, we will show that if A is a matrix operating on 7, then 


(AZ, ¥) = (€, ATY) . 


First, we expand # and y on the space basis €1, €2,--- Then by def- 
inition of the scalar product, the first factor is taken at its complex 
conjugate. Thus, in two dimensions, 
(z,¥) = ie? ye; = aye Gj) = wy! gi; : 
The matrix gj; is (€;, €)). 
If the basis €; were orthogonal, then g;; = 0;; = the Kronecker Symbols. 


For simplicity, let us work in two dimensions and in an orthogonal 


ye ds Q14101 1 A42%2 Y1 
Aq101 + A22%2 Y2 


In the standard scalar product, this is a real vector multiplied by a 


basis. 


column vector: 


Y1 
ok ok ok ok ok ok ok ok 
[( Ay + Ajo%y Ag Ly + AzoX} )| ; ( Yo 


(£, ATG) . 


or 


The generalization of this simple case to higher dimensions is straight- [a 
forward. As we defined previously, the adjoint matrix of A is AT (the 

conjugated transport). Now, if A is a unitary matrix which changes 

basis from €), €,... to é,&,..., we have the following equality: 


(Az, Ay) = (@, AT AY) = (2,9) - 


AtA=T is the unit matrix. 


152 CONTENTS 


Examples of Linear Operators 
1. For any matrix O%: 


O(# + 9) = OF + 0F 
2. If L = fPl(t) dt, then 
b 
Lz(t) = | I(t)#(t) dt 


3. If L =v, then 
Ey Stee 


Some Geometric Uses of Matrices 


Given two vectors, u and v, the area of the parallelepiped formed by 
them is the modulo or magnitude of of the vector product wu x v = 


|vi||@| sin 6. It can be expressed as a matrix: 


€, €) €3 
UXv= Uy U2 UZ 
Ur V2 U3 


where €; is the unit vector along the axis i(1, 2,3). If we have three 
vectors, including w along with v and w, then the volume of the par- 
allelepiped built by the three vectors is the base (|i x U|) multiplied 
by the projection of w on (wu x v): W- (ux v). This is equivalent to 
substituting €1, €2,é€3 by Wy 1, We, w3. Thus, the volume is 


Wi, W2 W3 
V= Uy U2 U3 
Uy V2 U3 


If we square this expression (the determinant of V = the determinant of V(the transpose 


Wi, W2 W3 Wy, Uy Ui ww wu w-v 
y= Uz U2 U3 Wa U2 V2 =) aw a ae 


Uz V2 U3 W3 UZ UZ 
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Substituting w, wu, and v by the unitary basis vectors €), é, and €3, 


gir Gi2 913 
V= 921 922 923 | = 9 
931 932 933 


Thus, the square of the volumen limited by €1, €2, and €3 is equal 
to the determinant of the metric tensor. 
In a change of basis, 


and 
Dae A Ana 
So the determinant of g/,,,, = g', the expression for which is 
g = |A,||Allg =A’g 


and 


Q 
4) 
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Group Theory: A Review 


The operators A, B, and Care said to form a group if 
e AB is another operator of the same group 


e A“! exists such that AA~' = A~1A = J, where A is any ele- 
ment of the group and I is the identity element. 


e Or ABC = (AB)C = A(BC) (The Associative Property). 


The order of the group is the number m of elements in the group 
such that for any element A, A” = J is the identity matrix. Any 
matrix or operator S that maps the space L into L’ and changes 
the elements of the group into another called an isomorphism and 
relates the elements in the following way: 


A= SAS , B= SBS , 


If A, B,C,... can be mapped toa group of matrices M(A), M(B),.. 


in the space L of n dimensions, we say that these groups of matrices 
are the representation of the abstract group A, B,C,... By chang- 
ing the basis of the space L, the form of the matrices will change, 
but there are some parts of them that remain invariant, such as: 


1. The sum of the diagonal elements of each matrix, which is 
called the character of the group 


2. The determinant of the matrices 


The classes of a group are the union of all the elements related to 
each other by a transformation created by an element of the group, 
say U. Thus, Aj, Ag,... form a class if 


A; = UA,;U~! . 


The following are various definitions associated with Group The- 
ory: 


0.8. 


. Conjugate Subgroup: equals aHa~ 
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. Subgroup: If H < G, and HF is a group itself under the same 


operation as G, then H is a subgroup. 


! where a is any element 


of G. 


. Co-Set of a Subgroup: equals aH (left co-set) or Ha (right 


co-set). 


. Invariant Subgroup: also called the Self-Conjugate or Nor- 


1 = H for any a of the 


mal Divisor. H is invariant if aHa™ 
G group. In other words, the left co-set of H equals the right 
co-set, or aH = Ha. That is, H commutes with any element 


of the G group. 


. Factor or Quotient Group: built by the co-sets of n invari- 


ant group. 


. Simple Group: a group which has no invariant subgroup. 


. Semi-Simple Group: a group that has invariant subgroups, 


none of which are Abelian. 


When the representation of the group G is of the following type, 


a1, a2 | ag 
P=) by 03: :bs 
0 0 C3 


then when applying D to a vector (x,y,z), we get 


and 


x x 
(D)| y J=I| y’ 
0 0 
0 a! 
(D)|} 0 |=] y 


R 
R 


156 CONTENTS 


We note that the component (, y) of the vector transforms along 
the same components. In other words, D(G) leaves the subspace 
(x, y) invariant. However, the space perpendicular to the plane (z, y) 
has a z-component that, under the transformation, is mixed with 
the rest of the space such that 


(0,0, 2) > (2’,y’, 2’) 


while 
(0) yO) 


If, however, D(G) is such that 


D(Gi) = 


then both spaces are invariant. (x,y) > (x’,y’'), z > 2’, and D(G;) 
is thus a reducible representation of G, which implies invariance in 
some spaces. 

If for any vector (x,y, z,u,v) the matrices D are such that the 
vectors transforms into a vector (2’, y’, 2’, u’,v’) then D is a reducible 
representation and can be written as 


a4, a2 413] O 
Az, 22 A293 | O 0 
D= 1 az31 a32 433 | O 


0 0 0 boy boo 


Thus, 
D = D(a) + D(b) 


and the space L is also separated as a sum of spaces, 
L=L(a)+ L(8). 


When D becomes a unitary matrix, we have reached full reducibility. 
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We will now show that a unitary matrix cannot be anything but 
completely reducible. Assume, however, that this is not so, and 
let us work in three dimensions. Thus, if D is assumed to not be 
completely reducible, then it looks like 


Qi, 12 | 413 
D= Q21 422 | 493 
0 0 033 


But if D is unitary, then its rows and columns are orthogonal to 
each other. So, the product of row one and row three must be zero. 
This implies that 

13 + 433 = 0 


and assuming that a33 4 0, this therefore means that a,3 = 0. The 
same principle holds for row two and row three, making a23 = 0. 
Thus, D must look like 


a1 a2} O 
a21 423 | O 
0 0 | 033 


which is reducible. 
For finite groups, there is always a unitary representation of the 
group (the proof for this is omitted). 
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An Abelian Group is one where all the elements commute with 
each other. Thus, if D’ are the components of the group (or their 
matrix representations) then 


DP Sy De 
If the eigenvalue of D’ is \; and its eigenvector Z, then 
Dae: 
Since D* commutes with any other D of the representation, 
DD 2] Dita NDE 
but since DD'# = D' Dz, 
d,(DzZ) = D*(D#) 


which makes D’, any element of the Abelian Group, just a number 
A;. This means that the dimension of an Abelian Matrix is one (just 
a number). 

In general, if the representation of the group is an n x n matrix, 
we say that the dimension of the group is n. As we stated before, 
the order of a group is m if for any element of the group a 


a:a-:a:...=a”"=T1 


where m is supposed to be the least integer. Example: For the 
rotational group a = e?, m= 4. 
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0.8.2 Review of Lie Algebra 


Lie Algebra was developed by Sophus Lie (1842-1899), a native of 
Oslo, Norway. 


Infinitesimal Transformations 


x 
a am da 
x aa xv’ +dzx 


Let us take the following example: one parameter, a, for the 
transformation that takes x > 2’ = f(x,a). The transformation is 
uniquely determined by the parameter a. 

We can go from x > x’ + dz’ by 1) f(a,a+da) or 2) by another 
transformation that starts at x’ and whose parameter is ), i.e., x’ zu 
x’ + da’: 


xv >a'+dz' by f(2’,b) 


or 


x’ +da' = f(a’,b) . 


Since this is equal to two successive transformations (zx to 2’ and a’ 
to «’+dz’), the new parameter b must be a continuous function of a 
since the transformations are assumed to follow a continuous path. 
Let is assume that b is a small variation, say some da. Then, b © da 
may be expressed as a function of a, 


bx da=x(a)da. (0.8.1) 
By using a Taylor Expansion of f(2’,a) with variation b = da, 


x’ +da' = f(2',da) » f(2’,0)+ (Car) _ ot =2'+u(2')da. 
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Of (2, a) 


We have used the fact that f(2’,0) = 2’ and set ( a ) equal 
a=0 
to u(z’), i.e., 
(’) = (cea) 
7 Oa a=0 

and 

da =z oa. (0.8.2) 
From (0.8.1) and (0.8.2), 

da’ 

rele x(a)da . (0.8.3) 


Thus we have separated and related the differential variation of 
the ”space” (2) which is being transformed and the element of the 
group that is specified by the parameter a. Next, we will see how a 
function F(x) changes as x > x + dz: 


OF 
Pee 
d a dx: 
but from (0.8.2) 
dx = u(x)da 


SO ‘ 
dF = (u(0)a5-)F 


6) 
In this way, we introduce the operator X = u(x)— so that 


Ox 
dF =6axXF . 


O 
X = u(x)— is called the differential operator of the transformation 


a 
x’ = f(ax,a) that allows one to determine the differential variation 


of any function F(x) as the original transformation f changes due 
to a variation of its parameter a. The transformation f forms a 


group, each element of the group characterized by a parameter a. 
The corresponding differential operator X also forms a group where 
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the operation among them is the commutator. They are called gen- 
erators of the transformation. 
If F = 2’, then 
xv =(1+X6a)X . 


Another example is the two-parameter transformation: 
xv =ar+b. 


The identity is a = 1,b = 0. The infinitesimal transformation is 
represented by the parameters da and 6b, so dx = 2’ — x is 


dx = roa + 6b 


As we have seen, the function u(x) in the expression of x is the 
coefficient of the increment of the parameter, so 


tite) =e 
and 
tts 6) ls 
The corresponding infinitesimal operators are 7 = uae 
es 
O 
Dye ee 
Ax 
and F 
Ky Se 
° Oa 


The variation of a function F' due to the variation of both parameters 
is 


dF — (OX + OpXp)E . 


The commutators of X, and Xz are 


i, Me) = he, = 
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and 5 
X,) = -— = —X,. 
[Aas o| Ox b 


Thus, the commutator brings back another infinitesimal operator. 
That is why the infinitesimal operators form a group under the op- 
eration/commutation. 

Another example is the one parameter rotational group of trans- 
formation in two-dimensions: 


x’ = xcos® — ysin ® 


y =xsin®+ycos®. 
For the limit @ — 0, expansing the functions of ® gives: 


Ocos ® Osin ® 
‘xs (1 DP — ®))=x—ydo®. 
ve xa(lt ( 5 )_8 y(0 + ( a0 8 ))=ax—-—y6d 


y ~xcdO+y. 
So, 
dx=2'-x=—yd® or u(x) =-y 
dg=7 y= 208° or uly) 2 
and 5 3 
or u(x) a =U ae 
0 0 
Ries iy) ee 
The sum of both operators is 
O 0 
ean one RE 
“Oy aye 


which is the angular momentum operator around the z-axis. Also, 


0 
dF(az,y) = 5g + ay =|- yore ey F=(6®)XF. 


0.8. LECTURE 8 163 


In the last equation, we have substituted dx and dy from the above 
expression. 

If F = x? +’, the norm of a vector (x,y), the operation of X 
acting on F' leaves F’ unchanged (the proof for this is straightfor- 
ward). 

As another example, let us analyze the Lorentz transformation 
in two dimensions, t and x. This transformation may be expressed 
by the following matrix expression: 


La _{ cosh(g) —sinh(é) Xo 
i —sinh(€) cosh(€) Ly 
where Xo = ct, x; = x, and the hyperbolic functions of the parameter 


€ for the case of special relativity is given by 


cosh(€) = y = me 


and ; 
sinh({) = ot 


It is easy to see that the transformation leaves xj — x7 invariant. 
For three dimensional space we have: 


The commutators of these operators are: 
[Xo3, X13] = X12 ) [Xi3, X12] = X93 ’ [X12, X23] = X43 . 


These operators ”X” are sometimes indicated by —7T and the trans- 
formation Aj2, as an example, may be written as 


cosh —sinb 0 
Ay2(8) =] sinB cosB 0 | =exp(—iBTi2) . 
0 0 1 
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For a small 6 (63): 


or 
eo 0 2 O 
& \- 1T1260 =1 —i 0 O 6B 
v 0 00 
or 
0 -—2z O 
Ti. = 2 0 O 
0 0 O 


In the language of group theory, the previous transformation is 
called an SU2: the ”S” stands for special (the determinant is equal 
to 1) and ”U” stands for unitary, that is, it leaves the scalar prod- 
uct invariant. The unitary matrices have inverses equal to to adjoint 
(aj; = a%;) and its rows and columns are orthogonal to each other. 
The number 2 stands for the fact that they are 2 x 2 matrices. 
The SU2 groups are homomorphic to the three-dimensional real 
rotation which forms what is called the O3 group. To see it, let us 
start with a general complex 2 x 2 matrix u which belongs to the 


SU2 representation: 


a ab ag fp a ab 
we ( 4 = and wu a5 | 


Since a and b are complex, there are four parameters, but only three 
are independent because of the fact that the determinant should be 
equal to 1: 

aa* + bb° =1. 


Now let h be one matrix of the SU2 representation; let’s expand it 
as a sum of the Pauli Matrices with coefficients x, y, z: 


ee ce e294 acipoee 1 
NTT! Gi ake Sh SES es a aya de OR <RE AN ail 
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So, 


hay hy z xr ty 
h = = x z= : 
( ho, hye ia ae he Try —z 


Since the determinant of h is equal to 1, we see that 27+y?+2? = 1, 
which expresses the conservation of the scalar product 7-7, where 
r= (a,y,z). The values of x,y, z as a function of h;; are 


= hie + hoy hay a Aig 

B= — oH SOC Hh = hee 
2 2% 

Since x,y,z are real, the matrix is Hermitian. Thus, h is given as a 

function of x,y,z. If we change coordinates (x,y,z — x’, y’, 2’) we 


generate another h’ matrix such that 7’. ¢ = (h’); then, 


ipl ces sah y= "athe Z =n, =—h,. 
Now we introduce the matrix u to go from h to h’, h’ = uhu"!. 
Then (2’, y’, 2’) will be expressed as a function of (x,y,z), and the 
parameters a and b of u. The result is a rotational matrix that 
changes (x,y,z) into (2’,y’,z’). Thus, from a 2 x 2 SU2 matrix 
we generate a 3 x 3 matrix of the O3 group representation that, of 
course, preserves the scalar product 7-7. Let us assume a simple 
unit matrix u, diagonal only and a function of one parameter, a. 


The most general one will be 


166 


CONTENTS 
From this, 
gz! cosa —sina 0 x x 
y’ sina cosa 0 (Ta ese 8 
ae 0 0 1 ys z 


So the rotation R is in the (x, y) plane around the z-axis. 
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Lorentz Generators 


By expanding the hyperbolic function of € and reducing it to a small 
variation fo €, d€, the Lorentz transformation becomes: 


(%)=((5)+(2 a) a) (2) 
vee ec. 


The infinitesimal operators are 


or 


O O 
0 ss "Axo “Ot 
and 5 
Xi = anya = a : 


In general, we can write the differential transformation due to an 
increment of the parameter € as 


The matrix (lo1) is called a generator of the transformation group. 


QO) sh 
m=(0 0] 


In a case of many variables x and one parameter €, 


In our case, 


dxo Xo 
=(lij) | 2 | 0g 
dz, T 


and the matrices / are independent of x and €, and again they are 
called generators of the transformation group. 
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It is easy to see that the function u;(x) and the (J;;) generator 
are related by 
U(x) xo 
= (lis) 


where (/;;) is a generator of the group and the infinitesimal operators 
are 


a Via 
mi i 
The commutators of X; are equal to the commutators of the gener- 
ators (1;;). 
In the previous case, the generator is 


ln = ( ae 
1 O 
but if we work with the rotation in four-dimensional space, there are 
six generators in the Lorentz group corresponding to the six inde- 
pendent rotations which leave the six independent planes (291,71), 
(X01, 2), (01, £3), (41, £2), (2,3) invariant under these rotations, 
and the main invariant may be written as (x? + y? + 2? — t?). 

The generators are such that their sum and product as well as 
their commutator belong to the same set. These properties charac- 
terize what is called Lie Algebra. The number of commuting genera- 
tors is called the rank of the algebra. The Lorentz group has a rank 
equal to 2. In general, its commutator obeys the following equation: 


Ue bvee| = Suv ne oe Gratin = Gualyy 3 Orv wes 


where the g’s are the metric tensor which for special relativity in 
rectilinear, orthogonal coordinates are are g,, = 1,—1,—1,—1 for 
the diagonal elements and zero for non-diagonal elements. 

For example, take uy = vy = 0, A= 3, and a = 2. This means the 
j and v correspond to the time variable while \ corresponds to z 
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and a to y. Thus, 
Le big = [lo3, loa] = g32l00 + Gool32 — g32loz — g3olzo - 
The only term not equal to zero is gool32 = 132 = ley. So, 
[lo3, loz] = Use « 


In the notation of infinitesimal operators, the commutator for this 


transformation is 


[X3, Xe] = Ai 
where 
a oO Oe oo a 
5 ater age: 
Dae gpg Se = Var ag, RS ea 
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The transformation f(z — x',a) may have more than one pa- 
rameter a and the function u(x) as in the previous case may refer 


to many variables. So let x represent the set: 21, 2%2,...,%p and a 
represent the set a1, d@2,...,a,. So equation (0.8.3) may be written 
as 

dxo = wiz(2’)Xp(a)da; (0.8.4) 
where? =D iniy 8, 22n, 7 ond 1 and ares 1s 202s Biwes lice ey Pe 


’ 


The infinitesimal operators are expressed in their general case as 


0 


Xs = PA . 
i 


It may be shown that, as in the previous example, their commutators 
are linear combinations of all infinitesimal operators. So if repeated 
indices imply summation, then 


[Xp, Xo] = CE Xs 
where the coefficients Ce are real constants and obviously: 
Ch = Cr. (0.8.5) 


Using the Jacobi Identity for the commutators of X,, X,, and X;, 
we find that 
ChCn + C8 C7, + C4hCi, =0. (0.8.6) 


po pr oT pp Tp~ wo 


Lie showed that by knowing the C’s that fulfill equations (0.8.5) and 
(0.8.6), we may eventually solve for (0.8.4), which will eventually 
unfold the transformation of f(z — 2’). The C coefficients are 
called structure constants for Lie Algebra. 

Let us look at the previous analysis regarding the Lie Algebra 
coefficients of the Lorentz group in a more general manner. Let us 
see how a metric tensor of the group comes about in relation to the 
Lie Algebra coefficients. 

The fundamental operation for the group of infinitesimal opera- 
tors is the commutator. Let A and B be two elements of the group. 
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One element of the group S may generate the other members of the 
group by using the commutators of A and S and B and S. Thus, 
let us call P4S' the commutator: 


(Pa)S = [A, 5] . 
If we use a set of Xo,...,X,,...,X, to expand A and S: 
A= OGih, Bnd. Oo = syXy 


then 
(Pa)S = ay8,[X,, Xi] = OSC Xa : 


With a similar expansion for B = b,X,, we calculate [B,[A, S]], 
which may be called (Pg P,4)S. Thus, 
(PaaS = bg Gy SvC i [Xo, Xa = GisiC yO ean: 


pv~oa 


Let us sum over ju, 0, and a, so that me may generate a matrix with 
indexes of of v and ¥: 


(PsP Oi =a,b,C2C 


pu~ ac * 


The s,, coefficient is part of the operator S. 
The trace of (PgPa)vy is 0,,(PaPa)v7dv7, 80 


Trace( Ps Pa) = Qubc9uc 


where 


GeS C2 CY 


pv ~ ao 


and the sum over a and yp is understood. 

This trace resembles the scalar product of two vectors, whose 
components in the system are a, and b,, which are precisely the 
components of the infinitesimal operators A and B in the base 
SC pepere. ereaet. See 

As we have seen before, the C’s are real and antisymmetric in 
the interchange of the lower indexes. But from the expression of 
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Juv, the interchange of js and o does not vary the value of g,,,. So, 
this metric matrix is real, symmetric, and leaves invariant the scalar 


product: 
(Ay B) = GuitnGp: : 


If a matrix is real and symmetric, it could be diagonalized by a real 
transformations. 
If the determinant of g,, is zero, then the scalar product 


(A,B) =0. 
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Next, we state some results on the Lie Algebra of a semi-simple 
group without proof: 


1. In order for the Lie Algebra to correspond to a Semi-Simple 
Group (that is, a group with invariant subgroups, none of which 
are Abelian), no element of the infinitesimal group should be 
orthogonal to all other elements of the group, or, in other words, 
the determinant of the g;; must not be equal to zero. This is a 
theorem proposed by Elie Cartan (1869-1951). Here it is worth 
observing that an Abelian Lie Algebra is of no major interest 
because the commutator of its infinitesimal operators vanishes. 


2. For a compact Semi-Simple Lie Algebra, g,, must be nega- 


tive by definition. Thus, we may choose g,,, = —d, , where 


LV) 
6 = 1 for zp =v, and otherwise 6 = 0. In this case, the struc- 
ture constants Ci, are antisymmetric under any interchange 
of indexes, and of course (A, A) < 0 for every elements of the 


algebra. Also, 


Trace(PgPa) = Qube6uo = —Gyby 
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Lorentz Generators and Lie Algebra 


Let us review some developments of our previous lecture on rela- 
tivity. The Lorentz transformation that changes the coordinates 
(t,x, y,z) to (7,€,7,¢) is a rotation in the plane (xo, 21) where the 
new coordinates are ro = ct, ©] = ©, Xo = y, and x3 = z. Thus, in 
matrix form: 


aN CT he es OD ct 
<a (eed Wa fe = y 0 0 r 
/ 

z t] 0 Or ahd y 
v3 G Oc ie eA z 


From this expression, we obtain again: 
v 

cT =cly — -yx 
Cc 


or 


and 


= —7- (et) ey Sale = Ob) 


The inverse rotation is done with the inverse matrix. That is, 


Xo ct 7 1 0 0 CT 
Co rn ae y 0 0 E 
mz) |y] | o 0 10 y 
v3 oO: 105204 ¢ 


Whatever the coordinates, the fundamental metric tensor is the 
same: 


Gij = 


Oo oo. = 
j=) 
| 
— 
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This metric guarantees that the velocity of light in systems moving 
one with respect to another with constant velocity is always the 
same. That is, 


dx* =c’dt? , dé? =c*dr’ etc. 


for any relative velocity v. The transformation from (t, x, y, z) coor- 
dinates to (7, €,7, ¢) will, however, depend on v and it may be seen as 
a rotation in the planes (29,21), (%0, 2), (0,3), (%1,%2), (11, £3), 
(x2, #3), or any rotation in the four-dimensional space. This last one 
will be expressed as a product of the six fundamental rotations by 
some constant angles. 

These rotations form a group, the continuous Lorentz Group with 
six dimensions. They also form a Lie Group, so any differential rota- 
tion is expressed in the six fundamental rotations, called generators, 
of the Lorentz Group. For instance, the previous rotation, for small 
angles €, may be expressed by the first two terms of the Taylor 
expansion 


aoi(&) & 5:3 + (lo1)ij€o1 


where do; is a component of the rotational matrix, and 


010 0 
a 1000 

in = [ggan(8)] = 0000 
000 0 


Knowing /o;, we may generate the entire transformation in the plane 
(aoa): 

The previous transformation from (29, 7) to (xj, x) may be ex- 
pressed by using hyperbolic trigonometry of an angle (3: 


CT cosh3 sinhG 0 0 Gt 
€ | | sinh@ cosh8 0 0 

n |. 0 0 10 y 
¢ 0 0 0 1 Z 
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The inverse transformation corresponds to an angle —3. We define 


) 1 
uh ==: So: oS. 7 
p Cc ? 1 sinh? 8 
cosh? 3 
Since cosh? 3 — sinh? 8 = 1, 
sinh? 8 _ 1 = v? 
cosh? B cosh? 8? 
Thus, 
1 _ vu nl 
cosh? 8 eC 4 
or 


2 
v 
cosh? 8= 7? and sinh? 6=77-1= ia . 


The rotations in the planes (2, 2), (#9, 23), (%1, 2), V1, 43), and 
(x2, 23) are expressed by the following matrices: 


cosh3 0 sinh@ O Oo Ok 0 
a= 0 i! 0 0 Ajet 00 0 0 | 
sinh’ 0 cosh@ 0 1 00 0 
0 0 0 1 00 0 0 
cosh3 0 0 sinh 0 nO eG ew 
| 0 1 0 0 Jou-[B584) 
0 0 1 0 Os se 208 20) 
sinh’ 0 0 cosh hse OF OG 
1 0 0 0 Oo Qk. reg 
0 cosh@ sinhG O 0 0 1 0 
ai2 = . > lhea= 
0 —sinhG cosh O 0 -1 0 0 
0 0 0 1 | | 0 0 0 0 
1 0 0 0 0 0 00 
aS 0 coshG O sinh§ | ees | 0 0 01 
0 0 0 0 0 0 0 0 
0 —sinhG O cosh6 0 -1 0 0 
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1 0 0 0 00 0 0 
ign 0 0 0 0 Sik 00 0 0 
0 0 cosh@  sinhf 00 0 1 
0 0 —sinh( cosh 00 -1 0 


It can be shown that the commutator of any two of these generators, 
say l,, and I,,, is given as a function of the other generators and 
the metric tensor. That is, 


Le boa = Guelisp = Guelvo a Qu plus ae Geli 


where 
10 0 0 
Got Wt 30 
2 |G iG: - St! 
Oe Ge: JOP. at 


The set of the generators is called the Lie Algebra of the Lorentz 
Group. 
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0.9 Lecture 9: Fundamentals of Today’s 
Cosmology 


0.9.1 Differential Operators in Curvilinear Co- 
ordinates 


Call ut, u?, u® the curvilinear coordinates in three dimensions. €}, é5, €3 
are the “unit vectors” on the three axes, which do not necessarily 

correspond to unit values. Thus, a vector F which begins at the 

origin P is expressed in curvilinear coordinates as 


> 


F=wé+uv76+u°6 


or 
dF = duté; + du7é + dué . 


In general, if F changes by dF , then its extreme Q will describe 
a curve whose length ds is such that 


ds? = gijdu'du! 


0.9. LECTURE 9: FUNDAMENTALS OF TODAY’S COSMOLOGY 1 Oo ———! 


where repeated indexes imply summation. g is the metric tensor, 
and gi; = e; - é;. From now on, we will be working with orthogonal 
coordinates such that g;; = 0, as long as i ¥ 7. 

An example in spherical coordinates: 


esl 


<y 


ul =r, u? =96, and u? = W. With unit vectors of #, 6, and w, 


dF = drf + (r0d0)0 + (rsin Odw)y 


and 
ej = 6 =(r0)6 , &=(rsind)). 


> 


So, 91 = 1, g22 = r7, and g33 = (rsin@)?. Since f, 6, and 2 are 
all orthonormal, all the other components of g are zero, for i # 7. 


Thus, the g matrix is diagonal, and the determinant is 
lg] = 911922933 - 
The element of volume within e€7, €5, and e3 is dv, 


du = /g du'du"du? . 
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The length of an arch along the coordinates u; is 
(dl); = hidu’ 
where h; = \/gii- Thus, in these spherical coordinates 
dr = ejdr + e9d0 + e3dy 


and 
ds? = (€j - €j)(dr)? + (3 - €3)(d0)? + (€3 - €3)(dwy? 
or 
Gi = Gi EG 


The differential lengths along each axis are 


dl, = |éj|du’ 
diy = |é3|du? 


diz = |é3|du? 


Thus, h; = |é|. hjdu’ is a length, and h; = \/g;; makes each of the 
terms of dr’ have the same dimensions (length). 
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The Gradient V of the Function ® 


_ 0b ,, 8%», 08 4 
~ ou h2Ou? ' h30u 


V® 


where hy = 1, ul =r, hg =r, u? = 0, hg = rsind, and ue = w. So 
in spherical coordinates, 


O® . 08 a | O® 7 
"580 ~ rsinddw” | 


The Rotational Vx of the Vector F 


p 6 ob 
3 O O O 
Veo h,Ou! h2Ou2 h30u3 


Be Fy F; 


where F\, Fy, and F3 are the components of F on the three orthog- 
onal directions with unit vectors f, 6, and w. We could resolve this 
matrix for V x F' and find that 


vxF = 4 OF; OF, -4( OF, OF; ) -a( OF, OF, , 
= h2Ou? h3du? h30u® h,Ou! h,Ou! h2du? 


where, as before in spherical coordinates, hy = 1, hg = r, and 


h3 =rsin@. 


The Divergence V- of the vector F 


From the Gauss Theorem, we may derive the divergence of a vector 
as the limit of the flux of that vector on an infinitesimal surface, 
containing the point of interest, per unit volume. 
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side2 


Let us calculate the flux through two opposite sides formed by 
(side 1) [rd6, (rsin #)dy| and (side 2) [(r + dr)d@,r + dr].The unit 
vector perpendicular to side 1 directed out of the volume is (a) x 6) 
and for side 2 it is (9 x w) = —(~) x 6). The total flux through the 


two sides is 


[(r + dr)? sin 6 — r? sin 6]dyd 


F-(6x #) 
6 x q) sin O[(r + dr)? — r?|dwdo 


-(6 
ey ra 
Since (4 x ¢)) =, and F'-# = F,, we find that the infinitesimal flux 
through two faces per unit volume is 
— —[F,(r? sin 0)|drdédy 
F: 


where dv = (r?sin@)drdédw. Thus, the contribution to the diver- 
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gence from the two sides is 


~ 10 ‘ 
(Ver. = 725, Err ) 3 
Similar calculations for the other sides reveal 
1 1 1 
O F O O 


— —(F,r?) + ——(Fosin 0) + —__—F, . 
2 Op vane ap er ate Oo u 


V-F= 
The Laplacian A of the Function ® 


Let us take a function ® = ®(u', u?,u?). The Laplacian of ®, A®, 
is the divergence of the gradient of ©, or V -(V®). Since 


Om pla?» 1 O® 
dr ' rd6' "rsind dw 


Pies, POPOV. OL OS > ae OL 08 
Ve or (Ge ) aC oo sin#) “sin 0 Oy a) 
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0.9.2 Potential Theory and the Hamilton-Jacobi 
Equation 


Assume a differential first order linear equation of the following type: 
we look for the unknown y(t, q) if it comes that 


F(t,q,p,y) = 0 (0.9.1) 
where g represents a set of independent variables qi, q2,... and p 
represents a set of independent variables p,,po,.... p is related to 
y(q,t) by 
_ oy 
p= aq . 


p and q are called canonical variables. 
It is possible to solve (0.9.1) by solving two equations called the 
Hamilton-Jacobi equations. In fact, 


dF OF OF. OF. dF dy. 


de oR oye On Oh on. 


Oy 
Ot - 


O 
By definition, 5, si 
qd 
In simple cases, F' and y do not depend explicitly on t, such that 


Bq! Op? a 
a O OF OF 
F 
T — 2 
(5, me ap? (0.9.2) 
If we make F ; Ae 
F F 
—p= d —e= 0.9.3 
rane pom (0.9.3) 
then we have 
pq—qp=0. 


Thus solving the canonical equations (0.9.3) we solve (0.9.2). 
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If y is not part of F, ie., if F' is independent of y, then (0.9.3) 
becomes OF OF 
—=-p d= 0.9.4 
pe a (0.9.4) 
Assume F' is the Hamiltonian of a dynamical problem. Let us 
take the example of a simple harmonic oscillator, where H is the 
total energy of the system: 
Sa 


Im 2 


Directly from this dynamical equation, assuming that the Hamil- 


tonean is constant, we have 
— =0= p— —- Aqq. (0.9.5) 
m 


Now, ee q and Aq = p, so (0.9.5) becomes simply 
m 
pg—pq=0. 


Besides, 
OH OH. OH . 


Oe ba ee 
Comparing with (0.9.5), 
OM. <P i, OH : 


which agrees with (0.9.4). 


Relationship of Potential Theory to the Poisson Bracket 


If we have any function F(p,q,t) = 0, then 


OF (OF, OF\_ 


Ot 
and using the Hamilton equations 
iF F OH FOH 
a= -(3 OH  OFO =the}. 
Ot Oq Op Op Og 
where {F, H} is called the Poisson Bracket of F and H. 
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Differential Equations for Potential Theory 


Newton’s Law of Gravitation says that the attraction potential be- 
tween two masses M, and Mz is directly proportional to the product 
of the masses in indirectly to the distance between them: 


MM 
a 


ee 


Forces are related to their potential by the fact that F’ is propor- 
tional to the gradient of U, VU. In three-dimensional Cartesian 


space, ; f P 
UL, w_ wv 
VU = Ae x ay Yo Ay 


where Z, y, and Z are the unit vectors along x, y, and z. 
The Laplacian of U is AU such that 


a 


_ UU, PU 


A 
- Og? Oy 


AU may be expressed as the scalar product of the gradient by itself: 
La VV, 


1 
The Laplacian of — is zero, 
r 


Ke 0 
rT 
so the Laplacian of the Newton Gravitational Potential is zero. 
The gravitational force F' between a uniform-density sphere M 
and a point mass m which are some distance r away from each other 


in space is 
Mm 
re 


F=VU=-G 


But if m is inside MW, some distance r away from the center of M, 
then the gravitational force between them is 


Far. 
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Some Typical Potentials in Classic Physics and the Corre- 
sponding Differential Equations 


The work of any force F' when applied from point 1 to point 2 is W 
such that 


2 2 
w= F(F)-dr = [ 9G ap = (pS) 


Conservative forces are those in which the work W done by them 
does not depend on the path taken. Obviously if VU is integrable, 
i.c., has a primitive, then this condition is held and the force is 
conservative. 
In addition, a conservative force implies that U is differentiable, 

that is 

OU U 

dxdy Oya * 
In fact, by Stoke’s Theorem applied to VU, 


gvu- a= f ds [V x (VU)] . 


The first integral is take on a closed curve which starts in point 1 
and ends at the same point. Since this integral is independent of 
the path taken and it is always zero, then V x VU = 0. 

If V x VU = 0, then it means that the following is also zero: 


0 Od O 
dx dy Oz| (8U BU > 
eu au aul” Ce Svan) i 
Ox Oy Oz 
Thus, ale must equal os etc. Therefore, a conservative force 
OxOy OyOu 


means a differentiable force or V x VU = 0. 

If this is not the case and V x VU 4 0, the vector f= Vx VU 
is generated by the circulation of VU around a closed curve. 

The Laplacian of some potential functions may or may not be 
zero, independent of whether they generate conservative or non- 
conservative forces. 
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Examples 
1. 
VU =27 7+ y+ 0k 
Then, 
i gj ik 
VxVU=| Ox Oy Oz | =2k(x-y) 
y~ 2x? 0 


which is different from zero unless x = y. 


Its Laplacian is zero. In fact, 
V-VU = 0,U + OU + O7U =0. 


In summary, this example corresponds to a non-conservative 


force, V x VU £0 with a null Laplacian, V?U = 0. 


2. If instead of specifying the force (i.e., VU) we specify U as a 
simple differentiable potential, we know that if 


Vx VU =0 


the force VU will be conservative, but its Laplacian may be 
different from zero. For example, take the potential 


U=2+y?. 
OU OU 
For this, V x VU = 0. However, since —~ = 6x and —~ = 
Ox? Oy? 
6y, then 
AU = 6(x+ y) 
which is only zero for x = —y or x =y=0. 


Let us take another example, 


Ua(e?+y?4+2)?. 
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Then, 

_ 9 21 21 42\-3 0 2124 42\-3 0 
AU = F(a)? +2?) Hts (way t2) Alt 
and also 

VxVvU=0. 


Thus, the potential is conservative and its Laplacian AU is 
Zero. 


3. Let us take a Hook-type of potential, where F' « r, or 


pe) 


wy 2 2 2 
Thus, 
VU =ai+yjt+ zk 
and 
VxVU=0. 
However, 
a a @\/x2  y? 2 
Au = ( | | )( t — 4 )=3 
On? OOF CFP 2 2 2 


Therefore, the Hook force does not hold a potential whose 
Laplacian equals 0. 


1 
If the force is Newtonian, where F’ x —, plus a Hook-type 
r 


force, say . 
F = ~3 + br 
- 
then the potential will be such that 
FP=—-VU. 
And so, 
<6 bie 
or 2 


and the Laplacian of U will be 
au = 55 (F2) =-a 
Or 
AU = —3b is a Poisson-type equation. 
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0.9.3. Simplifications Regarding the the Space- 
time ”Length” Invariant, the Geodesic Equa-Jggaaaagy 
tion, the Red Shift, Black Holes, and the 
Einstein Approximation 

Next, we will be looking at different approximations that will be 


indicated by $1, $2, ..., etc. and their application to: 


1.) The Space-time ”Length” Invariant 
ds? = G(a" >: (0.9.7) 
S1: Assuming that the prevailing term is go9(x°)?, then 


de =e (dr)? = gists ys. (0.9.8) 


io\2 
S2: Let us use only two coordinates; x? and x? are neglected. Thus, 
ds? = goo(dx°)? + gu(dz')? 
or 
Ce = goo(v")? + gu(v')? : (0.9.10) 
2.) The Geodesic Equation 
— = Thue? . (0.9.11) 


As before, let us use only x° and x! (neglecting x? and x?); thus pu 
is 0 or 1. So, 
dv® 
dr 


0 


= —27),v°u' = —g goov°v" . (0.9.12) 


Now, 
OGoo Ogoo OT 4 Ogoo 1 


Oka Saal OT U1 


0.9. LECTURE 9: FUNDAMENTALS OF TODAY’S COSMOLOGY =_=_™_! 


Substituting into the previous equation, 


dv? 1 0900 6 
=— De 
dr Joo OT 


Integrating this equation, we get 


Inv? =—In G06 
or 
v° goo = €c (0.9.13) 
‘ ' ; ‘i dt 
where € is a constant and c is the speed of light and v" = cr 
T 
dt 
Observe that for € = goo, oo = 1. In a flat space-time, where 
f 
dt 1 
dr a v2 


ae 


it means that v < c. So if our body is accelerated with an acceler- 


dv? tak ; 
ation _ then € may correspond to its initial velocity. 
T 


For ps = 1, the geodesic is written as (neglecting coordinates x? 
and x?) 


dv! 1 1 
irs c= —Too(v°)? a Pies ca (=9'"goo,1)(v°)? = (<g" 11,1) (v')? 
dr 2 2 
which can be simplified to 
dv! 1 
ae => 59 (Goo.1(v)” = gua(v')?) ‘ (0.9.14) 


To derive the Newton Law, we neglect the second term, assuming it 
to be much smaller than the first one. Thus, 


9° Joo(v°)* . (0.9.15) 
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But 
dot. .dohdg’’ dot, 


dr dxr®dr ano 


so substituting into the previous equation, 


da “Thy 6 
ano 94 Joo,1v - 
Now using equation (0.9.9), 
GOs 2. iad c 
dro I Daa : 


Since g!” = 0 for v £1, then using the lower indexes technique 


9°" Jo0,1 = 9" Goo, = —900,1 - 
The minus sign is to ensure that g™ brings a negative sign: 


dut € Jo0,1 
do. 8 /G00 = —c(/G00),1 - 


Since dx° = cdt, we find that the acceleration along x! is 


du! ; 
“de ='<-¢€ (./900) 1 : (0.9.16) 

Equation (0.9.16) resembles the Newton equation of motion if 
,/goo is proportional to the gravitational potential function of z?. 
Since goo for small curvatures will be close to 1, we can write 


goo = (1 + 20) 


where ® is a small quantity. Next, using the Taylor expansion we 
can express the acceleration for a unit mass as 


a= —(VGo0) a2 = —2(V1 + 28) 1 & —2(1+ 6), = -—2V16 
(0.9.17) 
using the approximation: /1+2®@+1+4+ 9. 
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3.) The Red Shift Due to Gravitational Effects and the 
Expansion of the Universe 


Using the simplification of (0.9.7), we can write 
(As? Ser? = gone T” 


where 7 is the period of the light wave at the point where light 
is emitted and T' is the period observed after the light has been 
affected by a gravitational potential which is implicit in go9. For a 
low gravitational effect, 


| U 
/ Joo ~ 1 “ae 
6.17 


where U is the gravitational potential U = “or Thus goo < 


band: = > rv. Thus the frequency vy = — < —; then 
T 


T 
the original hater at the source is affected by the gravitational 
field, and the observer far from the source sees the light at a lower 
frequency (towards red). The parameter used to characterize this 
phenomena is called z: 
VE — % 


Vo 
where Vg is the frequency of emission at the source and Vv is the 
observed frequency. 

Observe that if T = t; to and T = 71 —7, and assuming that the 
light is emitted at 7, and received at t,, then the distance in classical 
or special relativity between source and observer is c(t; — 71) = 10 
or c(t9 — T) = To. This is not the case if go) # 1, as in general 
relativity. 

The red shift can be affected by more than just gravitation: the 
Doppler effect as well as the expansion of the universe play large 
roles as well. The expansion of the universe refers to how, as time 
goes on, the distance between two points in the universe will change 
as galaxies drift away from each other and as stars within a galaxy 
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travel away from each other. If we assume that the space coordinates 
x’, x”, x change as a(t)x', a(t)x”, a(t)x?, then we may see that the 
red shift is proportional to a. In order to see this, we use the geodesic 
equation for the x° coordinates: 


0 
Co. betes 
ae = =p VU 

da! 
where v4 = —. 


= 
In accordance with the above assumptions, the invariant ds? 
changes as 


ds? = gooc7dt? + 5,,dx'dx* > ds” = gooc?dt? + g;dax'dz' 


Where gj, = —a?(t)d, and 6, = 1 for i = 1,2,3. In other words, the 
metric tensor gy, is 


1 0 0 0 
0 —a?(t) 0 
= 1q 9 a(t) 0 
0 0 0 a(t) 


Then the only Christoffel symbols different from zero are 
Tb = Tig — and Vii _< 
ca c 


(a refers to the derivative of a with respect to t). To prove this, let 
us calculate ['?, and T'jp. 


1 1 
Rai = 59 (Sik + 9k0,1 — 910,k) = 59 Ju10 : 


The only g!* 4 0 is g™ and the only gm, 4 0 is gu. Now, gt! = 
1 . 
~ wt) and gio = nae (t)):= ~ 704. Thus, 
a 


eee (0.9.18) 
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In a similar way, 


1 1 
Be = 59° (Sos + 910,1 — 91,0) = 59 (9110) . 
go = 1.80 
ea (0.9.19) 
Cc 


Therefore, the acceleration of a photon according to the geodesic 
equation is 


A= I}, v'v' = iyi (0.9.20) 
c 
for 7 = 1, 2,3. On the other hand, from the four-velocity U = y(c, v), 
EB 
PES iil = (me, 9) = (—>P) (0.9.21) 


E hp 
For photons, the momentum p = — = —. 


Cc Cc 
2 


E 
The magnitude of P” squared is P“P, = —, — (p)?, but since 
C 
E? 
ie Pe then P?P, = 0, or 


22 ee i 
9; ede dx? q 
ince v’ = = y— and my=™m™, 
dr’ dt Pu 
i da’ i 
Mov’ =m = 
0 ai P 


Thus, (0.9.20) may be written as 


and from (0.9.22) 


— = E?. (0.9.23) 
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Also, 
o dU d HE. 277: 
A _— — Th = — —( — Th i“ 
ap = ql ais ae ae yu) 
For a photon, 
° E dé 
A= —, 0.9.24 
mc dt ( ) 


a _aE 
Ba pay 2 ca 
a dt 
d dE 
lda 1 
ae a ee 
or 
pal 
a 


where @ is a constant. Since for photons F = hv, then as a(t) 
increases (universal expansion) the frequency of the received light 
decreases (shifts towards red). 


4.) Classic Visualization of Black Holes 


Gravitational attraction might be overcome by some high initial 
velocity. Assuming no attraction at a large distance (oo) from the 
central attracting body with mass /, then by conservation of energy 


we have 


mv? GMm 14, 
— = Mv, - 
2 r ae 


The minimum escape velocity corresponds to v., = 0. Thus, at some 


distance r from the center of the body M, a body with velocity vu, 


M 
will escape from M’s gravitational field at vu, > gout Forr=R 


2 
(the Earth’s radius) and M being the mass of the Earth, 


ve = /2|U| = /2(6.17 x 107) = 1.16 x 104 
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Since the velocity of sound is v, = Bayes 
S 


ve = 33.7 Mach 
Us 
For a body with very large M or very small r, v. may reach the 
velocity of light. Places where the escape velocity is greater than 
the velocity of light, so that no energy is released in the form of 
mass or light, are called black holes. 
Assuming that v. = c, then 
us = ae x 19279 
r 2G m- 
Thus, if M is the mass of a galaxy (= 104! kg), r+ 10“ m. If M is 
the mass of the sun (= 10®° kg), r ~ 1 km. If M is the mass of the 
earth (= 1077 kg), r+ 1m. 


Simplified Black Hole Derivations Using Equations (0.9.10) 
and (0.9.13) 
dt dt i 
The idea is to find the ratio or ee = as a function of the 
r v 
main metric tensor goo which we associate with a potential similar 


2 dt 
to the Newton one. Since gop & 1 — oule we will attempt to find A 
r f 


2 
as a function of 1 — aks 


F 
Let us take the equations 
v° goo = Ec (0.9.25) 


and 
C = goo(v’)? + gulv') . (0.9.26) 


Eliminating v° between (0.9.25) and (0.9.26), we get 


C’ goo = £7" + giigoo(v')? - (0.9.27) 
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Assuming 911900, Which is the determinant of |g], is equal to -1 as it is 
in the Lorentz transformation and the Schwarzschild tensor matrix, 
then 

C' go = Ee" — (v")? 


or 


pire 
goo = €* — (=) (0.9.28) 
Thus, 
dt 
oat “dr _ v 
dr dr vi 
dr 


0 
v 

where 7 is the proper time of the light. For —, we can use equations 
v 


(0.9.25) and (0.9.28), so 
v? € 


yi gooV &" — Joo , 


2m 
Now we assume that goo = 1 — —, so 
r 


or 


For r = 2m minus some infinitesimal value €, r = 2m — € or 2m = 
r+e, then 


Orica ay a ce Og 
r r r 
and 9 
m € 
@-(1-—)x@-- 
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Thus, 


62 — ce 1 r | = 2 te 
r *(e a 276° 
r £=0 


Ce ein, ae 
7( area) ws r 
Thus 7 i 
t 1 prdrr lr 
ae i = 0.9.29 
| cJo € & ce€? 2 ( ) 


For r © 2m, € > 0 and t > oo. Thus, a particle at this distance r 
will take an infinite time ¢ to leave the attraction of the black hole. 


5.) The Einstein Relationship in Different Approximations 


The Einstein assumption for space with gravitational forces is to 
constrain the curvature tensor in the form given by the Ricci tensor 
to a value of zero. Thus, 

Ry =0 
or 


Fiyg = a ao oa ee Leg = OS 


om) Vp,o vo” pa vp oa 
Let us now assume that (ER1): 

1. gu is diagonal. 

2. Jyw,0 = 0. 


oe Opps = UtOF t= 2 


200 CONTENTS 

Thus, the only R,, different from zero are Roo, Ri, Roo, and R33. 

Roo = Toot a Pool 1 = Pool h = Vola i Tools 

Ry = Moa Vio4 Tia + Pig to + TisP is + Molto — Pula — Pulte — Tul is 

Ry = T53,9 = Too + VieP 59 + Pool Yo + P3033 — Poe ty — Vel 01 — Poel To — Tael Ss 

R33 = —T 331 = 133.9 +1T3,053 + 133093 — T3301 — TasPn — P330 = 


Subapproximation for the Einstein Relationship 


Let us assume the following (ER2): 
1. Only 2° is relevant. 
Boi ye Pee 


Therefore only goo is considered. Thus, the Christoffel symbols 
which are relevant are 


1 
Ea oa = 9g” 00,1 


2 
and 1 
ile = 59 (—9o0.1) : 
Thus 
: eo ee ee 1 4100 2 
Roo = —Pooa + Polo = 5 (9 900,1),1 = qe 8 (900,1) 
and 


1 1 
Ry = To, +T%ol'to = 5 (9"" 900,1).1 of: z(9"" 900.1)” . 


Roo = 0 means 


1 saul 


Roo = 900,11 — 59 goo)? =0. 


Ri, = 0 means 


1 
Ri = 91’ Jo0,1 + 9°°go0,11 + 5 (9) (G00.1)° =0. 
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hige ee Roog™ + Ru, then 


S = 29°°g00,11 + 9'1'900,1 = 0. 


Now, 
g® = at. fet HOA 
ot Joo (goo)? 
so S becomes 
2 
co 2g ona = (goo.1) 9 
(goo) 
Thus, 
S=0 if Joo,11 = 0 (0.9.30) 


which is a type of Laplace equation for the potential goo. 
Equation (0.9.30) might be put into scalar form by extending the 
subscript 1 to pv: 


Joo,11 — Jood,uv - 


Let us introduce a generalized Laplacian UO, =, wv, so that we 


might write the scalar form of g” goo, as 


ge pv Goo = On: (0.9.31) 


Simplifications for the Ricci Tensor and the Gravitational 
Einstein Equations in a New Approximation 


Now let us make the following new assumptions (ER5): 


1. Gu is diagonal. 
1 Sih i 1 
2. goo = ——, which is a function of x* only. 
911 


3. go2 = 933 = —1. 
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According to the simplifications of (ER1), the only relevant Christof- 
fel symbols are T3,, Pg), and ['j,. Thus, 

Roo = Too, + Dolo: = DoF i 

Riu = Toi sb TYol'to a Plo, 


Ro = 0 
R33 = 0 
Now, 
1 
ie = 59 900.1 
0 1 © 
Poy =, 94 Jo0,1 
1 
Pe = 59 HA 
sO 
1 1 11 1 
Rn as (3am) (as) + (H-) 
00 5 (9 ( Joo,1)) 1 5( Joo) 900,1 ogo Te BA Joo) 900,1 
1 
— ~ 5 900900,1,1 , 


Thus, Rj, = 0 means that gooi1 = 0. Roo = 0 and gi; = 0 mean 
the same thing. 
If goo = y(x*), then y” = 0, or 


TN aan 


which is the Laplace equation for y. Therefore, y could be any 
function of x1 (let us call it r) such that Ay = 0. So goo may be a 


combination of a constant and —, such as the following: 
is 


B 
goo = A+ — 
: 


with A and B being constants. 
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0.9.4 Calculation of the Christoffel Symbols Un- 
der Some Simplifications 


Assume 
1. A diagonal matrix. 


2. Static conditions: g,,,9 = 0. 


1 
Then, since gy = <7 and gi, = 0 for p# v, 
g Vv 


a aa aa i 
Da =Jg a =9 5 (Guaw + Jav,p Quien) . 
1. For a = 0: 
00 if 
5 (9u0.v + Jovy) : 


Then, only = 0 and/or v = 0 will make T° 4 0. Thus, only 
Teo. Vb1> Po, 03; may be different from zero: 


Ot a 
Dw 9 


1 1 
Ea = 59 (Goo.0 + goo,0) = 0 le = 59° 900.1 
1 1 
is = = 9g" Go0,2 ie = = 9" goo,3 
2 2 
2. Fora = 1: 


1 
ine = 5 (Gu - Div. Guyi) : 


The first two terms are zero if uz # 1 or v 1, and the last 
term is zero if v € yw. The values of jz or v are 0, 1, 2, or 3. The 
couples (u,v) are (0,0), (0,1), (0,2), (0,3), (1,1), (1,2), (1,3), 
(2,2),(2,3) (333). Lins: 


1 
To = 59 (~900.1) P= 0 Pe=v 
1 1 
ag = 0 Dy = 59 ua Es 7 59 gu2 
1 1 
ie = 59 913 in = 59 (—922,) ee =0 


1 
Le = 59 (9384) 
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3. For a = 2: 
= 507 (ure iid = Gane) 
Thus 
Too = 59°*(~90.2) PO 420 
T%3 =0 i, = 597-9112) Re = 59729201 
Ts =U Tp a 59°9002 T53 = 5979009 
[33 = 50° (—I552) 
4. Fora =3: 
Ts, = 9 (use + Osi = Gans) 
Thus 
Too = °° —g003) Pe 0 Ey 6 
To3 = 50° 9s30 Ty = 50°(—g3) Th = 0 
l3 = 500551 T. = 50° —g223) [33 = 50832 


133 = 590585 
In the previous expressions of the Christoffel symbols, we assume 
that each matrix element g,,,, is a function of x’, x*, and x3. For the 
rest of this lecture we will assume that goo and gi, depend only on 
x', and that gg. and g33 may depend on x! and x? but not on 2°. 


With this in mind, the only [: which are nonzero are: 


1 1 1 

bie 59 900.1 Eos 59 (—Io0.1) Lo 59 9a 
1 iI 1 

ies = 59 (9221) P35 = 59 (—9384) ee = 59 924 
1 1 1 

lies = 59 (933.2) Ee = 59 9331 cs. i 59 933.2 
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The Schwarzchild Metric Tensor 


The coordinates used are spherical ones (x1 = r, x? = 6, x? = W). 
Thus the space-time invariant ds? has the following form: 


ds? = c” goodt? = Grr(dr)? — r?(dQ)? 


where (dQ)? is the differential solid angle, (dQ)? = (d@)?+(sin 6)?(dv)? [x 
Schwarzchild took goo and g; as functions of r only: 


goo = eZ” go = gut 
ee git = —e 
1 
2 2 
922 = —T GQ F=7-o 
: 
933 = —r? sin” 6 i : 
r2 sin? 6 


where v and A are smooth functions of r. In other words, the matrix 


is 


et 0 0 0 
se 7 0 
0 0 —r? 0 
0 0 0 —r? sin? 6 
Also, 
goo,1 = 2v'e” 911 = —2X'e 
922,1 = —2r 933,1 = —2r sin? @ 


933,2 = —2r? sin 6 cos 6 
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Then the Christoffel symbols different from zero are 


1 1 
Ea = = 9 G001 = 6° (Oye?) = py! 


2 2, 
1 1 
Poo = 59" 900.1 = —5(-e™)(2v'e™) = p/e2v-r) 
1 1 
Lie 59 91 = 5 (-e)(—2Xe"") =)’ 
1 1 
Poo = 59" 922.1 = Spe eh) = —re 
1 1 1 1 
Fr == e SS Sl SS —?2 = — 
12 549 922,1 a “a) r) a 
1 1 
T2, = —~9"*9339 = ——(——-)(—2r? sin 8 cos 9) = — sin @ cos 0 
ee 2 yaa pe: 
1 1 
a= 59 933.1 = Seine at sin? 0) = —re~* sin? 6 
1 1 1 1 
3 _ +33 = Soa 
TO 90 Ol 3 ene 
1 1 1 . 
33 = 59 933.2 — pl mperrecy yas sin 0 cos @) = cot 6 


Now, let us write the Einstein relationship using the Schwarzchild 


matrix: 


y) / 
Rag =. (eg = py st yr =) oY )e2e) =i () 
ifs 


Rea 2 NS 0 
Rigg = -1 +e 4 (Lr = 7X) S0 
R33 _ Roo sin? 6 


For Roo = 0 and Ry, = 0, we get 


or 


y+ =0 


vy+.A= const. 
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Assuming a flat space for r > oo, the constant must be zero. Thus, 
yv=-—. 
Now, eliminating A in Ra: = 0, we get 
—1+e"(1+2rv') =0 
or 
—(re*”) = 1 


or 


re’ = r + const. 


Let us call the constant —2m. Thus, 


Therefore, 


The same form that we obtain with the coordinates x° and x!, —— 


. 
is the potential due to gravity which deviates goo from 1 (flat space). 


Now, 


1 1 
2X —2v 
goo ie 2m 
Ts 


Let us show now that this solution which establishes that v’ = 
—X' and e?” = 1 — — satisfies the Einstein relationship Ro = 0. 


: 
For that, we must have 


+N + <0 


or 
9 / 
Pipe S= 6, (0.9.33) 
r 


From (0.9.32), 
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or 


F m 


ae r2 — 2mr 
and 
he< 2m(r — m) 
(r2 —2mr)? © 


Substituting the above into (0.9.33), 


—2m(r — m) 2m? 


2m 
| | 
(r2 — 2mr)? (r2 — 2mr)? (r2 —2mr)r ~ 


Multiplying and dividing the last term by (r? — 2mr), we find that 


the above is equal to zero. Thus, Rop = 0, and e?” = 1 — ll makes 
all Ricci tensors equal to zero. This is condition fniosed onto the 
fourth dimensional curvature by Einstein. 

With the metric tensor 


2m 
1 — — 0 0 0 
r 
1 
0 [oa 0 0 
0 0 —r? 0 
0 0 0 —r? sin? 6 


the fourth dimension invariant length, after using the Einstein rela- 
tionship (i.e., the Ricci tensor equal to zero), becomes 


2 


cao 
‘a 


2m 


ds? = (1 — —)ed¢? - — 77d6? — r? sin? @ dy? . 
- 


By using the variables x? and x°, such that dx? = rd@ and dx? = 
rsiné dw, then 


0.9. LECTURE 9: FUNDAMENTALS OF TODAY’S COSMOLOGY20). =! 


The metric tensor then becomes 


2m 

1 — — 0 0 0 
Me 1 
0 eT 0 0 
0 0 —1 0O 
0 0 QO -l 


and the determinant |g| = —1, just like in the Lorentz matrix. 
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0.9.5 The Hubble Law 


From observed red shifts, in 1928 Hubble found that the galaxies 
recede at a velocity proportional to the distance, 


Noe ae sk ge (0.9.34) 


Let us assume that this law is valid for all times t, or 


or ar(d) 
r(t 
= H(t)r(t 
ST ORO 
Integrating, 
t 
In = / H(t) dt 
TO to 
or ; 
Eee (0.9.35) 
TO 
=o Why 8 = oe 
The ratio is called R(t). As time passes, the galaxy’s dis- 


ro(to) 


tance increase by this factor R(t), which is related to the Hubble 
parameter H by 


or 


H(t) = —[In(R(2))) . (0.9.36) 


The function H(t) may take on different forms depending on what 
parameters are used to define the universe. Next, we will discuss 
particular cases of H as a time-dependent parameter. 


1. H(t) is constant and independent of time, so 


u(t) = = H(t) (0.9.37) 
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or by integrating, 
r it 
In—=H | dt =H(t—to) 
TO to 


or 
rO=nigne (0.9.38) 


The acceleration of the galaxy will be 


dr (t) 
-— H 
dt dt 


= Hv(t). (0.9.39) 


or 


or 


Therefore, 
a(t) = Huo(to)e#* . (0.9.40) 


Also, from (0.9.39), 
a(t) = Hv(t) = H(Ar(t)) = H?r(t) . 


Thus, 
a(t) = H?r(t) . (0.9.41) 


A force F(t) on the unit mass is 
F(t) = a(t) = H’r(t) 


or a potential U such that F = —V,U is 
2 


U(t) = =i . (0.9.42) 


2. Let us assume that 


212 


d 
The Hubble Law says that v = rH, but v = = so 


dt 
dr 
t)=—=r(t 
ot) = T= r(t) 
Integrating, 
r t 
In — = In — 
ro to 
wr r t 
R Ss 
ro to 
Thus, 


v(t) _ r(t)to a 
vo(to) ro(to)t 


CONTENTS 


(0.9.43) 


which means that v = vo, or v is constant and independent of 

time. In this case, the acceleration a = 0 and there are no 

net forces acting in the galaxy. We observe that as t > 0, 

H — oo, which corresponds to the initial "big bang” where all 

the distances r = rp were zero. Also, as t = to, R = 1. Since 

R= = the number N(t) density which is proportioanl to r? 
0 


are proportional to t?. That is, 


. Let us assume that 


The Hubble Law says 
r(t) 


u(t) =r(t)H = = 


dr(t) 


Since v(t) = ——, then by integrating we get 


dt 
B 
In a) =n(-) : 
to 


(0.9.44) 


(0.9.45) 


(0.9.46) 
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Thus, 


ie oe ee Ge (0.9.47) 


From (0.9.48) 


v(t) ( t < 
vo(to) to 
For 8 = 1, we have the same result as before: v is a constant. 


For 6 #1, 


U =o B-2 
ie he 7 a (2) Jat = (B- 1) gato) 
(0.9.50) 


Let us substitute vp by ro using equation (0.9.49): 


ee cc ee _ 4yTo(to) \ 62 
aoe = (4(3 5 )i (0.9.51) 
Again for 8=la=0. 


Let us get the acceleration (which is the force per unit mass) 
as a function of r(t), assuming rp and to as given values. Since 


or simply 


Thus 


214 


CONTENTS 


and 
B-2 


p-2 r\z, 12 r\ 8 ,g-9 

2 —|(~)"] = (=) af. (0.9.52) 
TO TO 

Using equation (0.9.52), the acceleration becomes 


B-2 2 


Tr TN OTB A= = B-2 
a= (6 -1)3(=) 7 8 = (8 1) 2rer 
to \To 
In other words, 
a= er (0.9.53) 
where ; 
€ = B(8-Ir§tp? . (0.9.54) 


2 
For 8 = ~, we have a force inversely proportional to r?, which 


is a Newtonian type of force. 
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0.9.6 Expansion of the Universe Under the Fried-3qggia 
mann Fundamental Metric and a Tensor 
Force Acting on the Material of the Galaxy 


In a simple fluid dynamic, classically the balance of matter in a 
given element of volume dV is given by the continuity equation 


~+4+V-7=0 (0.9.55) 


where 7 is the current density associated with the velocity by 7 = 
pu. If the element of volume suffers stresses from the rest of the 
continuous media, then the equation much change to include those 
internal forced plus the possibility of external forces. 

In four dimensions, we define the energy-current density J“. Ob- 
serve that we have added the word energy” to the current density, 
because after the equivalence between mass and energy we want to 
include those terms which classically may be studied only as energy 
contributions different form the mass-related terms. By definition, 


=> 


J" = poU" = poy(c, t) = (pe, J) (0.9.56) 


where 
j=pi and p=poy. (0.9.57) 
y accounts for the length contraction when going from system S' to 
the Sp or proper system (see Lecture 7). 
The new equation of continuity in four dimensions is the follow- 
ing: 


JH, =0 (0.9.58) 


’ 


where the coordinate variables are: vy = 0, 2° =ct;v=1,2'=72; 


PRI ayy S36 Se: 

This equation of continuity is not complete because our test el- 
ement of volume may be affected by stresses acting on the element 
of volume dV through the walls of the element dV. 
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In analogy to the classic studies of continuous media, we assume 
that the following tensor describes those actions in the rest frame 
So: 

0 0 0 O 


cal) 
0 tH ty tf 
0 
0 


py 

Tp = ¢21 $22 $23 
0 0 0 

31 32 33 
to t t 


The tensor T, may be understood in the following way. Assume a 
surface enclosing a differential unit of cubic volume dV of the fluid 
within the three orthogonal axes (x,y,z). The stress to the fluid 
on each of its sides is described by the two index tensor t? where j 
stands for the axis which is normal to the face of the cube and 7 the 
component of the stress on that surface in the direction of the axis 
i. The unit vectors are €1, €2, €3. 


> 


€3 


€1 


Let us apply the Gauss divergence theorem to the cube. The 
force on the face ds; perpendicular to é; is t! ds;. The flux on 
the forces with direction of é; is ft’? ds;, which is equal to f ae dV 
where j = 1, 2,3. Thus, by the theorem we identify the force fi due 
to the stress as the divergence of the tensor t’’: 


oy 
fi=t; . 


Also, the tensor expressed in the rest frame as i, may be trans- 
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formed to the laboratory frame as 
tii tl2 tis 
¢21 $22 $23 
p38 $32 $38 


where t” = tj, tj, = toy, and tf; = of The expression for ¥y is 
the following: 
1 
yu) = 5 
u 
‘2 


where wu is the velocity of the fluid in the direction of €;. The 
above transformation may be derived from the fact (see Introduction 
to Special Relativity by Wolfgang Rindler) that the force a on the 
element of area dsj, f= t? ds,, is equal in both systems S and So 
when the velocity u is parallel to f and the mass mp measured in So 
remains constant. Also, we must take into account that when ds; 
has a side parallel to €; (or ui) then as it moves ds, will suffer a 
contraction. For instance, dsy = Ax - Az, so dsg2 = dsgy(u). Thus, 
dsot!? = dsoati? = tj?y(u)dse, and therefore t!? = t}?y(u). 

For the other two faces perpendicular to €; we have ds = dso, 
as the motion of them will not have any relativistic effect. For 
example, the force on the face perpendicular to €; in the direction 


any wv yu!) 
f t4ds, = y(u)(1 — —)t24dsq = t24dso,. In general, 
of €3 is tds: = y(u)(1 — )to dso ve oe g 
/ 
(ul) = y(u)(1- = In our case, u’ (the velocity with respect to 
yu) c 
ie 
the rest frame) is zero, so y(u’) = 1 and ¢?! = (u)’ since ds; = dso 
y(u 


(no length contraction). 


For the diagonal element of the stress tensor like ¢! 


we have 
no variation in the area nor in the force. In fact, according to 
the transformation of the forces when the force is parallel to w, 
t)'dso, = t''ds,;. But as the surface moves with the fluid, dso, = dso 
(no contraction). Thus, tj’ = ¢1!. 
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By changing to the frame S, respect to which So moves with 
velocity wu, this tensor should include the energy due to motion. We 
may write it as 


pw | OY PO Pp” Pp 


The energy density is pc”, the total momentum density or energy 
current is cg’, and the total stress on the unit area as seen from the 
laboratory system is p’. 

Let us see next how g' and P” are related to the original tensor 
t. Let us assume that f’ is the external force on dV which is 
contrasted by the stress t’’ which creates a force in the 7 direction 
equal to t' per unit volume. Thus, the equation of motion is 

ie i i 
FAC dV) = (f'—t)dv . (0.9.59) 
The first side the equation implies the variation of g’, as well as the 
variation of dV. 

In order to see how t is part of g’ let us imagine a surface 
through which the fluid moves across with velocity uw. The energy 
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transmitted across in the direction of é; is t/’é; -d = t/'u;, where 
repeated indexes implies summation. This power is transformed 


into momentum density by dividing by c* in accordance with the 
ji 


relationship m = —. So the momentum density is = . As we 
have said before, the total momentum density will then be 
; pst 
g =pu + Uy (0.9.60) 


If heat or other types of energy had been crossing the imaginary 
unit plane, they should have been added to the previous expression 
of the momentum density. As we will see, all this extra energy will 
never have shown up in the equation if it were not for the equiva- 
lence between energy and mass. Observe that this extra energy is 
crossing the area instantaneously, so that the motion of the fluid is 
unimportant. The motion is taken into account in the first term, 


) 


pur 


Returning to the equation of motion (0.9.59), we observe that the 
element of volume will change with time. In fact, in the laboratory 
frame with coordinates (ct, x', x”, x”), the face a (see figure) in an 
interval dt sweeps a volume proportional to u'dx?dx%dt, while the 


aii 
face 3 sweeps a volume of (u' + at da") dada. 
xc 
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Thus, the change in volume in dt is 


dV du 
— = ; 0.9.61 
dt Oa! ( ) 
So equation (0.9.59) becomes 
d. Ous — 
~ gi + gi—_av = (fi—t%)dV . 9.62 
dV 9 + ga GW (f' — t)dV (0.9.62) 
The first term may be explicitly written as 
d , 0g'  Ogidzi dg' dg’ , 
—— ; = -Uu . 0.9.63 
dio ~ Ot * Oni dt Ot Oxi” ee) 
Then (0.9.62) becomes 
Ogi Ogi 5 OW aay 
— pa ft _ 44 
att Gad 19 aga Pts 
or as P 
g ind 4 489 i 
=~ eo ae 0.9.64 
fy gin +19) = J (0.9.64) 


The four-dimensional current which is given by 


=>, 


J = poU = poy(c,u) = (ep, pt) = (cp, 7) 
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has a divergence of zero. That is, 
Jt = 0. (0.9.65) 


This expresses the continuity of matter in four dimensions, similar 
to the classic continuity equation (0.9.55). 
The four-dimensional force is defined as 


FE = mo = 14 mol) = 1-(mole,) = (7, F) 
(0.9.66) 
In the case of a pure force (ame =O 
OF = (Fi 
and so aed 
Fao 7) (0.9.67) 


When an external force is acting on the element of volume of unit 
value, then the force Fi’ is related to the tensor T”” by the equation 
of continuity 

Ty = Fy. (0.9.68) 


For » = 0, we have T” = (c’*p,cg', cg’, cg?). For up = 1 to 3, we 
have T’” = (cg’, p"', p”, p’). Thus, there are two types of equations 
as a result of (0.9.68): 


2 
Le eT (0.9.69) 
Ot 
and Ais bass 
ig. Ope. 
at Gar at (0.9.70) 


The equation (0.9.69) implies the balance of energy in the unit 
volume of the fluid. The first term expresses the energy rate inside 
the volume, while the divergence of c?g taken to the second side of 
the equation is the energy current density that is received by the 
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fluid through the surface of the element of volume. Finally, ( fe u) 
is the power delivered by the external force of f into the element of 


volume. 
Comparing (0.9.70) and (0.9.64), we see that 


PLSgwy te . (0.9.71) 
By using (0.9.60), P’ is expressed explicitly as a function of t¥. 


A Particular Case of the Tensor Force JT” 


Let us assume that, just like in an ideal gas, 


pe -0i- 20 20 
0 -P 0 0 
TH = 
0 0 -P 0 
0 0 0 —P 


It is a diagonal matrix with p dependent on time, but P is constant 
and isotropic. We would like to work with the matrix T7”” and the 


Friedmann metric tensor gy. 


1 0 O 0 10 0 0 

rt @  Q ie ol inte bie 4 
Tw Tig | 09 —P 0 0 e a = ar tp 
: 0 0 -P 0 0 0 -+ 0 00 5 
0 0 0 -P - 3fl : 

0 0 0 -5 O:Q: 1: Se 

a a? 


In Lecture 4 we showed that the Ricci tensor has an absolute (or 
covariant) derivative equal to zero. We also showed that the Einstein 
equation under the action of a tensor TJ” is 


1 
RM = — gl” R= xT (0.9.72) 


where R#” is the Ricci tensor, and since the absolute derivative of 
the first side of the equation is zero, then 


Dol 
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This conservative equation can be written as 


oTw 


DES Ox oF eae 1 ie =0. 
Then for \ = p it becomes 
Oly Vv OL Va 
el eal ae Os: (0.9.73) 


Ox! 


The following table summarizes the terms different from zero in the 
previous equation, given that the only non-zero Christoffel symbols 
for the Friedmann matrix are 


G00 fe pa (0.9.74) 
Cc ca 
where 7 = 1, 2,3. 
vy |p| A= | a] Non-Zero Term 
0 | 0 0 Ty =cp 
OQ % 4 EN = 4p 
Oras), Ses 0 ere Fae 


Summing all of these contributors for v = 0 and z = 1, 2,3, we have 
the following: 


Bag? pao =p 43-0: 
ca a a 
Particular cases: 


1. For P = 0, we have 


p_ 34 
p a 
Integrating, 
[mn ply = —B[n aly. 


Assuming that the limits for p and a for t = 0 are p(0) = 1 and 
a(0) = 1, we have 


p=. (0.9.75) 
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2. For P=? =, 
2 a4» 
3-(= = 
cp t+ 3 (sec) 
and : 4 
eee cee 
p a 
Integrating, 
In p = In(a~*) 
Thus, 
1 
p= at 


CONTENTS 


(0.9.76) 


with the same assumptions of the limit values of p and a for 


t= p(0) Sand a0), 


